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ABSTRACT 



Quantum Chromodynamics predicts radiative terms which, at low transverse momen- 
tum Qt of an observed multi-particle system, can be large enough to threaten the 
legitimacy of a fixed-order perturbative treatment. We present a subtraction-style 
calculation which sums these terms to all orders for a direct photon plus jet final 
state, and show that the kinematic effect of such a resummation is significant 
even after integration over the jet. An improved matching algorithm is prescribed 
for the transition to the high-Qj- perturbative regime, and re-examination of the 
nonperturbative regime indicates that simple Gaussian parametrizations are to be 
favored, with coefficients that depend not on the photon-jet mass Q, but on the 
color structure of the initiated subprocess. Agreement with data is significantly 
improved, and the method is expected to have applicability in general for processes 
with steeply-falling spectra. 
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CHAPTER 1 



INTRODUCTION 

1.1 QCD as a Quantum Field Theory 

In particle physics we are continually searching for the most fundamental de- 
scription of matter and the forces which govern its behavior. Those who wonder 
about the properties of macroscopic objects may now look to molecular and atomic 
physics; those who inquire about the structure of the atom may now be answered 
with the details of electrons and nuclei. In each case a multitude of types of matter 
has been replaced by a more fundamental set, the variety of types being replaced 
by a variety of ways these types can combine to form bound states. Concurrently, a 
unified description has arisen for many forces once thought disparate. 

The situation is no different for the atomic nucleus or its inhabitants. In addition 
to radiation observed naturally, data on nuclear structure came from experimental 
probes of the nucleus, usually involving beams of particles known to be themselves 
structureless. As more energy was made available, the same method was used to 
probe individual nucleons (protons and neutrons). The mid-20th century saw a 
proliferation of new strongly-interacting particles (collectively known as hadrons) as 
the energy barrier for their discovery was crossed, and again it became natural to 
look for common properties ("good quantum numbers"), clues to a symmetry on a 
deeper level. 

In 1963, motivated by the observation that the known hadronic families seemed 
to exhibit symmetries of a particular branch of mathematics known as Lie groups, 
Gell-Mann [1] and Zweig [2] proposed a model of hadronic structure based on a 
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limited assortment of sub-nucleonic constituents called quarks. Each quark came in a 
particular flavor, had fractional electric charge, and had a corresponding antimatter 
partner of opposite charge. Together, they were able to be combined in ways that 
reproduced the known subnuclear spectrum. 

However, there was a spin-statistics problem. Systems of fermions (particles 
of half-integer spin, like the quark) must have wavefunctions which are totally 
antisymmetric under interchange of any two particles, while total symmetry is 
required of boson (integer-spin particle) collections. The quark model posited that all 
hadrons were quark (fermion) combinations, those of half-integral spin (the baryons, 
such as the proton and neutron), being composed of three quarks, and those of 
integral spin (i.e. mesons) composed of a quark and antiquark. Unfortunately, some 
of these combinations required symmetry in spin, space, and flavor, leading to a 
totally symmetric wavefunction. If the above compositions of baryons and mesons 
were correct, a new degree of freedom was required in which to properly perform the 
antisymmetrization [3] . 

A new type of charge, posited by Greenberg [4] and which eventually became 
known as color [5], thus solved the statistics problem. Again, group theory was 
instrumental (in this case a color SU{3) fundamental representation). With three 
colors, hadrons could be properly antisymmetrized if each quark in a baryon was 
a different color, but those in mesons were of the same color. The quark model 
seemed to be rescued, despite the fact that no individual quarks, or any other 
"colored" states, had been seen in the laboratory. The lack of experimental evidence 
for anything but colorless states could be postulated as a rule {color confinement), 
but without an understanding of the forces present among quarks, such rules would 
remain baseless. 

Classical theory described interaction at a distance in terms of fields, with 
forces that depended on the charges of the particles involved and the distance 
between them. Modern quantum theory, although maintaining use of the term 
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"field" , had replaced the concept with the (macroscopically equivalent) exchange of 
quanta, specific to the type of charge involved, which carried energy and momentum 
across the intervening distance. One theory of this type. Quantum Electrodynamics 
(QED), had already been successful in remodeling the electromagnetic interaction 
between charged particles (as well as all the known properties of light) in terms of 
intermediate, long-range bosons called photons. 

It had also long been known that symmetries of field theories gave rise to 
conserved quantities, like charge and momentum. In particular, conservation of 
electric charge in the coupled system of photons and charged particles was shown to 
be a consequence of the symmetry of the theory with respect to certain local gauge 
transformations, performed simultaneously on all fields. As the quarks had not only 
electric but color charge, both of which were expected to be conserved, it was natural 
to look for gauge bosons of the color force. However, unlike the electromagnetic 
interaction, this force was known to become weaker at higher momentum transfer; 
that is, the individual quarks behave more like free particles [6-7]. This property, 
asymptotic freedom, was searched for in the known quantum field theories, and 
found eventually to be associated exclusively with those gauge theories in which 
the generators of the group symmetry do not commute; i.e. are non-Abelian [8-11]. 

Fritzsch and Gell-Mann proposed, then, that the group symmetry of a non- 
Abelian gauge field be associated with color symmetry, and Quantum Chromody- 
namics (QCD) was born [12-13]. The gauge bosons of this theory, called gluons, have 
eight members, corresponding to the number of elements in the adjoint representation 
of the SU{3) color group. They couple to the color charge of the quarks, binding 
them together within the hadron. Being non-Abelian in nature, they also carry the 
color charge themselves and thereby couple with one another (unlike the photons 
of QED, which are Abelian). Not only does this explain the strengthening of the 
strong force with distance, and consequently why we do not see any free quarks or 
gluons (collectively known as partons), but also has implications for the question 
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of color confinement. ^ Indeed, Wilson [15] has shown that, for sufficiently strong 
coupling, the only finite-energy asymptotic states of QCD are colorless. Any attempt 
to separate one parton from another results only in the creation of a quark/antiquark 
pair somewhere in the intervening distance, the resulting colorless states propagating 
outward as everyday hadrons. This is the process of "fragmentation" or hadroniza- 
tion. 

1.2 The Parton Model 

Asymptotic freedom brings with it benefits and constraints. The collision of 
two hadrons is an enormously complex problem, and cannot be solved exactly, 
as the component partons can interact in infinite combinations. First, as we 
are not able to directly observe individual partons, we do not apriori know the 
precise distribution of momentum among the quarks inside the nucleon, nor can we 
predict from first principles the number of each type represented at a given energy. 
Just as an electron is at all times surrounded by its own cloud of photons and 
electron-positron pairs produced by these photons, which tend to alter the effective 
charge seen by a neighboring charge, so too do the quarks constantly emit and 
reabsorb gluons, themselves producing a cloud of quark-antiquark pairs. Increasing 
the interaction energy only serves to probe deeper into the cloud, thus changing the 
effective distribution [16]. Thus we find that the best description we can achieve 
of the nucleon at high energy is, probabilistically, a set of continuous densities, 
one for each flavor/ ant ifiavor, plus a gluon density. These parton distribution 
functions give the probability of finding a certain type of parton (quark flavor or 
gluon) at a certain momentum fraction and interaction energy within the nucleon. 
Although prediction of these densities is ultimately a non-perturbative endeavor, we 
can perturbatively derive the relationship between densities at different interaction 

^The term "parton", coined by R.P. Feynman [14] before the advent of QCD (and in a purely 
dynamical context), has since been identified with its basic building blocks. 
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energies. Measurement of the densities at one energy thus provides predictions for 
the rest. 

Secondly, if the final state of interest contains particular hadrons that are to 
come from partons involved in the scattering process, instead of attempting to 
describe quantitatively the mechanism by which this hadronization occurs, we can 
currently account for it only phenomenologically, by introducing a fragmentation 
function for each of these partons. In a manner analogous to the parton distribution 
functions, these give the probability that a certain hadron will be produced from 
the fragmenting parton, with a specific fraction of that parent's momentum. Both 
distribution functions and fragmentation functions must ultimately be measured by 
experiment. 

The good news of asymptotic freedom, however, is that at small distances 
(and correspondingly high collision energies), the strong force decreases and the 
partons behave more like free particles [17]. This allows us to factorize what we 
can perturbatively describe (the interaction of free partons) from what we can't 
(distributions and fragmentation). In the parton model of high-energy hadronic 
collisions, the initial reaction takes place over such a small time scale that it becomes 
plausible to treat each hadron as a collection of free partons, each traveling in the 
same direction as the parent hadron, and only one of which takes place in any 
given "hard scattering". That is, one may calculate basic subprocesses involving 
one parton from each hadron, then sum incoherently over the contributions of 
all partons via convolution with the parton distribution functions. Similarly, the 
process of fragmentation is assumed to take place over a time scale much longer 
than the hard scattering, and can therefore also be treated independently. Thus a 
description of hadronic collisions is necessarily built up from both perturbative and 
non-perturbative pieces. 

One of the conveniences the new quantum model allows is the schematic represen- 
tation of a subprocess as a Feynman diagram, each fermion or boson being represented 
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by a line. Line intersections (vertices) can represent instances of absorption, emission, 
or pair production, depending on the direction of time. To each element of this picture 
a mathematical expression is assigned, in such a way that one can easily derive the 
correct expression for a subprocess: one simply draws its diagram, starts at one end, 
and follows the lines and vertices, writing down the corresponding factors in order. 
The result is called the amplitude for the reaction, and its magnitude squared is a 
measure of the probability that the incoming particles will interact in precisely this 
way to form the outgoing particles. Since it is the vertices at which the force is 
felt, the overall strength of the reaction and its probability for occurring become 
dependent upon the number of vertices in the diagram. 

Of course, for a given set of external lines, there is an infinite number of diagrams 
that could be built. Each makes a contribution to the total probability for the 
reaction, in relation to its order, or exponent of the coupling strength. It becomes 
clear that if we hope to make meaningful predictions through this perturbation theory, 
this relation must be inverse, i.e. the coupling constant (unit of charge) used must 
be less than one. Otherwise, the diagrams with the largest contribution would be the 
ones with the most vertices, and we cannot analyze an infinite number of diagrams. 
With a small coupling constant, we can usually stop at first or second order, confident 
that the higher-order diagrams are of negligible contribution (an important exception 
will be discussed in Section 1.4). 

1.3 The Perturbative Domain 

Having a small coupling constant is not necessarily enough; if the coefficients of 
successive orders grow, we may lose the applicability of perturbation theory. Hence, 
while in many cases we are justified in calculating only to leading order (LO), to 
confirm the convergence of the series and provide corrections to the leading-order 
result, we must often look at the next-to- leading order (NLO) terms as well [18]. 
Loop diagrams, having the same external legs but more internal vertices, will of 
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course contribute, as will the release of bremsstrahlung radiation which inherently 
accompanies acceleration of charged particles. Thus inclusive data is considered 
(all events are kept which contain the desired final state particles, regardless of the 
detection of additional debris), and all contributing diagrams up to a given order are 
calculated. 

As an example, we'll show the leading order (LO) and next-to-leading order 
(NLO) contributions to the production of a muon pair in hadronic collisions. At 
lowest order, there is but one contributing subprocess, proposed initially by Drell 
and Yan [19-20]. A quark from one hadron annihilates with an antiquark from 
the other to produce a vector boson (hereinafter chosen to be a photon), which 
subsequently decays to a muon-antimuon pair. The expansion so far consists of one 
term of second order in the electromagnetic coupling a, zero*'* order in the strong 
coupling Os (see fig. 1.1). Order a'^Og subprocesses include gluon bremsstrahlung 
from one of the incoming quarks (fig. 1.2), and quark-gluon scattering (fig. 1.3), 
which produces an additional quark in the final state. The exchange of a gluon 
between incoming quarks is a contribution of order a'^a^ (fig. 1.4); its interference 
with the Born term, however, leads to an a] quantity. Photon bremsstrahlung or 
exchange is also a possibility, but these order processes are suppressed by the 
small size of a/a^. ^ 

As a rule, bremsstrahlung and loop terms contain various mathematical singular- 
ities as masses and momenta are taken to certain limits, but these divergences can all 
be removed consistently in the full theory, as they must to avoid the absurdity of an 
infinite cross section [21-24]. Infrared singularities arise from both bremsstrahlung 
and loop diagrams, in the limit of zero radiative energy, and cancel if these cou- 
sin the diagrams, solid lines denote quarks, wavy lines photons, and curly lines gluons. The 
coupling constants e (electromagnetic) and g (strong) are more commonly traded for the quantities 
a = e^/47r, = 5^/47r. 
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Figure 1.1. Lowest Order (Born) Drell-Yan diagram, lAif ~ [e^]^ ~ a^. 




Figure 1.2. Gluon Bremsstrahlung diagram, lAif ~ [e^f?]^ ~ a^a, 

8 



Figure 1.3. Compton diagram, ~ \^9^ ~ c?ol, 




Figure 1.4. Virtual Drell-Yan diagram, \M.f ~ ~ a^a^. 
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tributions are both included. Loop diagrams also contain ultraviolet divergences at 
infinite loop momentum, and these are factored into a redefinition, or renormalization 
of the (unobservable) bare masses and couplings. The remaining class of divergence 
arises in the limit that bremsstrahlung radiation becomes parallel with the emitting 
particle. These collinear singularities are factored into a redefinition of the parton 
distribution functions (or PDFs, the probabilities that the colliding hadrons contain 
the relevant partons). In both renormalization and collinear factorization, the 
singularities (or poles) are accompanied at order a" by a power series of logarithms 
involving an energy characteristic of the subprocess, e.g. the invariant mass Q of the 
muon-pair. In the leading log approximation, the largest of these logs at each order 
is summed over all orders and factored into the couplings and PDFs also, making 
them energy-dependent. ^ As it has been shown that (for each parton species) the 
divergent logs arise in the same way for all processes with a large scale [25-27], once 
the new coupling and distributions are calculated, they are transportable in this 
approximation to other reactions. A lot of work is saved in the rough calculation of 
a new process, as only the Born subprocess need be derived, and the most dangerous 
higher-order terms accounted for by convoluting the Born term with the new coupling 
and parton distributions [28-29] . Further improvement comes, of course, with the 
inclusion of subleading higher-order terms, but as we shall see, even this may not be 
sufficient for processes which also involve a small scale. 

1.4 The Need for Resummation 

During the late 1970's, the leading log approximation was applied to a number 

of processes [30], although we will focus here on muon pair production. In addition 

■^These energies are historically called scales (the collinear in particular being the factorization 
scale, which we'll denote as Mf). One of the assumptions of the early, or "naive" , parton model was 
that the parton distributions should "scale with" the momentum of the parent hadron; i.e. should 
depend only on the momentum fraction involved. As we are here taking our picture of the hadron 
as a hierarchy (or cloud) of partons which we sample at an effective resolution l/Mf, and building 
it into a redefinition of the parton distributions, this is no longer a valid assumption, though we 
transfer the terminology. 
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to single-scale observations such as the invariant mass Q of muon pairs [27, 31- 
34], physicists looked at correlations between the muons, such as the transverse 
momentum imbalance, that is, the total Qt of the pair (here and throughout the 
text, we'll use "Qr" to denote the transverse momentum of a system of two or more 
particles, and "/c^" oi ^^px" for single particles). To a good first approximation, at 
relativistic speeds along a single beam line, the constituents of hadronic interactions 
present no transverse momentum relative to each other or the beam, and thus at 
Born level (just the Drell-Yan subprocess), there should be no imbalance, and we 
get a delta function at total Qt = 0. ^ The available data [37] indeed showed the 
cross section to be largest at Qt = 0, falling off gradually (although still steeply) 
with increasing Qt- This data being inclusive, the contribution at Qt 7^ could 
be due to "recoil" of the pair against bremsstrahlung from the initial state partons, 
recoil against a final-state quark in the quark-gluon "Compton" subprocess described 
above, or some nonperturbative, "intrinsic" Pt of the incoming partons (see figure 
1.5) . The former having not yet been calculated, it was unclear how much of a 
contribution was to be ascribed to each of these effects, and the approach varied 
among researchers. 




Figure 1.5. Intrinsic transverse momentum. 

"^Indeed, the "naive" parton model takes this as a starting point; the current "effective 
constituents" interpretation recognizes that there are no natural cutoffs that would limit the 
transverse momenta of partons relative to their parent hadrons, and hence we're led to the idea of 
"intrinsic" parton pr [35-36] . 
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A leading log calculation with simple Drell-Yan as the Born term admitted only 
intrinsic pt, and thus some attempted a description solely through pj'-scaling of 
the parton distribution functions [38-39]. This would "smear" the zero-Qr delta 
function into a more physical, Gaussian shape, and it was hoped that the amount of 
intrinsic px needed would be small. This turned out not to be the case [40]. When 
in fact the next-to-leading order bremsstrahlung and Compton subprocesses were 
calculated, although the large-Q^ region seemed to be sufficiently accounted for, an 
unacceptable 1ow-Qt prediction remained. Some authors simply left their results 
as valid only in the high-Qr region [41-42]. Others, notably Altarelli, Parisi, and 
Petronzio, were able to devise a model of intrinsic pt which also brought the Iow-Qt 
distribution under control [29, 43]. The average intrinsic px required, however, 
was still a matter of phenomenology, not derivation. In the same way that the 
factorization and renormalization scale dependences of a leading-log calculation are 
reduced when next-to-leading (and higher) contributions are added, higher-order 
calculations were needed which focused on transverse momentum effects. 

It was recognized that the incoming quarks could radiate any number of unde- 
tectable gluons, if the "jets" of hadrons thus produced had energies less than the 
trigger threshold, and that this additional, soft radiation could collectively influence 
the transverse recoil of the observed flnal state. An initial attempt to attack these 
contributions was undertaken in 1978 by Dokshitzer, D'Yakanov and Troyan (DDT) 
[44], who showed that in the soft limit, a pattern emerged as one went to higher 
orders. Careful inclusion of n gluons into the flnal state shows that even after the 
cancellation of the infrared divergence at Qt = 0, there remains a series of logarithms 
at order n in the strong coupling of the form ^ 

^Tlie "+" designation can here be thought of as a simple reminder that the pole at Qt = is 
removed. Later, in Section 2.4, we'll introduce "plus-distributions" more formally. Neither these 
logs nor the leading order delta function at Qt = are sufficiently physical representations of an 
actual Qr-distribution; their combination into such a representation is the goal of resummation. 
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1 2n-l 

E [in(QVQT^)]: • (1-1) 

At first glance, more gluons seems to make the problem worse. At high Qt (close 
to the mass Q), these contributions are small. At low Qt-, however, we have a two 
scale problem, the contributions of which may easily overwhelm the smallness of a^". 
The convergence of a perturbative expansion is thus thrown into doubt. Furthermore, 
as was shown by Altarelli, Ellis, and Martinelli (1979) [33], this is the case no matter 
how you define the parton distribution functions. 

That a pattern existed, however, with at least the largest (m = 2n — 1) logs 
sharing a common coefficient across orders, allowed DDT to resum this tier of 
logarithms over all orders, yielding a finite result at Qt = 0. Their method, 
flawed by a lack of regard for the conservation of momentum among the emitted 
gluons, was improved upon by Parisi and Petronzio in 1979 [45-46]. The new 
procedure used the momentum-conserving delta function to Fourier transform the 
cross section into impact parameter space, perform the sum there, and transform 
back to momentum space. The amount of intrinsic pt needed to reproduce the data 
was drastically reduced, though the extremely low-Q^ region still required some such 
nonperturbative help. 

A later analysis (1981) by Collins and Soper [47] for the reverse Drell-Yan 
process uncovered a systematic way to include subleading logs, thus improving the 
perturbatively-tractable region, and provided a model for non-perturbative effects, 
allowing the prediction to extend to even smaller transverse momenta. In 1984, 
Altarelli, Ellis, Greco, and Martinelli [48], along with Collins, Soper and Sterman 
[49], reapplied this formalism to the general case of intermediate vector boson 
production in hadronic collisions, and it became evident that reduction of the 
required intrinsic px was linked in large part to systematic addition of further tiers 
of subleading logs. This reflected the fact that ultimately both intrinsic and recoil 
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contributions derive from parton interactions, whether internal or external to the 
unobservable "dissociation" of the hadron. 

1.5 Goals and Outline of this Work 

Although the bulk of resummation work has been devoted to Drell-Yan and other 
intermediate vector boson processes, the methods have recently been shown beneficial 
in the prediction of other processes with steeply-falling transverse momentum spectra 
[50]. These include double jet (dijet) production [51] and direct double photon 
production [52-53] . 

The study of hadronic photon production is important for a number of reasons. 
Single photon production, due to its leading-order dependence on the gluon content 
of the incoming hadrons, serves as a good test of our knowledge of these densities. By 
comparing the results of single and double photon production, the strong coupling 
may be evaluated. Finally, in the ongoing search for the Higgs, theory predicts a 
two photon decay, making precise calculation of the direct reactions necessary for 
elimination of background events. 

However, the focus of most of the literature has been on observed systems of 
two or more particles; less well-known is the effect of resummed soft radiation on 
single-particle px spectra. This is not surprising, as the resummed logarithms arise 
naturally only when asking questions about the properties of multi-particle systems. 
Any transverse "kick" from unobserved radiation, though, should distribute itself 
among the observed particles, and it is only natural to ask about this distribution as 
well. In this dissertation we apply the current resummation techniques intermediately 
to a direct photon plus jet final state, and show that the kinematic effect on the 
photon pt is significant even after integration over the jet. 

This requires that we keep track of more variables, and, unlike many previous 
calculations, refrain from integrating over them analytically. This is beneficial from 
an experimental standpoint as well. In order to extract useful signals in practice, 
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it is sometimes necessary to make kinematic cuts on the detected particles. These 
cuts may not be expressible analytically at the parton level, where the integrations 
need to be done, as there is usually a complicated mapping among the associated 
variables. A multitude of convoluted holes may appear in the partonic phase 
space, thus making of little versatility formulae in which the relevant variables have 
already been integrated over a continuous, predetermined range. In addition, jet 
definitions, detector acceptances, and kinematic cuts are often changed, necessitating 
full recalculation of analytic predictions. 

We will handle experimental cuts by performing the required integrals via a Monte 
Carlo approach. By taking advantage of the method's ability to calculate Jacobians 
automatically, changes in applied cuts are easily accomodated, and added benefits 
include the ability to calculate any number of observables simultaneously [54-55]. 
Precedence has been set for the application of Monte Carlo to photon physics, as 
Bailey, Owens, and Ohnemus [50, 56] performed such studies in the early 1990's, 
incorporating a synthesis of NLO analytic and Monte Carlo methods to predict the 
photon pt distribution. Inclusion of soft gluon effects is the focus of the current 
work, with the expectation that the resultant method will have applicability to other 
single-particle predictions. 

In Chapter 2, the fundamentals of QCD perturbation theory are presented, start- 
ing with the development of the Lagrangian and the Feynman rules for calculating 
subprocess matrix elements. We discuss the singularities that arise in different types 
of diagrams, and the manner in which they disappear in the full theory via the process 
of renormalization. Finally, we describe the QCD hard scattering formalism, which 
combines the perturbatively calculable matrix elements with the experimentally 
determined parton distributions to arrive at predictions for observable quantities. 

Chapters 3, 4, and 5 together present the evolution of current techniques for the 
resummation of those large residual logarithms which arise naturally in perturbation 
theory, but are problematic for its convergence. Chapter 3 motivates a statement 
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of the problem and offers the first attempts at a solution. Chapter 4 describes the 
efforts of Collins, Soper, and Sterman to derive a more comprehensive picture, and 
solution via a transform to impact-parameter space. Chapter 5 reviews both the 
manner in which this formalism may be analytically simplified and its connection 
with the idea of intrinsic parton transverse momentum. 

Using the diagrams and corresponding matrix elements of Appendix A, in Chapter 
6 we begin the process of calculating the pieces of a NLO direct photon (plus jet) 
cross section. Contributions from three-body final states are approximated in the 
singular regions of phase space, and these approximations are both subtracted from 
the original (to provide one finite piece) and added to the virtual two-body terms 
(which are singular but of the opposite sign) to provide another finite piece. Along 
the way, we perform a number of analytical integrations, collected for reference in 
Appendix B. Casting these finite pieces in a form congruent to that of the NLO 
expansion of our resummation scheme, we then simply read off the coefficients of the 
original resummed form. Detailed calculational results are presented in Appendix C. 
Finally, we present a method for matching the resummed and perturbative results 
at the boundary between their regions of applicability. 

Crucial to the versatility of this calculation is the use of Monte Carlo techniques 
for performing remaining integrals, and a short treatise on these techniques is 
included in Appendix D. With the resulting FORTRAN program (available upon 
request), we produce physical predictions for the single direct photon px spectrum 
at a variety of energies; these results, and general conclusions, are presented in 
Chapter 7. Observations on the form of the nonperturbative parametrization are 
also discussed there. 
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CHAPTER 2 



QCD PERTURBATION THEORY 



2.1 The QCD Lagrangian 

We begin by constructing the QCD Lagrangian, from which the equations of 
motion and scattering matrix follow. The Einstein summation convention will be 
used throughout, in which repeated indices are to be summed over. 

The quarks will be represented by Fermi-Dirac fields ipf^i{x) of flavor /, color i, 
and mass rrif. Such fields obey the Dirac equation, and so for a system of free quarks 
we would have but 

C, = ^Pf,Srd,-mj]^pf'\ (2.1) 

Here the flavor index / runs over the number of flavors Nf, while the color 
index i (and later, j) runs over the number of colors N^. Once we impose local 
SU{Nc) gauge invariance, the existence of gauge boson fields (the gluons) A^(x) 
becomes necessary, as we must replace the gradient with the covariant derivative 

= 8^6"^^ + z(7(T(j)*"'y4^. There are — 1 such fields, corresponding to the number 
of elements in the adjoint representation of the symmetry group, and we'll use the 
indices {a, b, c} to distinguish them. The generators Ta form a basis for the Lie 
algebra of the SU (Nc) symmetry group, and obey the commutation relations 



i fabcTc 
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(2.2) 



{Ta,Tt} ='-f + dabcTc. (2.3) 

The fabc are totally antisymmetric, and are called the structure constants of the 
group. The dabc are totally symmetric. 

Now that we have gluons, we need to include a kinetic term for these fields which 
preserves local SU (Nc) invariance. With field strength tensors given by 

F^u = d,Al - d^Al + gfabcA^Al , (2.4) 
a suitable construction is simply 

= -^Fr^;, . (2.5) 

In contrast to the photons of QED, an Abelian theory, the gluons of QCD couple to 
each other, as indicated by the third term in equation 2.4. As in QED, however, the 
boson fields as currently written have extra, unphysical degrees of freedom. We can 
get rid of these by imposing a gauge-fixing condition, and there are many options, 
not all of which preserve relativistic covariance. The Lorentz condition, d^A^^ = 0, 
does however, and leads to the Lagrangian term 

Cf.. = -^{d,A^)\ (2.6) 

Use of a covariant gauge is not, unfortunately, a sufficient constraint. The 
longitudinal polarization of the gluon fields remains, and if not dealt with can destroy 
the unitarity of scattering amplitudes. Fortunately, as was shown by Faddeev in 1967, 
one can couple a set of scalar ghost fields (pa to the gluons in such a way as to cancel 
the unwanted degrees of freedom. As there are — 1 such degrees of freedom (one 
per gluon field), there must be the same number of ghost fields, although the coupling 
is not purely color-matched, as the required addition to the Lagrangian shows: 
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■ghost 



(2.7) 



The final QCD Lagrangian in a Lorentz gauge is thus given by the sum 




(2.8) 



2.2 Feynman Rules for QCD 



By quantizing the fields in the above Lagrangian, we can derive quantum- 
mechanical amplitudes for particular field interactions and propagation from one 
point to another in spacetime. From these, physical predictions such as cross sections 
and decay rates can be found. Quantization is not a unique procedure, but all 
methods lead to the same physical predictions. In the traditional canonical operator 
formalism, the fields are identified with quantum operators which obey certain 
commutation relations. The Green functions from which a particular amplitude are 
built are calculated as vacuum expectation values of products of these operators 
at particular spacetime points. In the functional-integral approach, the Green 
functions are obtained by integrating the product of the fields over all their possible 
intermediate forms, with a Lagrangian-dependent weight. A further method is that 
of stochastic averaging, in which the fields are identified as stochastic variables and 
the Green functions averages of field products in equilibrium. 

In momentum space, the Feynman diagrams for QCD and their corresponding 
Green functions are shown in figure fig:feynrule. Indices i and j run over the Nc 
colors, indices a through d (as well as all Oj) will denote the — 1 members of the 
adjoint representation of color SU{Nc), Greek letters fi and u (as well as all /Xj) are 
4- vector indices, and momenta will be labeled by pi and fcj. Spin and polarization 
indices are given as A. In the denominator of the propagator factors, it is customary 
to add a small imaginary amount ie in order to satisfy causal boundary conditions, 
but this will be ignored for clarity. There is some freedom in assigning factors of i 
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and minus signs with these diagrams; the particular convention used here is that of 
Muta [57]. 



"a(p) 



■"a(p) 



fc) 



gluon 
propagato 



!.d„„(fc) 



propagato 



quark 

propagato 



• •« • 




" "f<lf<2f3A»4 vortex 



gluon — quark 



Figure 2.1. Feynman Rules for QCD. 



The functions d^^, V^^f,^^^, and W^^l^l^f^^ are given below: 
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y^a^a^a^a^ ^ / .13,24 _ H4,32n i f f 12,34 _ .14,23n 

+ - f''''')9,.,49,3,2 , (2.9) 

in which denotes the combination 

and the Minkowski metric (in D dimensions) is taken to be 

D-l 

1=1 

There are additional requirements when loops arise in diagrams. First of all, a ghost 
loop must be included whenever a closed fermion loop appears. Secondly, for all 
loops, the loop momentum k must be integrated over: 

d^k 



i2.f ' 

with fermion and ghost loops carrying an extra minus sign. Lastly, additional 
symmetry factors must be included for the types of diagrams shown in figure 2.2. 

2.3 Renormalization 

2.3.1 Motivation 

In the calculation of the amplitude corresponding to a particular Feynman 
diagram containing internal (virtual) particles, there will necessarily be unobserved 
degrees of freedom which must be integrated over. For example, a diagram with 
two external legs and one loop will have two characteristic momenta, one (p^) for 
the external legs and one [k'^) for an internal leg (or propagator) of the loop, the 
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Figure 2.2. Symmetry factors. 

other loop propagator momentum being fixed by momentum conservation. Tlie loop 
momentum k'^, being unobserved, must be integrated over, but by the uncertainty 
principle, it can be arbitrarily large as long as the "size" of the loop is correspondingly 
small. Thus the potential arises for the integral to diverge. Each propagator will 
contribute a power of 1 or 2 to the denominator of the integrand, thus lessening the 
degree of divergence, but if the dimension of the numerator is greater than that of 
the denominator, we will still have a problem. By the same token, if the rest mass 
of the loop particle is zero, there will be an infrared divergence at the low end of the 
integral, k'^k^ = 0. 

As noted in Section 1.3, there are divergences even in purely tree-level graphs, 
and regularization is the general term used to describe the process by which 
divergences are parametrized and the amplitudes written as explicit functions of 
these parameters. Being arbitrary, such parameters must ultimately disappear in the 
prediction of any physical observable, but in the mathematical calculation thereof, we 
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can temporarily excuse their presence. One main class of regularization scheme is the 
cutoj^ method, in which ultraviolet integrals are stopped at some high scale or infrared 
divergences are avoided by introducing small masses. However, such techniques tend 
to destroy either the Lorentz invariance or the internal symmetries of the Lagrangian, 
and so are rarely used. More common is the method of dimensional regularization, 
in which integration over what would be a 4-dimensional momentum variable is 
analytically continued to a space of arbitrary dimension D, usually parametrized as 
D = 4 — 2e. The integral is now formally finite for e > 0, and the amplitude can 
be expressed as a function of e. In the end, ultraviolet poles in the regularization 
parameter e must either cancel or be absorbed into redefinitions of the observable 
masses and charges. This latter procedure, regardless of the means of regularization, 
is given the name renormalization. 

The Lagrangian of a theory will be a function of the fields and "bare" parameters, 
such as masses niQ and couplings cq- Derivation of the Feynman rules from such a 
Lagrangian will result in a set of "bare" Green functions, the simple propagators 
and vertices of the theory. Usually, the Lagrangian is written as a sum of terms, 
each of which corresponds to one of these functions. More complicated diagrams 
are built up from these, and can be grouped into two classes - those which have 
the same external legs as the bare diagrams and those that do not. Of the ones 
that do, some are single-particle-irreducible (simply-connected diagrams which have 
external propagators removed and which cannot be split into two by cutting one 
internal propagator). We denote these by the term "SPI", and for vertices, this is all 
we have in this sub-class. For propagators there will be other diagrams, but these 
can be shown to be just combinations of the SPI diagrams. We'll now define "full" 
propagators and vertices as sums of all diagrams with the same external legs as their 
corresponding bare diagrams. For vertices, this will then be the sum of all SPI vertex 
diagrams, while for propagators, we have to first sum the SPI propagators, and then 
link these sums in a geometric chain (see figure 2.3) [58]. We now have a set of full 
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Green functions, each in one-to-one correspondence with its bare counterpart. The 
remaining class of diagrams, those which don't have the same external legs as the 
bare versions, were still constructed from them, and so can be calculated now via the 
full Green functions, expanded to the appropriate order of the coupling. 



-iT,{p, mo) 




O. CO 



i 



^ - mo - ) 

Figure 2.3. Building the full propagator. 

In general, these Green functions will be UV-divergent, as they contain "clouds" 
of virtual activity. However, in a renormalizable theory, to each simple propagator or 
vertex there corresponds a characteristic quantity, potentially observable in the "full" 
case, unobservable in the "bare". For propagators this is the mass; for vertices it is 
the coupling. The correspondence between the full and bare Green functions gives us 
a means of dealing with these divergences, as well as a physical picture of how higher 
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orders of virtual activity "shift" the masses and couphngs from their bare values. 
The full Green functions emerge as divergent functions of the bare parameters, but 
if we define observable quantities such that the bare parameters are themselves 
divergent functions of the observables, then the Green functions can be rendered as 
finite functions of the observables. 

2.3.2 Canonical Procedure 

The exact process of renormalization depends on the theory under consideration, 
but (without getting into these details) one can describe, in a general way, the steps 
involved. We use one mass, coupling, and regulator for simplicity. 

1. Compute the full functions Fj to a given order of the bare coupling cq, using a 
regulator A when necessary. One then has the functions rj(eo, ?Tio, A, {pj}) to order 
A^, and these are divergent as the regulator approaches some limit A ^Hm- Here 
{pj} are the momenta of external legs. 

2. Impose renormalization conditions to define the form of the renormalized 
Green functions at a particular scale /i = f{{pj})- These forms should contain a 
dependence on the renormalized couplings and masses {e, m} and should be finite 
as A — * Kiim- Suitable forms are usually suggested by the bare Green functions 
themselves, e.g. the full propagator should be able to be brought into the same form 
as the bare propagator, the new, observable mass (now a complicated function of 
the bare parameters) taking the place of the bare mass (see figure 2.3). Imposition 
of such constraints on the Green functions Fj results in relations e(eo, jtiq. A) and 
m(eo,'mo. A) between the renormalized and bare quantities. 

3. Solve for the bare parameters {eo,'mo} and substitute into the functions Fj. 
The result is Fj(e,m, A), which will be finite as A ^ ^um- 

We'll call this the canonical renormalization procedure, and it will always work 
in a renormalizable theory, yet is cumbersome at high orders A^. An equivalent 
but better organized approach to the bookkeeping is provided by the counterterm 
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method, which anticipates the shift in parameter space by building it into the 
Lagrangian at the start. 

2.3.3 Counterterm Method 

1. Rescale the bare parameters (including the field strengths ipo) by introducing 
new, temporary unknowns Zi, which relate them to their observable counterparts; 
i.e. ijj = Z^ipo, e = Z^eo, etc. 

2. Trade in these multiplicative unknowns Zi for additive ones 6i in such a 
way that, upon substitution of these expressions for the bare parameters into the 
Lagrangian, one reproduces the original Lagrangian (now in terms of the renormalized 
observables) , plus additional counterterms. We'll refer to these pieces as Cqg and 
CcTi respectively. The counterterms will have corresponding Feynman rules of their 
own. 

3. Specify the renormalization conditions, again as desired relations of form 
between the unrenormalized and renormalized Green functions at scale /i. 

4. Calculate the full functions using the rearranged Lagrangian. To a given 
order A^, one will find the usual divergences from Ccgi plus new divergences 
from the counterterms. However, we can now get these to cancel, as the set of 
unobservable shifts 5i are free parameters. By adjusting these to maintain the 
renormalization conditions (which by construction define nonsingular amplitudes), 
we can consistently absorb the infinities of the theory into the bare parameters, 
leaving Green functions that are finite and dependent upon only the renormalized 
charges and masses. 

The counterterms in Cqt-, with suitable renormalization conditions, are commonly 
said to "subtract off" the divergences of the "original" Lagrangian Ccgi but there 
remains an ambiguity in this process, as one can always subtract an extra, arbitrary, 
finite amount as well. The size of this piece is, in general, related to the particular 
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renormalization conditions, and so together these two choices define a renormalization 
scheme. We will discuss the most commonly used schemes in a later section. 

As an example, we take the renormalization of QED. The bare Lagrangian is as 
follows (we denote bare quantities by the subscript "0"): 

C = -^{F^' + 4(^^ - mo)tP - eo^TM^o^^ • (2.12) 
We now impose the following sets of rescalings and variable changes: 





ryl/2 , 






= zl'^A^" 




ry ryl/2 


= eZi 


(2.13) 


8m. = 


Z2mo — m 




5l = 


Z,-l 




52 = 


Z2-I 




h = 


Z3-I, 


(2.14) 



and thus transform the Lagrangian into a sum of two pieces, one of which resembles 
the original Lagrangian, the other provides the counterterms. The bare quantities 
and intermediate scaling factors disappear: 

^ = ^CG + ^CT 

CcT = -^53(i^n' + ^(^52^ -U^-e5i^7/.^A^- (2.15) 

For renormalization conditions, we use the following: 
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Hij) = m) 











d'^ p=m 

n(g2 = 0) 







r^(p = p) 



(2.16) 



where the amphtudes F^, S, and 11^^ = {g^'^q^ — q'^q'^)Il denote the amputated full 
vertex, single-particle-irreducible (SPI) electron propagator element, and SPI photon 
propagator element, respectively. When calculated using the rearranged Lagrangian, 
which includes counterterms, each will, in general, be a function of the quantities 
{e,m,Si,S2,S3,6m,^ ,q}, and the regulator. Enforcing the above conditions then 
determines {6i, 62,63, 6m}- 

As we have said, the precise conditions used are a matter of convention; the crucial 
point is that we define the full Green functions (and thus the renormalized masses 
and couplings) to (1) be finite, and (2) to take a certain form at some renormalization 
scale fi or set of scales. In the QED example above, we used a momentum subtraction 
scheme, in which the scale used depends on the Green function being renormalized. 
We could, however, have chosen a single scale, common to all amplitudes, i.e. 

= — and q^ = —fi^. The counterterms then subtract only the divergent poles 
[Minimal Subtraction, or ^^MS"), or these plus standard constants Modified Minimal 
Subtraction, ^^MS"). In this dissertation we will use the MS convention. 

Either way, we are left with renormalized Green's functions G and couplings A 
which, though finite, seem to depend upon a new parameter (fi) not present in the 
original, bare versions [Gq, Aq). We know that these must be related by the (now 
calculable) rescalings Zf, for example, there must be a calculable function Zq of the 
scaling factors Zj such that 



G(/i,A) = ZG(/i,A)Go(A) 
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(2.17) 



However, we have calculated this relation to a fixed order only. How can we 
generalize this to account for not only the poles but the associated scale dependent 
pieces In^/i/p found at even higher orders? Specifically, if /i is taken large enough, 
how do we justify a perturbative expansion, unless the coupling ~ A^^ becomes 
correspondingly small? There must be a Renormalization Group Equation (RGE) 
which expresses the fact that a change in scale /i needs to be compensated in order 
to keep the bare Green's function (and the bare coupling) invariant. From equation 
2.17: 



or 



)i— + (3—+^]G{^i,\) = Q, (2.19) 



with 



I3(\) s f,^ (2^20) 

-'(^) 4f ■ <^-^^* 

Equation 2.19 is known as the Callan-Symanzik equation [59-60]. The functions 
(3 and 7 are unitless and must depend on the coupling only. Yet they describe the 
compensating shifts in coupling and field strengths (respectively) for a change in 
scale jji. So these strengths must themselves be functions of /i, and equations 2.20, 
2.21 lead us to the perturbative form of this dependence. Here we simplified to a 
theory involving one field and one coupling; in general, there will be one 7 for each 
field and one (3 for each coupling. 
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2.3.4 Application to QCD 

Now that we've described the basics of renormalization, let's list the relations 
appropriate to QCD. As expected, they are somewhat more complicated. We'll use 
a superscript "0" to denote the bare quantities. 



rrii 



2F 
[1% 



}^{g^Ai,)\x) 



(2.22) 



The running coupling = g'^/4:TT is a solution of the Callan-Symanzik equation 2.20, 
and comes out (to two loops) as 



a^(/i) 



27r 



Mn(/i7A2) 



2c ln[ln(^VA^)] 

b ln(/iVA2) ^^/''^ 



(2.23) 



in which 



llNr 



2Nf 



1 

86 



34, 



13, 



-NnN 



Nf 



As expected, this exhibits the property of asymptotic freedom 
coupling as decreases. 



(2.24) 

at large /i, the 
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2.4 Factorization 

With a renormalized theory, one can begin to calculate physical quantities and 
compare these predictions with experimental measurements. For a particular set 
of initial-state (incoming) reactants {a,b} and final-state products {c,d} (e.g.), one 
uses the Feynman rules to calculate an amplitude Aii for each topologically unique 
way this reaction can occur. Quantum mechanics then dictates that the probability 
for the reaction is proportional to the squared sum of these amplitudes: 

E|M|' = E-^*-^- (2-25) 

i i 

In addition, the squared amplitude refers to (and is a function of) particular 
values of the particle momenta {PayPb^Pc^Pd}- That is, the reaction rate per unit 
time into volume element d^pcd^pd is given by 

E I-^^I'7^%1^7^%1J(2-) Vb:^ +Pt-P'.- P'd) 
i (27r) 2p0 (27r) 

and one must integrate over this phase space in order to obtain the full reaction rate. 
Note that a delta function is included to conserve momentum. 

In the case of hadronic interactions, the initial state particles will be partons, one 
from each hadron. By the collinear geometry of the interaction, the parton momenta 
are known up to fractions Xa and Xi, of the incoming hadron momenta, and the cross 
section (an experimentally measurable quantity) for two partons, of fiavors a and 6, 
with these momenta, to produce the given final state, is found by dividing the above 
reaction rate by a fiux factor 2w(s,ma,mb), in which s = [p^ + PhY 

w{x, y, z) = \J x'^ + y"^ + z"^ — 2xy — 2xz — 2yz . (2.26) 

For all but the heaviest quarks, we can take the parton masses to be zero (relative 
to the hadronic center of mass energy S), and w{s,ma,mb) —>■ s = XaXbS. The 
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partonic cross section for a particular 2 — > 2 particle subprocess is thus given by the 
product 

da = — d T2 , 

2s 

in which the 2-body phase space factor is 

d'r^ = ftof2o (^-)'^'(P'^ +Pt-P'c- P'.) . (2.28) 
(2vr) 2p0 (27r) 2p^ 

The partonic cross section is perturbatively calculable to any order if one is 
willing to calculate the required amplitudes. However, it is not yet useful in practice. 
First, the above cross section pertains to partons of definite flavor, color, spin, 
and momentum fraction. In any real hadronic interaction, a multitude of partonic 
flavors will contribute, at all momenta and with all possible color and spin values. 
Spin and color averages are abbreviated by a line over the squared amplitude sum, 
e.g. However, we do not apriori know the momentum or flavor distribution of 
partons inside hadrons, nor do we have a deterministic model for the hadronization 
of final state partons. We are forced to find a way of factorizing the reaction into 
component steps, some of which are dealt with perturbatively, the rest requiring 
effective parametrizations and experimental input (see figure 2.4). 

In the parton model of high-energy hadronic collisions, the initial reaction takes 
place over such a small time scale that it becomes plausible to treat each hadron 
as a collection of free partons, each travelling in the same direction as the parent 
hadron, and only one of which takes place in any given "hard scattering" . That is, 
one may calculate basic subprocesses involving one parton from each hadron, then 
sum incoherently over the contributions of all partons. Similarly, the process of 
hadronization is assumed to take place over a time scale much longer than the hard 
scattering, and can therefore also be treated independently. 
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(2.27) 



Figure 2.4. Factorization of the hadronic interaction. 



Thus, we can assign to each initial state liadron a set of parton distribution 
functions qiji ^, and, if we are looking at particular hadrons in the final state, to 
each final-state parton a given probability Dujh to fragment into a state containing 
the desired hadron. The physically realizable cross section for the reaction is then 
obtained by weighting the subprocess cross section da with these distribution and 
fragmentation functions, evaluated for particular flavors and at particular fractions 
of the parent (or daughter) hadrons' momenta, then summing over all flavors and 
momentum fractions: 



^Parton distributions are variously denoted in the literature by qf or G/ (where / is the flavor), 
or fq (where q is the flavor). 




(2.29) 
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If suitable parametrizations of the distribution and fragmentation functions are made 
at the outset, these quantities may be fit to known cross section data and used in 
subsequent calculations. 

In the process of calculating da to higher orders, two main classes of diagrams 
arise. First there are vertex and propagator modifications, in which extra connections 
are made between the legs of a diagram, without changing the number of external 
legs. These produce ultraviolet (UV) and infrared (IR) singularities (or "poles") at 
the upper and lower limits of loop momentum, and in the previous section we have 
seen that the UV poles may be regulated and absorbed into the masses and couplings 
of the theory. The second class of diagrams involve the emission of an extra particle 
into the final state, and also give rise to IR singularities. Some of these arise in 
gauge boson emission at zero boson momentum, and cancel the IR poles from virtual 
contributions [24] . As well, there are collinear singularities at zero angle between 
the daughter and parent particle. 

The collinear singularities, arising as they do from the inability to distinguish a 
bare parton from one accompanied by a number of collinear partons, get moved into 
a redefinition of the distribution and fragmentation functions, usually accompanied 
by logarithms of a momentum scale characteristic of the subprocess. The evolution 
of these quantities is now dependent upon the change in momentum transfer at which 
they are sampled, a result quantified for distribution functions by Altarelli and Parisi 
[28], and for the fragmentation functions by Owens and Uematsu [61]. 

At order relative to the leading order, there can be as many as A^ extra 
radiated partons, any or all of which may be collinear. We thus find a power 
series of logarithmic divergences at the NLO and higher orders. In the leading 
log approximation, only the logarithm of highest power at each order is retained. 
Iteration of the evolution equations effectively sums the collinear logarithms to solve 
for the distribution and fragmentation functions, with the original momentum of 
parametrization as reference point. The leading log approximation dictates that we 
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can now simply use the modified coupling, distribution, and fragmentation functions 
with the Born matrix element. A number of calculations have been carried out in 
this way [63-65]. 

To show how singularity cancellation and factorization works, we'll look at 
a simple example, that of double photon production. We won't be looking for 
particular hadrons in the final state, so no fragmentation functions will enter into 
our calculation. 

The "Born' (leading order) subprocess cross section is calculated using the 
diagrams shown in figure 2.5. All contributions have already been regulated using 
dimensional regulation, with e = (4 — D)/2. To order a^, we have 

rf^o = ^ — I' '^-^H{x.,x,)d^T,nM\l , (2.30) 

2b Jo Xa Jo Xb 

where the luminosity function H (which contains the parton distributions gj), the 
squared, summed and averaged matrix element S|M|q, and the 2-body phase space 
factor d^r2 are given as follows: 



H{Xa,Xb) 



<lf/A{Xa)qf/B{Xb) + if ^ f) 

1 Mtt 
8^VT 



8iVc(47r)V^aVf(l-e) 
v{l - v)]' 



rfi 



-dv 



(2.31) 
(2.32) 

(2.33) 



Note here that the phase space factor has been simplified by use of the momentum- 
conserving delta function, and becomes dependent only upon a single variable 
V = I ( 1 + cos 6'3 ) , where 9^ is the angle makes with pi in the partonic center-of-mass 
frame. The Mande/siam variables {s, t, u} are then given by s = XaX^S, t = —{l—v)s, 
u = —vs. 

In a leading log calculation, we would stop here, take e ^ (as there are no 
divergences), and evaluate the distribution functions at a factorization scale Mf 
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Figure 2.5. Leading order diagrams for qq ^ 

characteristic of the subprocess (perhaps the mass Q = of the photon pair). If 
this leading-order term had contained a dependence on the running strong coupling 

we would evaluate this, too, at a renormalization scale /i, usually taken to be 
equal to the factorization scale. In fact, we will be adding higher-order terms, so the 
strong coupling will show up in the final result. 

At next-to- leading order, two things happen. First, we will have contributions 
from tree diagrams in which a gluon is radiated from one of the quark legs ^ 
Second, we will have virtual diagrams in which a gluon is radiated and reabsorbed, 
so that there are still only two particles in the final state. Of course, emission 
and reabsorption amounts to two QCD vertices, whereas the tree diagrams only 
contribute one, so if we're working to NLO only, we keep only the interference terms 
between the Born and virtual amplitudes. 

The tree diagrams for q{j>i)q{p2) l{p^)'l{pi)9{'P^) ^"^^ shown in figure 2.6, and 
their contribution is as follows: 

da, = ^ /' — — x,)d^T,nM\\ , (2.34) 
2d Jo xi Jo X2 

where here 

^We will also have contributions in which a gluon initiates the process, and a quark appears in 
the final state, but we'll ignore these for the sake of this argument. 
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H{xi,X2) = qf/A{xi)qf/B{x2) + if ^ f) (2.35) 




(2.36) 



(2.37) 



Figure 2.6. Diagrams for qq — > 77(7. 
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The function Mi{p'r^,p^, xi,X2, e) is a complicated one, and this makes integrating 
over the extra gluon (pg) difficuh. As we are here interested only in the singularity 
structure, we will make an approximation. In the limit that the gluon becomes 
collinear with an incoming quark (psHpl, e.g.), its momentum, as well as that of 
the daughter quark (p^, which enters the subprocess) are related to the parent in 
terms of a simple fraction z. That is, = (1 — z)pi and = zp^. The matrix 
element simplifies, the integrals become easy, and the remaining phase space is that 
of a two-body final state. 



In short, the matrix element factorizes into a Born term describing a 2 ^ 2 
subprocess, and a splitting function P^^ describing the probability for a quark to 
split into a collinear quark-gluon pair. This latter function is independent of the 
particular subprocess, as long as it contains an incoming quark. Of course, we get 
a pole from the collinear singularity, and an important multiplicative factor. The 
result, in D = 4 — 2e dimensions, is: 



< 



Figure 2.7. Pi ■ P5 = collinear kinematics. 



H{Xa/z,Xb) 



qf/A{.Xa/z)qj/B{,X2) + {f ^ f) 



(2.38) 



EIM 




(2.39) 
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d^'T^—i — 



l/47ryr(l-e) s dz{l-z) 



l-2e 



e V s y r(l - 2e) (47r)' 4 



y3-2e 



(2.40) 



which gives the cross section 



da, = %r^r^rf-r,siMi 

25 JQ Xa Jo Xb 



X 



^47r/i2N 


/r(l-e) /•! rf^ 


z 




/ r(l - 2e) z 


1 - z. 



2e 



HixJz,Xb)P'{z,e) . (2.41) 



Here Pgg{z, e) is the aforementioned sphtting function, pertinent to the process 
q — > qg, where the quark on the right-hand side is involved in the QCD subprocess. 
This function will always appear in such a coUinear limit, regardless of the subprocess, 
and has the form 



1 + z' 



eil-z) 



(2.42) 



Unfortunately, there is also a soft pole here at z = 1. We'd like to separate 
this so we can focus on the purely collinear contribution. We start by adding 
and subtracting, under the z-integral, the z ^ 1 limit of H{xa/ z,Xh)Pqg{z,e). 
Remembering that all of this multiplies the collinear — pole, we have 



-1 /•! dz 

e Jxa z \1 



H{xa/z,Xb)P'{z,e 



2CFH{Xa, Xb) 



qq\ 

^ dz / z 



l-z 

2CFH{Xa,Xb) 



(2.43) 



^ z \1 — z J 1 — z 
This allows us to express the first integrand (which is finite as 2; ^ 1) in terms of a 
plus- distribution, a handy notational device which is defined upon convolution with 
a smooth function. We will have plenty of use for these in subsequent chapters, as 



well. Here, the plus-distribution (1 



is defined over a function f{z) as 



dz- r 

[l-z) 



1 — z 



(2.44) 



The first term in equation 2.43 then works out to 
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/ —H{XaZ,Xb)PZ,{z)^ 

e Jxa z e 



21n- 



1- Xr, 



+ 0{1) 



where is the plus-distribution form of the e-independent qq sphtting function: 



9<J 



1 + z' 



[l-z) 



(2.45) 



the second term becomes 



2CFH{xa,Xb) I I , l-a^a , 
— - In + C e 

e L2e Xa 

This operation sphts dai into two contributions, the first of which we'll call dacoii 
and the latter dasoft- The same terms arise, of course, when we look at the limit in 
which the gluon is radiated from the other incoming leg; in that case we'll be looking 
at terms involving H{xa,Xb/z). When we put everything together, we'll add in these 
contributions. 

Meanwhile, the (renormalized) virtual contribution is: 



da. 



virt 



E/ dXa 



2S Jo Xa Jo Xb 



''^'^'H{xa,Xb)d''T,f:\M\l 



where 



(2.46) 



S|M|; = S|M|; 



1C p "iC p 



+ o{i) 



T(l-e] 



as /4:nfi 

2^\~Y J r(l - 2e) 



(2.47) 



e 

The virtual diagrams are shown in figure 2.8. 

In all, then, we have one LO and three NLO contributions, which all add to give: 



where 



= 2S 



l^dxa f'^d^r,E\M\lH{xa,Xb 

Jo Xa JO Xb 

i^coll 



rfi- 



1 + 

e27rV s J r(l-2e) 



^soft + n 



virtl 



(2.48) 
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n, 



coll 



1 dz H{Xa/z,Xb) +, - 

„ z H(xa,x,) ^^^^^ 



1 - Xal - Xb 



dz H{Xa, Xb/z) 
Z H{Xa,Xb) 



Pliz) 



qq 



+ 3Cf + 0{e) 



n 



soft 



2Cf 
2Cf 



- 3Cf + 0{e) . 



(2.49) 
(2.50) 
(2.51) 



One can see by inspection that, as predicted, all the poles cancel except those 
associated with convolutions over splitting functions. However, if we were to redefine 
our quark distributions such that 



,K..M|, . (2.52) 

then we could absorb the remaining divergence into the distribution functions. The 
remaining NLO contributions would come solely from the finite 0{e)/e = 0{1) 
remainders, we could take e — > everywhere, and the distribution functions (as well 
as the running coupling, to this order) would then be evaluated at a factorization 
scale Mf, characteristic to the observed system (in our case, s). 



dcr 



dXa dXb ,4 



-d-'To_^M\„H(x,,xi„M,) 



(2.53) 



25* JQ Xa Jo Xb 

This reflects the interpretation noted above, that a parton which participates in a 
subprocess is in fact a component of (a component of a component of...) the parent 
hadron, and at higher energies, we resolve more and more components. Of course, 
we have here shown only the first term in such an expansion; at higher orders (more 
splittings) we must absorb terms with more poles, higher orders of the coupling, and 
more complicated splitting factors. 
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Furthermore, we have shown only the q ^ qg sphtting function (boldface type 
referring to the parton participating in the subprocess). In general, there are splitting 
functions for q gq, g qq, and g ^ gg. 

The above redefinition is in fact a combination of two separate reorganizations; 
the first absorbs the poles (and an arbitrary constant), the second relates to energy 
dependence. Expanding the order factor above, we find 



Aqfi 



a. 



X] 



1 dz 



qf{x/z)P^ 



-1 M? 

- - In 47r + + In ^ - 0{1) 



(2.54) 



2tt Jx z 

where we have included the —f^ = C(l) terms from 2.48; that is, the entire remainder. 
•je is the Euler constant, ~ 0.577216. 

We can now absorb the pole and as much of the finite remainder as we like into a 
new distribution function qf{x). In the DIS scheme, everything is absorbed but the 
energy log: 



01 dz r — 1 
q^ix) = qf{x) + / -g^(x/z)P+(z) ln47r + 7^ - 0{l) 



In the MS scheme, only the factor ln47r + 7^; is subtracted: 



(2.55) 



2^ 



dz 



qf{x/z)P^_ 



-1 



In An + 7£; 



(2.56) 



In this dissertation, we will be using the MS scheme exclusively. 

Either way, the energy scale logarithm which remains is just as ubiquitous as the 
constants — ln47r + 7^; we have already absorbed, so we continue with the further 
redefinition 



qj{x, Mf) = qj + Aqj{x, Mj) , (2.57) 

where 
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Table 2.1. Splitting function pieces. 



qq 
(19 
9Q 
99 



Cf 

1 

2(l-e) 

Cf 
2Nc 

Q) 



(1-^)+ 

Z^ + {1 



+ z(l-z) 



rij 
§5(1 



U{z) 



12 



-1 



AgKx,M;) = ^(InMj - In/i^) /' -g/(x/.)P+( 



(2.58) 



In the limit of infinitesimal change of scale, and taking into account the other 
possible splitting functions, we are led to the aforementioned Altarelli-Parisi equa- 
tions: 



dqf{x,M^] 

dqg{x,M^] 
dlnM^ 



^ -[gK^MM^)P+(^) + g,(x/.,M2)P+(.)] (2.59) 
as dz 

2lT Jx z 



qf{x/z, M')P+{z) + q,{x/z, M')P+{z)] . (2.60) 



For reference. Table 2.1 lists the terms involved in first-order splitting functions; the 
conventions used in this dissertation are: 



1. The most general version is Pj^' = p^j 
daughter, j the parent. 



ptj + pfj + ^p]j 



where flavor i is the 



2. Non-plus-distribution versions are obtained from the above by simply 
removing the "+" -signs. 

3. "Primed" functions include the e-coefficient p^, non-primed functions do 
not. 
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4. "Hatted" functions do not include . 



5. Plj = PijPij includes only the coefficient of e. 



Examples (flavor labels suppressed): 



P = p'p 

P = P°(P + /) 

P+ = p'^p^ 

P+ = P°(p+ + /) 

P' = p\p + ep') 

P' = p%p + p' + ep') 



,+ 



P' = p\p+ + ep') (2.61) 

Thus we can derive the scale dependence of the parton distributions, and take 
into account the leading dependence at higher orders by iteratively solving the 
Altarelli-Parisi equations. Just as in the case of renormalization, a full solution 
requires that at some scale Mq we measure the distributions by comparison with 
experiment. Use of the new distributions in cross section calculations now requires 
only that we calculate, but then drop (subtract) collinear poles (and attendant finite 
pieces appropriate to the chosen factorization scheme). 
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CHAPTER 3 



EARLY RESUMMATION SCHEMES 



3.1 Origin of Terms in the Unresummed Distribution 

Here we will look at how the relevant logarithms arise for the QED process 
e~^e~ — ^ fi'^fi'^ . We will then describe the evolution of a viable resummation method, 
and apply it finally to the more complicated QCD process pp fJ'^fJ''^ + X. 

As discussed in Section 1.4, if we define = {p'^+ + p'^-)'^ as the mass (squared) 
of the observed muon pair, and Q"^ = {pth+ +'Ptii-)^ transverse momentum 
(squared) of the pair, calculation of the leading-order contribution to a process of 
this type results in a delta-function at Q\ = 0. That is, the Qr-distribution is of the 
form 



S{Qt') (3.1) 



ctq dQ'^dQ'^ 
where (Jq = da/dQ"^. 

In the relativistic limit, adding a single radiated boson to one of the incoming 
fermions modifies the matrix element thus: 

[p — ky 

Here Aiu{p) is the Born matrix element, p is the incoming electron momentum, and 
k is the momentum of the emitted massless boson, which we take to be a photon 
(see figure 3.1). For a high-mass muon pair, we may take the soft radiation limit [k 
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Figure 3.1. A single radiated photon in the initial state. 

small relative to p), and ignore terms involving k in the numerator of equation 3.2. 
This leaves us with the factorized form 

A^-^-»A<^"(.). (3.3) 

By momentum conservation, Q"^ = {p ■ k) , and we find a radiative correction 

of order a to the Born cross section: 



'^'^ S{Qt') + ^{HQVQI) + 0{1)). (3.4) 



aodQ^dQ^ ' 2txQI 

This cross section obviously diverges as Qt — ^ (see figure 3.3, dot-dashed line), 
a result of the soft and collinear divergences encountered when, respectively, the 
photon is either emitted with zero energy or is emitted in the same direction as the 
incident electron. In either case, p - k = 0. In a full calculation, these divergences (or 
poles), will be cancelled by the contribution of other diagrams, and the remainder is 
usually described by means of a "plus "-distribution; that is, a distribution the value 
of which at any part of the domain is meaningful only when convoluted with another 
smooth function (see Section 2.4 for one example). 

In our case we have a divergent function u^Qt) given by 
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^(Qt) = ^MQVQt) + 0(1)) . (3.5) 
We wish to separate off the singularities at = 0, which will be cancelled, and 
retain a function u_^.{Qx) which is finite and integrable over the entire range of Q'^, 
i.e. < < Q^. That is, we want 



^{Qt) ~^ ^+{Qt) + singularities, (3.6) 
with the plus-distribution z/+ defined as follows: 

/"^ dQlu4QT)f{QT) = r dQliy{QT)[f{QT) - /(O)] . (3.7) 
Jo Jo 

At = 0, fiQr) ^ f{^) =0 and the integral is defined, even if the shape of z/+((5t) 
is non-physical (like the Born term delta function). We will thus rewrite equation 
3.4 in the form 



(To dQ'^dQj^ 271 



i^+{Qt) = 


Qt 


+ 

+ 


rO(l)i 
■ ■ 


+ 


(3.8) 


If we now consider multiple radiated photons i 


see fi 


gure 3 


•2), 


we find that the ma- 



trix element — and thus the cross section — also factorizes in this soft approximation; 
that is, for n-photon emission we obtain a matrix element of the form 



and a cross section 



{a/ 271 



(Jo dQ'^dQj' 



where z/_|_(A;t^) is the first-order contribution for a single photon (here 



i=i i=i 



and the factor n\ arises from the indistinguishability of the emitted photons. 
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Figure 3.2. Multiple radiated photons in the initial state. 

3.2 Development of a Resummation Scheme 

We have seen how the series of large logarithms (equation 1.1) arises at each order 
in a. A first attempt at an all-orders result would naturally be to sum over just the 
largest log at each order. Our cross section contains a momentum-conserving delta 
function, but given the presence of n!, ignoring this function will allow us to directly 
form an exponentiation in momentum space. We now have completely independent 
photon emissions. If, as we have assumed so far, each has a relatively small transverse 
momentum, limited by that of the muon pair {Qt), when we sum over all orders: 



1 da _ d {a/27ry 
d ^ {a/2Txy 

2 



dQ 
d 

dQi 



n 

i=l 



v+{kTl)dkT. 



V+{kT^)dkT^ 



T n=l 



(3.11) 



we find, for u+^kx'^ 



lll(s/fcT^)+0(l) 



the result 



da 



(Jo dQ'^dQj^ 



dQr^ 
a 



-^ln2(QVQi,)-#ln(Q2/Q|) 



In ^ + 0(1) 



e 



^ln2(Q2/Q?,)-^ln(QVQ?^) 



(3.12) 



in which we have used the fact that the coefficient of our plus-distribution is constant: 
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/ ^ V+{kT^)dkT^ = [ MkT^)dkT^ - [ vikT^)dkT^ 

Jo Jo JqI 

= / v{kT^){l - l)dkT^ - v{kT^){l)dkT^ 

Jo Jq%. 
= O-iln^^-O(ln^). (3.13) 

This method was first attempted in 1978 by Dokshitzer, D'yakonov, and Troyan 
(DDT) [44]. ^ It gives, as we see in figure 3.3 (dashed fine), an exponential Sudakov 
suppression of the Qt distribution as 0. This effect is a consequence of the 

neglect of momentum-conservation among the emitted bosons, effectively imposing 
the condition that all of them have zero kx in order for a zero-Qr niuon pair to be 
produced. The exceedingly small phase space for n bosons to be produced, all with 
near-zero kx at the same time, leads to the suppression. We know, however, that 
this cannot be correct. It should be possible for two or more bosons of non-negligible 
kx to be emitted and still maintain a zero total Qt- 

A subsequent analysis by Parisi and Petronzio [45-46] kept the factorized form 
for d^a, but also kept the full momentum-conserving delta function and first-order 
contribution. This time, the exponentiation develops automatically in impact- 
parameter (6) space. The delta function can be rewritten as a Fourier transform 
in this space: 

S\Qt -Y^kTi) = 7:^2 f d%e-^'-^^^--T.yT.) . (3.14) 
i (^tt) J 

The cross section for n-boson emission then becomes a Fourier transform and 
inverse-transform, with the exponentiation sandwiched inbetween. To start. 



^Actually, they rcsummcd only the leading logs, m = 2n — 1. We include the 0(1) terms here 
as well. 
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(27r) 
1 



n\ 27r 

J2^ „~*Qr [^^+(^)] 



e 



(3.15) 



where z^+(&) is the transform of u^^kx) in impact parameter space. 

The next step is clear. Summing over the contributions from all numbers of 
photons, and including the 5-function at Qt = 0, Parisi and Petronzio arrived at 



da 



1 



(To dQ'^dQ\ 



{2ny 
1 

(2^ 



.n=l 



E 



a 



n\'2n 



(3.16) 



where the last exponential on the right is sometimes denoted by a. The prescription 
is, then, a Fourier transform of an exponentiated transform of the next-to-leading 
order (NLO) perturbative result. For our simple QED example, this latter transform 



is 



= r,^rM(e«'_i)rHO!ZM0±^l, ,3,17) 

JO JO 2 ^ ^ 

where the upper limit on kx has been relaxed relative to DDT and the plus- 
distribution definition 3.7 again used. In the high-6 limit, this goes as v{h) ~ 
-\n^{l+h'^Q'^), giving 



da 



(Jo dQ'^dQ 



Jo 



3 7r/2o 



constant. 



(3.18) 



The Sudakov suppression disappears, and one regains a finite, non-zero cross section 
at Qt = (see figure 3.3, solid line). 
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-0.05' ' ' ' ' 

QT/ Q 

Figure 3.3. Resummation brings tlie low-Qj^ region under control. 

Resummation thus allows one to include the effects of multiple soft boson 
radiation to all orders, while keeping the Iow-Qt region under control. The original 
perturbative prediction, necessarily unphysical as a result of the delta- functions and 
plus-distributions which dictated its bookkeeping, nevertheless was integrable: 
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Resummation maintains the normalization while providing a more physical model 
of the shape over the allowed Qt range. Effectively, the delta function at Qt = 
is smeared by the combined recoil of multiple emissions, each of which may add 
or subtract from the transverse momentum of the observed muon pair. ^ It is 
worth noting, though, that although the perturbative and DDT results had an 
implied kinematic cutoff at Qt = Q (and in the case of DDT resummation a similar 
restraint on the of each emitted boson), no such requirement was built into this 
final resummed result, and one must integrate over all Qt to recover the correct 
normalization. In other words, while 3-momentum conservation was observed, energy 
conservation was not. 

Imposition of such a constraint, if done smoothly, would tend to reduce the 
smearing effect, but it is not the only effect yet to be included. There remains 
the possibility (especially in QCD applications) that the incoming reactants may 
have an "intrinsic" transverse momentum, a nonperturbative effect which would 
again broaden the spectrum. Finally, there is the expectation that, due to the 
approximations which allow the above-described exponentiation to work, the correct 
shape of the distribution at higher Qt (close to Q) would be better described by 
the original perturbative result, as the logs in that region do not overwhelm the 
expansion parameter. A proper description of transverse momentum spectra will 
necessarily be a delicate balance of perturbative, nonperturbative, and resummation 
effects (see figure 3.4), and in the chapters to come these will all be explored. 

3.3 Application to QCD 

Parisi and Petronzio applied their method to a QCD example, pp fi~^fi~ + X, 
and compared the result with data. Due to the differing structure of the theory. 



^The average transverse "kick" given to the system, < Qt >, will be a useful quantity to study 
in what follows. 
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Figure 3.4. Three regions of interest. 

there are details present here which do not apply to the QED case. Some of these 
are irrelevant to their argument and will be ignored. Others constitute new features 
of the procedure, or pitfalls that would otherwise constrain its effectiveness, and these 
will be discussed. 

The structure of the resummed calculation is the same as that of the previous 
section (equation 3.16), although here the desired quantity is 



E 



(Pa 
d?QT 



resummed 



y=0 



a 



■K dQx dQdy 



y=0 
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(2vr)^ 



/ 



d% e 



ih-Qr „A(6) 



(3.20) 



where the free variables are not only the mass Q and transverse momentum Qt of 
the muon pair, but also the rapidity y (although they look only at the y = Q region). 
^ The "Sudakov" exponent A is still a transform to 6-space, but as it contains the 
effects of higher-order corrections which now depend not on the QED coupling but 
the running QCD coupling asikx), the latter cannot in general be pulled out of the 
integral. There is instead 



in which z/_|_, like the Born cross section ctq, is defined by the unresummed perturbative 
result after cancellation of the poles: 



The structure of the subprocess diagrams is the same as well, in that to leading 
order there is but fermion-antifermion annihilation, leading to an intermediate vector 
boson, which subsequently decays to the muon pair. At next-to-leading order, 
radiation from the initial state is included. The relevant diagrams have already 
been shown in the Introduction, figures 1.1, 1.2, 1.4. ^ 

Here, the fundamental interacting particles are the quarks which make up the 
colliding protons, and so parametrizations of these parton densities will necessarily 
be required. In the low-Qr approximation to which the resummation formalism 
applies, it can be shown that these partons' momentum fractions depend solely upon 
the rapidity of the muon pair and its mass relative to the center-of-mass energy of 



with a gluon from the other, will not be considered here. The absence of a soft pole in this case 
leads only to subleading terms which are either also present in qq annihilation or which the current 
formalism has no mechanism to express. These terms will be discussed later in this chapter, and in 
Chapter 4. 






(3.22) 
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the proton-proton system v^- In this center-of-mass frame, the two momentum 
fractions are 



= -^e^y , (3.23) 



" vs 

and (Jo will necessarily include sums over contributing flavors of the form 

in which Qf is the charge and qf the density of flavor /. The sum extends over flavors 
and antiflavors, and color and spin sums must, of course, also be taken into account. 

The details had also been calculated by Altarelli, et.al. (1978) [29], among others, 
and showed the expected logarithms of Q/Qt- The difference, relative to the rough 
model of the previous section, is that in practice the coefficients of these logarithms 
are particular to the subprocess under study. For Drell-Yan production these are 
2Cp and —SCp, where Cp = 4/3 is a. factor arising from the color structure involved. 
Specifically, the details are: 



z^+(fe 



3C, 



- 1 - 



2Cf[ 



(3.24) 



Before this calculation could be compared to data, a potential problem in the 
transform to 6-space needed to be circumnavigated. This transform involves an 
integral over the range < < Q"^, but the running coupling as{k^) diverges at 
the small scale A: 



as{kp) 



(33 - 2Nf) ln(/t2 /A2) 



(3.25) 
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This problem was avoided by replacing k\l with {k\ + M"^)/ A?', where M is some 
scale greater than A, but small enough that at large kx the replacement doesn't affect 
the result. 

Parisi and Petronzio then dealt with the possibility that the annihilating partons 
might carry an intrinsic transverse momentum. This would be a small effect, 
observable only at very low Qt of the pair, and would tend to "smear" the resulting 
Iow-Qt distribution further, as the additional momentum could either add to or 
subtract from that provided by the radiative effects. This is most easily incorporated 
in impact space by tacking on a multiplicative factor to the 6-space integrand: 

gA(fe) ^ ^m^-\b^<PT\..> . (3.26) 

The Gaussian width is controlled by < PT^int the average intrinsic 
assigned to the incoming partons. By damping out the high-6 region, the effect 
of the "high-frequency" components is reduced and the momentum distribution 
softened somewhat, as necessary. This is a simple example of a "nonperturbative 
parametrization" which was to be more fully developed in later years, as the following 
chapter will show. 

Utilizing the scale-dependent distribution functions of the CTEQ collaboration 
[66], and with the small values < PT^mt > = 1-2 GeV^ and = 1.25 GeV^, one 
can achieve close agreement with the data of reference [37] (see figure 3.5). Here 
S = 750GeV^, = 56GeV^, and a conversion factor C ~ 389nbarnGeV^ was 
applied to show the results in units nbGeV"^. At higher center-of-mass energies, 
nonperturbative effects are dwarfed by perturbative ones, the results approach 
independence of the phenomenological parameters, and the amount of intrinsic pt 
needed to reproduce the data is considerably less. 
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Figure 3.5. — /i^/i X resummed distribution. 
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3.4 Subleading Terms 

One final note: not all next-to-leading order terms have been included in this 
calculation. For both the gluon-radiative and Compton contributions, there exist 
remnants of collinear pole cancellations which go as (772-) times convolutions of 
parton densities with splitting functions. These, along with other finite terms, are 
not handled in the method of Parisi and Petronzio, which is focused primarily on 
consolidating the most singular terms. 

In 1984, Altarelli, Ellis, Greco, and Martinelli (AEGM) [48] revisited the 
6-space formalism with the goal of including these terms. To start, they recast the 
calculation within the framework of the subtraction method, which not only performs 
the resummation operation upon Iow-Qt approximations of the NLO terms, but then 
subtracts these approximations from the NLO terms to create additional finite pieces 
without double-counting. 

In addition to these, AEGM found a way to include the collinear remnants 
mentioned above. Specifically, after transforming to 6-space, there are terms of the 
form 



g/(a;,/i 



In 



Q 



2 

Tmax 



+ ^[g/(x,/i^)oC;(^)] 



(3.27) 



in which gj(x, /i^) is a parton density of flavor /, evaluated at momentum fraction 
X and scale /i. Pff{z) is an appropriate splitting function, and Cf{z) contains other 
finite terms. The symbol [qf o P] refers to a convolution of the parton density with 
the function P (similarly for [qj o C]): 



qf{x,fi^)oP^^-{z)=J2l^ ^g/(^,/^')i'//(^ 



(3.28) 



By defining a new factorization scale 
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In Mj ^ In + d^Qre^^'^^^ (^) ^ , (3.29) 

they were able to absorb the remaining plus-distribution terms into the parton 
densities via rescalings 

qf{x,Mf) = q,{x,^^^) + ^^\n^J: f ^qj(^ ^^')P,^{z) . (3.30) 

ZTT Li ~ Jx Z Z 

/ 

The remaining [qj o C] convolutions were left alone (and not exponentiated), 
although they too were now evaluated at the new scale M|. As it turns out, 
Mf ~ ho/h, which is the value suggested almost concurrently by the work of 
Collins, Soper, and Sterman (CSS), who used renormalization methods to solidify 
the theoretical basis of the 6-space formalism, and extend it in a systematic way to 
include subleading tiers of logarithms. We will study the CSS work in more detail in 
Chapter 4. 
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CHAPTER 4 



THE CSS FORMALISM 

4.1 Derivation 

In 1984, Collins, Soper, and Sterman (CSS) used renormalization group methods 
to arrive at an improved resummation scheme, one which not only preserves momen- 
tum conservation, but which can be made increasingly more precise as additional 
terms are calculated [49] . Although derived for Drell-Yan processes, the result is, in 
principle, generally applicable to all two-scale processes. 

It will be helpful here to recall a bit of the physical picture: at leading order, two 
partons (one from each incoming hadron) collide and produce a measured final state, 
of mass Q and transverse momentum Qt- In the case of Drell-Yan production, this 
final state is a pair of leptons. As the incoming partons are assumed to be coUinear 
with each other and the hadronic beam, if no other particles are observed, there is 
nothing for the system to recoil against, and one predicts that Qt will equal zero. 
Similarly, by energy conservation, the mass of the system of colliding partons (as 
measured in the same frame) must also have been Q, each parton contributing a 
fraction of this energy. 

In general, of course, other particles are generated from this partonic collision, 
and one can go to higher orders in perturbation theory to predict more accurately the 
QT-dependence as measured in the lab. As noted in Chapter 3, these predictions are 
plausible, without modification, at high-Qr? but in the limit Qt ^ Q, the biggest 
contributions come from cases in which the additional particles are either too soft to 
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individually detect or are themselves nearly collinear with one or the other incoming 
parton. Resummation is then required for a complete description. 

Fortunately, in this limit the long-distance (low energy) and short-distance (high 
energy) phenomena become roughly independent of one another, and one can expect 
a factorization of the cross-section as follows: 



i 

X J.(-,— ,,(/.)) J,(-,— ,,(/.)) 

X S{^M,))H{^AX9{,)) ■ (4.1) 

Here high-energy quantities such as Q and the parton energies organize 
themselves within a hard scattering factor H, while small, non-collinear transverse 
momenta krs contribute to a soft function S. Associated with each incoming 
parton (and the nearly collinear quanta k^i surrounding it) is a non-perturbative 
factor Jj, which remains dependent upon both low and high-energy quantities [kxi 
and respectively). For the transverse momenta shown, there is a delta-function 
included which relates their vector sum to the Qt of the observed system. The hard 
quantities are related to each other in a way: in the limit described, both and 
p^ are expressible as simple fractions of Q. Thus the functions Ja^b have an implicit 
dependence upon Q which will become important in what follows. Each function 
depends upon the QCD coupling strength g{fi) and its argument, the renormalization 
scale yu; the fact that the cross section itself does not depend upon /i will also be of 
import. 

Proof of such factorizations can be found, for example, in references [47, 
67]; the particular choice of variables which express the high and low-energy 
scales is different, but the general argument is the same. The crucial point 
is that factorization, coupled with invariance of the cross section under certain 
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transformations, leads to renormalization group properties of the factored functions, 
and thus to exponentiation. An all-orders resummation of the leading terms is the 
end result, the corresponding leading coefficients determined by comparison with the 
perturbatively-calculated result at fixed order. 

Equation 4.1 above is in the form of a multiple convolution; that is, only the 
matrix element, not the phase space, has factorized. Significant simplification arises 
after a Fourier transform to impact parameter space, 

J d?QTe'^-'^''^6\QT - E ^T^) = e*^' = H e*"^" , (4-2) 

i i 

in which the impact parameter is the two-dimensional axial vector h. Then, for each 
transverse momentum /ctj, one obtains 

j d'kT^e''■'^^'f{kT^) = f ip) . (4.3) 

This produces the simple product 



as(h, Q) = Ja{^, bfi, 9{f^))jb{^, bfi, g{fi))s{bf,, g{fi))H{^, ^, ^,9{f^)) , (4.4) 

where the tilde above each function is a reminder of the transform that has taken 
place. 

Taking the logarithmic derivative of both sides with respect to the large scale Q^, 
CSS arrived at 

dlnas dlnJa dlnJb d\nH 
dlnQ^ ^ dlnQ^ ^ dlnQ^ ^ dlnQ^ ' 

The logarithms of Ja and Jb now contain explicit dependences upon Q, and, as 
shown in references [47, 67-68], each derivative thereof contributes to two separate 
quantities, one a function of bfi and g{fi) only, the other a function of Q/fi and g{fi) 
only: 
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1^ + li^ = ^^^^^ + ^^^^^ • ^^-^^ 

As -ff doesn't depend on b, its derivative contributes only to a function of Q/fi and 

|j^ = GH(Q//i,^7(/i)). (4.7) 
Adding Gj and to form a new function G, one arrives at 

= gif^)) + G(g//x, . (4.8) 

CSS then made use of the fact that the cross section is independent of renormal- 
ization scale /j,: 

d dlnas d dlnas ^ .[^^^j^o, (4.9) 



d In Q'^ din jj, din fi dlnQ'^ din jj, 
and so, since K and G share only a dependence on g, it must be that there exists 
some function 7 such that 

dK dG w / N N / . N 

-A(5(/i)) . (4.10) 



dlnjj dlnjj. 

All the 6-dependence is in K, while all the Q-dependence is in G, and both can be 
independently scaled, since both satisfy their own evolution equations. 

It will be instructive here to stop and see why this is important. In a perturbative 
calculation, after dealing with the soft and coUinear poles, the cross section can be 
written as a sum of two series in as/2n: 

<^=i:[^^rK'+<^r], (4.11) 

N=0 '^'^ 

in which the coefficients are either integrably divergent as Qt — ^ (here (y^s^\ the 
"singular" piece) or are zero in this limit and thus pose no threat to the convergence 
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of the expansion (ajj , the "regular" piece). ^ The former series (u^ ) corresponds 
to the cross section we began with in equation 4.1, and we know something about 
its perturbative structure: 



^{N) 



Ti''\Q/^^,g{^^mQT) 



Qt 



m=0 

Here the "plus-distribution" [ ] , denotes a regularization such that 



(4.12) 



in"(gVg2 



Qt 



fiQr) 
or 



In" (gvg2) 



Qi 



fiQr) - /(O) 



ln'"(QVg^) 



Qt 



ln-+^(QVp^) . 



(4.13) 



m + 1 

CSS realized that, once the Fourier transform to fe-space is performed, as in 
equation 4.4, the large logs of Q/Qt become large logs of Qb, or equivalently, logs of 
Q/ fi and bfi: 



I 



d^QTe 



ln™(QVQT) 



Q 



2 

T + n=0 



m+1 



(4.14) 



Since these large logs occur, the expansion doesn't converge, and one can't even 
approximate das well by the sum of the leading log terms (m = 2N — 1). 

Fortunately, given the separability of these logarithms (as shown in equation 4.8) 
and the independent renormalization groups of K and G (equation 4.10), CSS were 
able to do something about this. ^ They scaled /i in up to order 1/6, and /i in G 
to order Q. Both the large logs of 6/i in K and the logs of QZ/i in G then tended to 0. 



-'^This structure results from the use of the subtraction method (see Section 3.4). Asymptotic 
hmits of the fixed-order result (that is, the logs of Q/Qt) arc subtracted to form c^''. and then 
added back in to form a'g^^ after pole-cancellation. 

^The derivation to be embarked upon is rather involved, and the casual reader may wish to skip 
to equation 4.35. 
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In practice, this need not be exact, and CSS allowed for variation by using ci/h and 
C2Q1 where Ci and C2 are both of order 1. The canonical choice for these constants is 
ci = 60 = 2e~^"^ ~ 1.123, C2 = 1. From equation 4.10: 



K{c,,g{c,/h))-K{h^^,g{^^)) = - '^^[g{^)) 
G{l/c,,g{c,Q))-G{Q/^^.g{^^)) = %{g{fl)) , (4.15) 



and thus 



[ir(ci, g{c,/h)) + G(l/C2, g{c2Q))] - [K{bfx, g{fx)) + GiQ/fi, g{f,))] 

rclQ^fi^^igifi)) 



Now, for any function F{b,Q), 

dF{b,Q) dF{b,Q) d\n{l/b^) 1 dF(b,Q) 



(4.16) 



d(l/62) 91n(l/62) d(l/62) (l/62)5in(l/62) 



.2\ ' 



(4.17) 



so 



/ (iF(6,Q) = F(ci/C2g,g)-F(6,Q) 

■JF{b,Q) 

d{l/¥) dF{b, Q) 



li/b^ (1/62) 5in(l/62) 
If one takes F{b,Q) = K{bfi, g{fi)) + G{Q/ fi, g{fi)), it follows that 



(4.1^ 



K{bfi) +G{Q/fi) = K{c,fi/c2Q,g{fi)) + G{Q/fi,g{fi)) 

r-lQVcl d{l/b^) d[K{bfi) + G(g//i)] 
A/62 (1/62) 51n(l/62) 

where 
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d[K{hii) + G{Q/^i)] d[K{cug{ci/b)) + Gil/c,,gic2Q))] l{g{c^/h)) 
91n(l/62) 91n(l/62) + 2 

by equation 4.16 and the theorem 



(4.20) 



d /"^ 

dx J a 

with a; = t = /i^, and a = clQ'^. Meanwhile 



f{t)dt = fix) , (4.21) 



^[K{c^,g{c^/b)) + G{l/c2,g{c2Q))] _ ,,,,, ai^(ci, (7(1/^^)) ... 

51n(l/62) -/^l^?U/'')j I4.22j 

for (3{g{x)) = and 



KMC2Q, ^(/i)) + G(g//i, ^(/i)) = Kic, g{c2Q)) + G(l/C2, gic^Q)) , (4.23) 

by equation 4.16 again, with ci/fe — ^ C2Q. 

Performing the change of variables fl = ci/b in the integral of equation 4.19: 

Jl/b^ (1/62) Jcl/b^ ^ ' ^ ^ 

CSS derived, in all. 



Kibfi,gifi)) + GiQ/fi,gifi)) = - -^A{cr,g{ji))-B{cr,C2.g{c2Q)), (4.25) 

where 



A( /-NN _ 7(^(^)) , ^-^^ ^^(Cl,^(/i)) 

^(ci,^?(^)) = -^ + 2^(^M) 
S(ci,C2,(?(c2g)) = -[Jaci,(?(c2g))+G(l/C2,(?(c2g))] . (4.26) 



Thus, returning to equation 4.8, 
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d\nasib,Q) f4Q^ da'^ ^ , . ^^ . „^ 

^if^ = - V i^-^'^""*^') - B(c„cM)) . (4^27) 

where, since K,G,'j, and (3dK/dg all have perturbative expansions in a^, so do A 

and 5. 

Treating equation 4.27 as an ordinary differential equation in C2Q with a param- 
eter 6, CSS found the solution, and generated the Sudakov exponent, as follows. 
Separating variables and integrating both sides from ci/h to C2Q gives: 



d\n.as{h,Q) = \nas{b,C2Q) -\nas{b,Ci/b) = -S{b,Q) 

ci/b 

or 

^s{b,Q) = a5(&,ci/c26)e-^('''Q) , (4.28) 



where 



/■ln(c2Q) _ p fC^Q" _ -, 

S{b,Q)^ dln(c2Q) / ^A(/i) + i?(c2Q) . (4.29) 



40' dfi^ 



Icj/b'^ /i 



This Sudakov exponent can be simplified as follows. The first term, 

/•in(c2Q) r4Q' diL^ ^, , f'iQ da fP- du', ,^ 

/ rfln C2Q / -V^ = 2 / — / —Aifi') , 4.30 

iln(ci/6) ^ Jcl/b^ fi^ Jc^/b Ji Jc,/b Ji' ' ^ ' 

followed by an integration by parts gives 



21n/i / —A(a) -2/ — In uAu 



/ci/fe /i 

^■C2Q d^ 



2\nc2Q r —A{fi)-2lnci/br ^A{fi) - 2 T ^\nfiA{fi) 

Jcl/b fl Jci/b /i Jcl/b /i 



while 
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ln{c2Q) 
ln(ci/fe) 



d\ii{c2Q)B{c2Q) 



2 hy,^ ^^^^^ 



(4.32) 



and thus 



S 



cl/b2 jj,^ 



MC2Q 



)A{fi)+Bm) 



(4.33) 



Evaluating 4.4 with Q = C1/C26 gives the coefficient of the Sudakov exponential 
in equation 4.28. Pulling out the parton distributions from ^ and rewriting the 
rest as separable functions C-^:-, CSS obtained : 



as{b,fi) = f r,—fa/Aixi,iJ.)Ca/dixyxi,b,ci/c2,giiJ.)) 

^ rl (1x2 

X J2 , fb/Bi^2,f^)Cf,r^{xyx2,b,ci/c2,g{n)) , (4.34) 



for which they, in practice, set the scale /i = Ci/b. 

To regain a useful expression in Q^-space, CSS performed the inverse Fourier 
transform on equation 4.28, added back in the finite remainders an from equation 
4.11, and arrived at their final result for the perturbative region 

^ = 1^ / d^be^'-Q^-dsib, ci/c26)e-^(''''3)| + , (4.35) 



in which 



r4Q^ drP r r^O^ . 
^ = /4l^N^)^(/^) + ^(^) (4-36) 



^ /-i dxi dx2 ^ rOsij^.^ iN), , ^ ^ xl 
^F = Y.L— L — Y.[^^] (^k\a,b,QT,Q,fi,-,-). (4.37) 

a,b ^1 ^1 ^2 ^2 N=i ^TT Xi X2 

CSS resummation allows for the resummation of tiers of logarithms two at a 
time. That is, given the ffist {A,B,C} coefficients, one can resum the ffist 2N 
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tiers of logarithms over all orders. Implicitly, there occurs a reorganization of the 
logs (L = \yi{Q/Qt)) such that successive orders grow by as opposed to a^L^: 



(J ~ 



as{L + l) 



+ at 



+ a': 



(L^ + L^) + (L + 1 



+ 



(L + 1) + + L^) + al{L'' + L') + ■ ■ ■ 

(L + 1) + a,(L-^ + L") + ^^(L^ + L'*) + . . . 
(L + 1) + a,(L=^ + L^) + + L^) + ■ ■ ■ 

+ ••• 

To obtain these coefficients, one expands the resummed form to the same order as 
one has calculated perturbatively, and compares this expansion with the asymptotic 
approximation of the perturbative result. Double-counting is avoided, as in the 
AEGM result, by use of the subtraction method. Furthermore, the C coefficients 
contain the unexponentiated subleading terms discussed in Section 3.4. 

4.2 Extension to the Nonperturbative Region. 

Strictly, the above result is valid only at low h {h <^ 1/Kqcd)- At h >~ a^t^; 
the coefficients A and B in the Sudakov exponent become dependent not only on 
as{p) but on the parton masses in the form mf/fi. A high enough impact parameter 
is reached that perturbation theory is not a valid description, and /i is then allowed 
to go small enough that both g{fi) and rrif/fi blow up. As Parisi and Petronzio have 
shown [45], as is dominated by 6 ~ a^(a^) ^ ^^^S^ Q/Qt, but 

we don't with current technology obtain large enough Q to ignore the b > 
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a. 



+ at 



a". 



region. At low b, Ci/b and hence /i stays high enougli that botli mj/fi and g{fi) are 
small, and one recovers the perturbative results above. 

It is clear that some arrangements must be made to parametrize the effects of the 
nonperturbative region, while keeping the perturbative formalism in use in its region 
of applicability. CSS suggested the following: define a function 

which goes no higher than b^ax- Use this b* for evaluations of the perturbative das, 
and rewrite the full das as 

dasib) = dasib*)[-^^^] , (4.39) 

where the second factor is approximated by a parametrization in b. 

Perhaps the easiest way to predict the form of such a parametrization is to return 
to equation 4.8 and directly integrate, starting at some minimum scale Qq at which 
one would deem a finite-order expansion in terms of as{Qo) "sufficiently accurate" 
{e.g. Qo = ci/bmax)- Dividing by the result taken a.t b = b* yields a ratio of the form 
[69] 



dasib) 

exp 



^^^^^^^ ^-hM HQ'/Ql) - hAixl b) - hsixl b)] . (4.40) 

We will hereafter refer to the exponent above as the nonperturbative function 
Snp- Like the parton distributions, the coefficient functions h are intended to be 
universal (generally applicable) and extractible from data. The one constraint is 
that exp{Sj\fp) ^ 1 as 6 ^ 6* ^ 0. 

To date, sufficient data have not been taken to study the flavor dependence of the 
functions hA,B, nor the dependence upon x° or x° individually. However, attempts 
have been made to fit simplified versions of the above form. Ladinsky and Yuan (LY) 
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[70] in 1994 and Landry, Brock, Ladinsky, and Yuan (LBLY) [71] in 1999 used the 
three-parameter form shown here: 



Snp = -g2b'\n(^) - g,b^ - g^g■sb\nimxy,) . (4.41) 

A previous analysis by Davies, Webber, and Stirhng (DWS) [72] in 1985 was 
made with a two-parameter form obtained by setting g^ = in the above. LBLY 
also studied the two-parameter form, and the resuhs of all these efforts are collected 
in Table 4.1 below. In Chapter 7, we will provide evidence that a single-parameter, 
Q-independent form may agree better with data. 



Table 4.1. Nonperturbative Parameters. 





DWS 


LBLY 


LY 


LBLY 


giiGeV) 


0.40 


0.24 


0.11 


0.15 


g2{GeV^) 


0.15 


0.34 


0.58 


0.48 


gs{GeV-^) 


N/A 


N/A 


-1.5 


-0.58 


QoiGeV) 


2.0 


1.6 


1.6 


1.6 


bmaAGeV-') 


0.5 


0.5 


0.5 


0.5 


PDF 


DOl 


CTEQ3M 


CTEQ2M 


CTEQ3M 


DATA 


E288,R209 


E288,R209, 


E288,R209, 


E288,R209, 






CDF-Z,E605 


CDF-Z 


CDF-Z,E605 



In lieu of equation 4.35 then, the final CSS result can be written as: 

rf(7 = 1^ J d'^y^-'^^'^dasib*, ci/c26*)e-^('*'^)e^'^^(''^)| + dap ■ (4.42) 

4.3 The Sudakov Exponent at very low h. 

Now, what about the very low-6 region, where ci/b > Here there are no large 
logs, but one does not want the Sudakov exponent to change sign. In practice one may 
consider cutting off the exponent, that is, simply taking exp[— S'(6)] exp[0] = 1 
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for this region. But will this preserve the proper normalization upon integration 
over Qt? Altarelli, et. al. (1984) [48] calculated the exact first-order result for the 
Sudakov form factor and arrived at 



S{h,Q) 



27r 



(4.43) 



which by inspection has the desired property of becoming zero as 6 — 0. 

In the proposal of Ellis, et. al. (1997) [73], the CSS form for S is maintained, 
but new scales \{h) and /i(fc), which never go above C2Q, are used as lower limits: 



S{h,Q) 



To agree with the Altarelli result, these scales must be defined such that 



•^aQ dx , C2Q 
— In 

\ X X 

rc2Q dx 

X 



/■CSV ax 

J a X 



1 i„2 caQ 

2 A 



In 



dx , 
— In 1 - Jo (to) 

X X 
caQ dx , 



r '^[i-j,(hx)] 

Jo X 



and so must be 



(4.44) 



(4.45) 



M6) = Qexpf-io — (l-'^o(M)} 

X{b) = Qexp{- / —ln{Q/x){l-Jo{bx)) }. 

At large 6, A and fi both ^o/^? in agreement with CSS. 



(4.46) 
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CHAPTER 5 



KT-SPACE RESUMMATION 



As successful as the CSS formalism has been, it is not necessarily easy to put into 
practice. Performing Fourier transforms numerically is a slow process, requiring many 
samplings of the parton distributions to achieve a prediction for a single Qt value. 
The impact parameter integral formally extends to 6 = oo, but of course one must in 
practice halt the integration at some hmaxi and either drop the rest or introduce an 
asymptotic expansion [74]. Moreover, the differing treatment of asymptotic pieces 
between resummed and subtracted versions can affect their cancellation, leading to 
unphysical behavior at large-Qr- One must take pains to smoothly yield to the 
ordinary perturbative result in this region. It would be nice if one could perform the 
Fourier transform analytically, or avoid it altogether. 

A number of authors have undertaken this task, building on the foundation laid by 
Dokshitzer, D'yakonov and Troyan (DDT) [44, 75]. Their central result is essentially 
that of equation 3.12, in a more general two-scale form: 



5q55qp^ = 5^ Esqi <^.)^■.w^<^^)A(-..<^.)^■.w^ 0^)1 . (s-d 

where the momentum fractions and Sudakov exponents are 



b 



b 



Q_ 

exp 



2n 
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b y b 



(5.2) 
(5.3) 



5.1 Kimber, Martin, and Ryskin 

In 1999, Kimber, Martin, and Ryskin [76] revisited the theory behind the 
DDT equation, and reinterpreted it in terms of /cr-dependent parton distributions 
f{x, kT,Q), e.g.: 



d 

fa{Xa, ka, Q) = -^[fa{Xa, A)Ta(A, Q)) 



(5.4) 

X=ka 



That is, the perturbative terms encompassed in the Sudakov exponential Ta are 
combined with standard distributions fa{xa,X) to form new, unintegmted parton 
distributions, describing partons with transverse momentum ka at samphng energy 

ka<Q. 

The DDT form is then obtained by apphcation of the strong- ordering condition 
to the convolution over fa and /{,. The Qr-dependence being now described by a 
single luminosity function Cah'- 

dQHQldy = Wd-y^, ^^'^""^ ""'^ ' ^^'^^ 

with 

r ^ ^ ^ d'^ka d'^kh 1 

^ab = / fa{Xa, ka, Q)fb{Xb, kb, Q)5'^{ka + h - Qt) ,n , (5.6) 

J ki ki TT 

application of strong-ordering allows only one of the off-shell partons to contribute 
at a time. That is, either ka <^ kb — Qt or kb <^ ka — Qt- In the first case, the 
delta function collapses to S'^{kb — QT), while in the second case one gets 6'^{ka — QT)- 
Summing over both possibilities yields 



^abiXa,Xb,QT,Q) = -^^\faiXa,X)Ta{X,Q)fbiXb,X)TbiX,Q)]^^ , (5.7) 

which reproduces equation 5.1. 
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Note that since Ta^b{xa^h,Q,Q) = 1, it follows that integrating the luminosity up 
to scale Q recovers the standard luminosity : 



fQ^ 2 

/ dQrpCab{Xa, Xb, Qt, Q) = fa{Xa, Q) fbi^b, Q) , (5. 

Jo 



and in particular, 



k 



^fi{xi, ki, Q) = fi{xi, Q) . (5.9) 



5.2 Ellis and Veseli 

In 1997, R.K. Ellis and S. Veseli (EV) came up with a similar extension of the 
DDT result, while maintaining the improvements of the CSS formalism [77]. They 
did so by successfully performing the CSS 6-space Fourier transform analytically, for 
the first three tiers of logarithms. Starting with the purely perturbative CSS result 



in which fa^b are the convolutions over finite, unresummed C^ji corrections 

/a(a;a,x) = H I —fa/A{Xa/z,bo/b)Ca/a{z,bo/b) 

- JXa Z 

fb{xb,^) = Y.[—fb/Bi^b/z,bo/b)Cy-^{z,bo/b), (5.11) 

~ J Xf, Z 

b 

EV defined a fixed scale ratio t = Q'^/S and used the rapidity y to write 



= I dXa I dXbS{t - XaXb) ^^^^'^j^j!^f^^^ . (5.12) 

JO Jo atayaQj. 
Then, after substitution of equation 5.10, they took the A^th moment with respect 
to t: 
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dt t 



N 



da 



dtdQ\ 



-S{b,Q) 



= E?r^^ bMQTb)e' 

~t 2 Jo b b 



The modified parton density moments fi evolve according to 



(5.13) 



l'NfiiN,fi), or 



din fi"^ 

MN,bo/b) = /.(iV,Q)exp 



(5.14) 



in which the are as(/i)-dependent anomalous dimensions. After defining x = Qrb, 
this allowed EV to reorganize the cross-section moment as 



E(iV) = I - dx xJo(x)e"-(^'«-«) 



2 JO 
where H^iQ) = Ea,b fa{N,Q)fb{N,Q) and 



(5.15) 



Un{x,Qt,Q) 



blQl/x-2 jl 



00 n+1 

EE 

n=l m=0 



^ 'Aln^ + fi + 27V 



2tt 



2 u2 ■ 





Qlb 



(5.16) 



That is, the new parameters nDm are directly calculable from the CSS parameters 
{Ai.Bi}. With the relation 



d_ 

dx 



[xJi{x)] = xJq{x) 



(5.17) 



equation 5.15 can now be integrated by parts, with the result 
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(Xo Hn{Q) 

2 Qjn 



xJi(x)e"^ 



dx xJiix) — ; 

dx 



(5.18) 



Due to the rapid damping of exp(WAr) and Ji( 
boundary term; and, using 



oo, one can ignore the 



X 



X 



Q'^dx bdQr 
exhibit the similarity to the DDT formula: 



dQl ' 



(5.19) 



(5.20) 



In fact, were it not for the logs of x/Bq in Un, one would have the standard Sudakov 
exponent, containing only logs of Q'^/Q'^, which are, of course, the important ones. 

For this reason, EV split these off and defined a remainder TZn- First, they 
defined an integrand TZn via 



Un{x,Qt,Q) 



oo n+1 

n=l m=0 
oo n+1 



27r, 
2^, 



In'" ^ D 



exp[E E 

n=l m=0 

WAr(6o, Qr, <5) + ilN{x, Qt, Q) 



ill ^2 nJ-^m 



+ TZn{x,Qt,Q) 



(5.21) 



where 



oo n+1 

nN{x,QT,Q) = exp[^ ^ 

n=l m=0 
oo n+1 

- exp[^ 

n=l m=0 



2^ 
27r 



n D^T-^ 

/02 1,2 " ™ 



(5.22) 



Then, since Z^a? no longer depends on x, and /q°° dxJi{x) = 1, equation 5.20 becomes 



s(iv) = -^{^fo//^(g)e"-('""«-«) + 7i^(gT,g)} 



(5.23) 



where 
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TIn{Qt,Q) = (foHNiQ) / dxJi{x)nN{x,QT,Q) ■ 

Jo 

Using equations 5.14 and 5.16, the Sudakov exponent is defined: 



(5.24) 



-5(fe-»bo/QT,Q) 



V[/ dXaX^ fa{Xa,QT)]\ f dXbX^ fb{Xb, 



-S(bo/QT,Q) 



(5.25) 



at which point the inverse transformation back to Q-space becomes easy. Comparison 
of equations 5.25 and 5.13 leads to Elhs and Veseh's final result in the perturbative 
regime 



da 



dQ^dydQl 



d 



dQ\ 



fa(,Xa,QT)fb{Xb,QT)e 



-S{bo/QT,Q)' 



dQl 



(5.26) 

where S^^ denotes the inverse transform N Q, and the Sudakov exponent becomes 
the same as in the 6-space formalism (equation 4.36), except that the lower limit feg/^^ 
is replaced by Q"^. By using known integrals of Ji{x) In"^ ^, Ellis and Veseli were 
able to prove that the remainder term TZ contributes no more importantly than the 
NNNL series of logarithms, that is, those corresponding to the CSS B^'^^a'^ terms. 

One still requires a prescription for dealing with the very lowest Qt values, for 
which perturbation theory doesn't hold. As in the CSS paper, this is handled by 
including a non-perturbative function -F/vp, this time of Qt, while sampling the 
Sudakov exponent and parton distributions at a scale which never goes below a 
limiting value QTHm- To maintain proper normalization, while affecting only the 
small-Qy region, Ellis and Veseli chose the following parametrization: 



Fnp{Qt) = l-e-^'^T^ 
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(5.27) 



Qt 



QTUrn 



The resummed piece then becomes 
das - d 



dQ'^dydQ^ 



E 



T 



dQ^ 



T 



where 



(5.28) 



(5.29) 



S{Qt,Q) 



dfl^ r Q 



(5.30) 



This can be used with the same perturbative remainder ap as in equation 4.37. 
Transforming to 6-space is avoided, drastically improving the speed of numerical 
implementations, although the method remains susceptible to pathologies at high- 
Qti especially at fixed-target energies, as explained further in Section 6.6. 
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CHAPTER 6 
DIRECT PHOTON PRODUCTION 

Here we begin the process of putting together a consistent, resummed, next-to- 
leading order (NLO) description of single photon hadroproduction. ^ However, the 
resummation prescriptions described in Chapters [3,4,5] all describe how to reorganize 
logarithms of Q'^/Q't-, where Q and Qt are properties of the kinematics of a pair 
of observed particles. In predicting the cross section for inclusive direct photon 
production, we are, in the end, concerned only with the properties of the photon, so 
if we want to include the effects of resummation we must, as an intermediate step, 
define a second potentially observed particle for application of these prescriptions. 
Fortunately, kinematics requires there to be a second final-state particle at leading 
order anyway. To NLO, the required Feynman diagrams involve up to two additional 
particles in the final-state (besides the photon). At the subprocess level these are 
either quarks or gluons; by the time they reach the detector they have fragmented 
into hadronic jets, so named because each has very little transverse spread about its 
original partonic vector. That is, the hadrons into which each parton fragments are 
concentrated within a narrow "cone" about the parent parton's direction. 

If one is ultimately unconcerned with precisely what kinds of hadrons are 
produced (as in our case), then even if one could identify each type within a particular 
jet, one would be ignoring (summing over) this information, so in our calculation we 

need only keep track of the kinematics of the jet itself (or, equivalently, the parent 
^AU relevant diagrams and matrix elements are collected in Appendix A. 



81 



parton); the fractional probabilities of particular fragmentations all sum to one by 
definition. 

That said, we are left with the following question: our NLO calculation consists 
not just of 2 — > 2-body leading order (LO) contributions (for which assignment of 
the jet is unambiguous), but also 2 3-body contributions involving a photon and 
two potential jets. To which do we pair with the photon and apply the resummation 
formalism? 

The answer is both, in turn, but in the following way. Experimentally, if the 
jets are far enough apart in direction as to be distinguishable, we can define the 
"pair" to be the photon plus the jet of greater transverse momentum. Given that 
the experiment cannot tell which parton initiated that jet, however, our calculation 
must consider both possibilities. And, in the instance in which the jets are not 
distinguishable, our calculation must treat the resulting kinematics as that of a single 
jet. Each experiment has its own criteria for distinguishability, the details of which 
are important only if the final predictions include explicit jet characteristics. As this 
will not be of concern to us, our calculation can consider any nonzero separation 
distinguishable. 

Thus, in addition to subprocess matrix elements and phase-space factors dT 
(both of which are given in terms of subprocess momenta {pf}), the ingredients of 
our calculation should also include a jet definition, that is, a set of delta-functions 
the job of which is to define the jet variables pj and pair variables = p^i^ + pj in 
terms of these subprocess momenta. 

The photon (p^), of course, can be unambiguously associated with one of 
the outgoing subprocess momenta {^1,^2)^3} (we'll choose Pi). However, after 
application of momentum conservation, many of the variables included in {^1,^27^3} 
will not be free but expressible only in terms of other variables in the set. Depending 
on the contribution being calculated, the most convenient remaining set for evaluation 
of that contribution may not include pi, and thus carrying around explicit definitions 
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for both the photon and jet would help keep track of their properties in terms of the 
free subprocess variables. In the end, integration over the phase space dT then yields 
a cross section in terms of the desired jet and/or pair observables. For the sake of 
clarity we will suppress these delta-functions in what follows. 

The result of this chapter will be a collection of singularity-free pieces which, when 
added, give an unresummed, NLO cross section for 7 + jet production, but which 
are also organized such that their assignment to CSS-style resummation coefficients 
is self-evident. It is then a simple matter to use these coefficients in a resummed 
calculation, and integrate over the jet numerically. Indeed a FORTRAN program, 
switchable between NLO and resummed output, has been written concurrently with 
this dissertation, and is available upon request. 

The attempt is made throughout to keep as many free variables as possible, 
consistent with the requirements of divergence cancellation and infrared safety. This 
is required in order to account for kinematic recoil of the photon-jet system against 
the soft radiation. 

A word on notation: the 4-momentum of particle i will be represented by 
throughout. The energy of such a particle is while the magnitude and direction 
of 3-momentum are \pi\ and pi, respectively. We use massless particles exclusively, 
for which \pi\ = p^. The transverse momentum, |pj|sin6'j will hereafter be noted 
as simply pi, while the rapidity is yi. If we need to refer to the vector transverse 
momentum of particle i, we will denote this as pTi- Dot-products of 4- vectors, 
usually represented using Greek indices and the Einstein summation convention, 
as in PiPjfj,, will for the sake of clarity be shown as pi ■ pj. Both describe the quantity 
PiPj — pjp] — pIp] — p^p^j- The following table, valid for massless particles, should 
help to clarify our notation: 



Pi 



\Pi 



Pi cosh Ui 
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p • = \pi I sin 9i cos = Pi cos 
p,^ = |pi I sin 9i sin 0^ = sin 0^ 

=|pi|cos6'i =pisinhj/j. (6.1) 

6.1 Two-Body Final States 

The two-body phase space in arbitrary dimension D = A — 2e is 

^"r^ = , , C% i2nf6-ip^^ +Pi-P^.- f.) , (6.2) 

2py(27r) 2p^(27r) 

where the incoming parton vectors are, in the hadronic CM., given in terms 

of the momentum fractions {xa,Xb} as follows: 



= a:,^(l,0,...,0,l) 

Pb = x„^(l,0,...,0,-l). (6.3) 

For a given matrix element the corresponding 2-body cross section is written 
as an integral over these momentum fractions, weighted by the squared matrix 
element (suitably summed and averaged over final and initial spins and colors), 
the phase-space factor, and a luminosity function C which gives the probabilities 
of finding incoming partons with those momenta inside the parent hadrons. Thus 
the cross section for subprocess i is written: 

d^'ai = dxa C dxbd''T2Ci{xa, Xb)Y.\M I J , (6.4) 
Jo Jo 

in which the luminosity is a function of the relevant parton densities {qf, qf} inside 
hadrons A and B: 

J:fjqf{xa)qf{x,) 
Ci{Xa,Xb) = ^ . (6.5) 

ZXaX}}0 

The indices / and / each run over all flavors of partons relevant to subprocess 
i (here the gluon, as well as all quarks and antiquarks, are considered "flavors"). 

84 



Thus, e.g., for the qg — > 7g subprocess the sum is simply J2f Qf{xa)Qg{xb), in which 
/ ranges over quark flavors and antiflavors. We will sometimes refer to this sum, for 
a given subprocess i, as Hi{xa, Xb). 

It should be noted that charge factors, denoted in this work by Qf, are present 
at photon/quark vertices, and, as they depend on the quark flavor /, need to be 
included in the sum. Each gives the value of the relevant quark's electromagnetic 
charge in units of the proton charge (+e). As they arise in the matrix element, each 
above should properly have a subscript //, and the summation symbol X)// 
should be placed outside C. For clarity, we take these abbreviations to be understood 
in what follows. Actually, we will no longer refer to the matrix element as a whole; 
each can be written as the product of factors 

T.\^\' = ^2,[K'T'^f{i%, (6.6) 

in which C1J27 = 2(47r)^aasCi?/i^'' is independent of the particular subprocess and 
K'T'j,j{v) gives the angular dependence as a function of v = e^Y(e^^ + e^^). This, 
of course, is dependent upon the particular subprocess, and upon the quark charges 
involved. The "primes" denote e-dependence; when e is taken to zero, the primes 
will disappear. 

We may now begin to simplify matters by integrating over some of the unneeded 
variables in d^r2. Breaking down the pi factor as 



X d^~^pTi6^~\pTi - Xt) 

= -^^{^a - Xa)5{Xb - Xb)d'^~'^PTl5'^~'^{pTl " ^t) , 

(6.7) 

where 
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X, = ^{e~y^ + e~y-) , (6. 

we can integrate over Xa,Xb, and pti, which has the effect that d^r2 becomes 



d ^2= o . n.. ,D-i ■ (6-9) 



Since we'll require (at most) a transverse momentum, rapidity, and azimuthal 
angle for each particle, we can further reduce our variables as follows. Again, 
note that here and throughout this dissertation, we will abbreviate single-particle 
transverse momenta as = \pi \ sm9i : 



jd\p2\sm''^'e2de2d''-'(f)2 



P2'^dp2dy2d^~^(j)2 
Of course, we could stop here, and leave our result in the following form: 



(6.10) 



d'^ai = £(x„ Xb)u;2^[K'T'^jiv)]d''T2 , (6.11) 

where 

,Dr dyidp2dy2d'^^^(j)2P2^'^ fr, , o^ 

However, to use this result for prediction of photon cross sections, it will be helpful 
to know what the associated photon variables {p^, y^, (p-y} are in terms of the above. 
Since we are assigning the photon to Pi and the jet to P2, we have the following: 



P7 = P2 Pj = P2 
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2/7 = 2/1 Vj = 2/2 



07 = 02 + Vr 0i = 02 • 



(6.13) 



For a leading-order (LO) calculation, this is all that would be necessary, and the 
reader could simply plug in the ah — ^ '-jd matrix elements given in Appendix A. ^ 

For a NLO calculation, however, in which the panorama of 2-body contributions 
includes virtual and counterterm pieces, as well as subtraction terms derived from 
3-body asymptotic approximations, all with poles in e, the constraint of infrared 
safety forces us to choose some common set of free observable quantities, in order 
that addition of these pieces will result in pole cancellation. We have four so far: 
{y^,Pj,yj,(f)j}, but in order to fully define two independent massless 4-vectors 
(as our formalism requires), we need six altogether. Specifically, in order to use 
the resummation procedures described earlier in this work, we will need the pair 
transverse momentum Qt', in order to incorporate smearing of the photon kinematics 
due to recoil, we will also need an additional angle, which we take to be the axial 
direction 0^ of the pair. 

For the two-body contributions, of course, Qt is zero, and 0^ is undefined. We 
will consequently tack on a factor 



for these pieces, to supply the additional two degrees of freedom required. 

Furthermore, there will be certain multiplicative factors common to all pieces, 
which will be convenient to express explicitly here. To that end, for our ab — > 'yd 
subprocesses, we take D = A — 2e, and define the following quantities: 



^QCD Brcmsstrahlung diagrams also contribute to order aas; these will be discussed at the end 




dy^dpjdyjd^ '^"(pjdQxd^ ^''05 



of the section. 
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through which one can show: 



daf°'^ = C{xa,Xb)u;2-y[K'T'^fiv)]4''r2 

^ debjKSiQx). (6.15) 

For a resummed result, the pair mass Q also requires definition; for these pieces 
it is Q = y^XaXhS = ^2pj{l + cosh(yi — 1/2)), as can easily be verified via the pair 
definition q^^ = pi^ + pj. 

Two-body Bremsstrahlung also contributes; here one final-state parton in a 
purely QCD subprocess (such as qq — > qq) fragments into hadronic material and 
an essentially coUinear photon. The probability for this to occur is thus most 
conveniently described by photon fragmentation functions D^/q{z) and D^/g{z), each 
a function of the fraction z of the parent parton's momentum that is carried away by 
the photon. The transverse momentum imbalance between the photon and away-side 
jet in this case is not zero, but instead given by Qt = p-y — pj = {1 — z)pj. 

There are 11 such subprocesses to consider (see Appendix A), and the cross 
section for subprocess i, in terms of the definitions 6.14, is: 

daf"'^ = 4[7J>.- / dzSiQr - (1 - z)p^)D,/,iz) 

= d^[yj]Ki ——D^/J^ -]0{{l - z^in)Pj - Qt) , (6.16) 

a pj " \ pj J 

in which D^/c{z) is the fragmentation function relevant to parton c in the subprocess 
{ab cd). Here a lower bound Zmin is included for the general case in which minimum 
photon energy constraints {isolation cuts) are experimentally imposed [50]. 

One would expect such processes to be of order act^; in fact, the D^/c are of 
order a/ as, and so one obtains overall 0{aas) quantities. We use the leading-log 
fragmentation functions of reference [78] : 



LL f 



F{fi) rg^(2.21 - 1.28^+1.2922)20-049 

z [ 1 - 1.63 ln(l - z) 

F(/i) 0.194 ^^ ~)i-°V0.97^ 



+ 0.0020(1 -2)'-Vi-5^ 

(6.17) 



where F(2,/i) = (a/27r) Inifi^A^) with A = 0.2GeV. 

At NLO, we will require first-order expansions of the parton distribution and 
fragmentation functions, in order to cancel coUinear singularities which arise in 



three-body final states (see Section 2.4, in particular equation 2.52). In the MS 
scheme, with common factorization scale Mf, these counterterms are as follows: 



where 



a 



P~^(z) 



dz H{Xa/z,Xb) Pa/di^) 



Xa Z H{Xa,Xb) 



P^Az) 



1 dz H{Xa,Xb/z ) ^b/b 
X, Z H{Xa,Xb) 



27r V M] J Til - 2e) 



(6.18) 
(6.19) 
(6.20) 



(6.21) 



In addition, there will be virtual diagrams, which contribute at NLO in the form 



dar'' = d,Ylj]K,f3,6{QT) 



Al A 



-if 



(6.22) 



The virtual diagrams are shown in Appendix A; their parameters A^, Al, A^, are 
given in Appendix C. 

6.2 Three-Body Final States 

As a general guide we use the work of Ellis, Kunszt, and Soper [79-80]. 
These authors arrived at a general subtraction algorithm for calculating infrared-safe 
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quantities in hadron-hadron collisions. The process is as follows: We begin with the 
3-body matrix elements for NLO single-photon production, as listed in Appendix 
A. These will have singularities in those regions of phase-space for which one of the 
final-state particles become soft or collinear with another particle. By the process 
of partial-fractioning, we can break these matrix elements into smaller pieces, each 
of which has at most two singular regions, one soft and one collinear. The potential 
singularities manifest themselves as quantities in the denominator which tend to 
zero in the singular regions of phase-space, but as the numerators simplify in these 
regions also, we can define approximate, simplified numerators which can be added 
and subtracted in all regions. Each partial-fractioned 3-body piece pairs with a 
subtracted approximate piece to form a function which is finite everywhere, while 
the leftover (added) approximations are dealt with separately. This is the essence of 
the subtraction method, as contrasted with the Phase-Space-Slicing (PSS) method, 
in which the singular regions of phase space are cut out and their effects added in 
later after simplification and analytical integration. A NLO inclusive single photon 
calculation using the PSS method, but without transverse-momentum resummation, 
was performed by Baer, Ohnemus, and Owens (1990) [56]. 

The leftover "subtraction" pieces contain all the singularities of the original 
3-body piece, but are simpler to work with, as many angular dependences fall away. 
What is required is that the singularities be made explicit so they may be cancelled 
against virtual pieces or absorbed into distribution and fragmentation functions. It 
should not matter, for a given singular quantity, whether this end is achieved via 
integration in Z)-dimensions or via extraction in terms of a pole and accompanying 
plus-distribution. In the present calculation, we use the former method for all but the 
required six free variables, and handle poles in these variables with plus-distributions. 

Once the Born terms and virtual contributions are added, and the poles taken 
care of, what remains is a completely finite NLO result. From this, one may 



90 



obtain parameters for use in a resummed version of the calculation by expanding 
the resummed form to NLO and comparing coefficients. 

Crucial to the success of any such venture is the recognition that one of the 
vectors we intend to observe, namely that of the jet, is not uniquely defined in terms 
of the vectors we begin with (call them {^1,^25^3})- We may declare that Pi be 
assigned to the photon, but it is an experimental reality that any jet we observe 
may correspond to either or both of the final-state partons emitted in our 3-body 
subprocesses and P3). Any attempt to force a one-to-one correspondence will 
lead to uncancelled poles, which is another way of saying that the result wouldn't be 
infrared- safe. One must define the jet kinematics in such a way that if the final-state 
partons are together within a certain angular region, this is counted as a single jet 
(see figure 6.1). One practical consequence of this is that singularity cancellation can 
only be done after the free variables have been reduced to those of the photon and 
jet, not those of {pi ^'^^ 

- OUT - 




Figure 6.1. p^ inside jet cone. 

6.2.1 Partial-fraction Three-Body matrix elements 

Three-body phase space, in analogy with equation 6.2, is: 
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2rf(27r)^-^2pO(27r)^-^2pO(27r)^-^ 
X (27rA^(p;^+pC-p^-p^-p|;), (6.23) 

while the squared, summed, and averaged matrix element (or partial-fractioned 
piece thereof) is denoted by ^^id is a function of all five momenta 

As written, the final-state 4-vectors are allowed to range over all space, and there 
are then up to 12 potential singularities: 

1. p'i, soft, 

2. Pi, P2, Ps coUinear with the beam in either direction, and 

3. Pi ■ P2 = 0, pi ■ Pa = 0, p2 ■ Pa = coUinear singularities. 

If we assign pi to the photon, and through cuts make sure it is never soft or 
coUinear to the beam, this reduces us to 9 singularities. If we go further and declare 
P3 to be that parton with the lowest transverse momentum (written as simply ps), 
then the P2 singularities disappear, providing we symmetrize our matrix elements 
with respect to P2 ^ P3 • We will consequently define a symmetrized squared matrix 
element (note the tilde instead of a bar) 

EI-MI' = EI-^I'(p^:>P^Pi>P2,p!;) +T.\-Mf{p'a,P',,P'uP'3,P^) , (6.24) 

and introduce a theta function 6{p2 — Ps). 

We are left with 5 singular regions, corresponding to P3 being either soft or 
coUinear with one of the other 4 particles. Since momentum conservation prohibits 
more than two of these from occurring at once (for nonzero photon p^), the divergent 
terms with which we must deal will be at worst of the form 
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1 

Pm ■ PsPn ■ PZ ' 

This arrangement is singular in two different "coUinear regions" . To simplify to only 
one coUinear possibility, we may partial-fraction these terms as follows: 



1 1 



1 1 



(6.25) 

Pm ■ PsPn -PS PS- [Pm + Pn) " " " 

The factor [p3 ■ (p^ + Pn)\~^ gives a singularity only when P3 is soft, but if this turns 
out to be the case, then the simpler quantity will be zero as well, and within an 
already partial-fractioned term, it will be simpler to express the pole with ^3. 

So, in all, we can split each (symmetrized, summed and averaged) matrix element 
into four terms, each a product of an explicit [pzPm'Pz\~^ singular factor and a 
general function of {^1,^25^3}- 

£|A<lL-.3= E I-^. (6.26) 

m=a,b,l,2 P^ P^ 

in which, by virtue of equation 6.25, itself (and thus the cross section) can be 
written as a sum of terms 

^n.- E .'^'"l (6.27) 

n=t^l,2 ■ (Pm + Pn) 

Refer to Appendix A for a list of the 3-body matrix elements which contribute. The 
"^mn, however, will not be listed. 

6.2.2 Reduce Three-body Phase Space. 

The structure of the cross section is analogous to the two-body case: 



d'^as = iyxi£ dx2C{xi,X2)^\M\]jd''T3e{p2-p3) , (6.28) 
with the phase-space factor 
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2p?(27r)^-^2p0(27r)^-^2j90(27r)^-^^ ' ^^'^ ' 

but after the manipulations given in equations 6.7 and 6.10, we now find 



XI = ^(pie^^i +p2e^^" +^36^^=*) . (6.31) 

in whicfi 



Vl = \lvl+Vl + 2P2P3 cos (02 - 03) 

0, = t,,-ifP^^i^'^2+P3sin03\ ^g32) 

VP2 COS02 +P3 COS 03/ 

We must now arrange our cross-section ingredients in a form which makes the 
function to be approximated exphcit. We'll call this function J-'mn, corresponding to 
the partial-fractioned matrix element numerators defined in section 6.2.1. The phase 
space and singular factors will not be a part of this function, while the luminosity 
factor and \E'm„ will. By analogy with equations 6.6 and 6.14, we'll also find it 
convenient to extract the overall factor uz = 4(47r)^Ci;'aa^/i^'^ common to all our 
3-body matrix elements, and define the new quantities 



(ie[123] = dyidp2dy2d(j)2 '^dp^dysacp^ 



Xe) ---- 



by which our partial-fractioned cross-section pieces can then be expressed as 

,D 4[123]x(e) ^ . 

d (Jmn = ^mnO{p2-P3) , (6.34) 

P3 Pm ■ P3 

with 
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•^mn(Pl,P2,P3) = >C(Xi,X2) ^^^^^ T • (6.35) 

VZ ■ [Pm + Pn) 

Keep in mind, of course, that we will often refer to dam and which are, by 
equation 6.27, sums over damn and J-'mn, respectively. 

Strictly speaking, the above is used only for the m = {a, b, 2} contributions. 
The 771=1 contribution contains poles in the limit pi, and we'd like to take 

advantage of the math already done for the p^ p2 case. Thus we begin by switching 
1 ^ 2 in the phase space element dV^ only (after approximation, we will switch back 
in the subtracted piece, in order to use the same Monte-Carlo event generator). Thus, 
in place of equations 6.33, we wind up with 



(i,[213] = dy2dpidyid^ "^"(pidp^dysd^ ^'03 



in which xi and X2 retain their usual definitions, but this time pi and 0i are fixed, 
while P2 is a function of {y2,Pi,P3}- 



-<P3 



P2 = \Ipi+pI + 2piP3 COS (01 

_irpiSin0i +p3Sin03n 

02 = tan ■ — + TT . (6.37) 

Lpi cos 01 +P3 COS 03 J 

6.2.3 Define and Subtract Approximations 

We'll now define the additions and subtractions of approximate versions of these 
functions J-'m- Each multiplies a factor [p^ Pm 'Ps]"^ which becomes singular in the 
regions in which either disappears altogether {p^ ^ 0, the soft region) or becomes 
collinear with particle m (the collinear region, pm ■ Ps 1). The essence of the 
subtraction method is, as previously stated, to subtract off a simpler asymptotic 
approximation which retains the same pole structure. This "subtracted piece" (or 
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asymptotic piece) pairs with the original three-body partial-fractioned piece to form 
a function rfcr^*"**^, finite everywhere in phase space, and which goes to zero as the 
singular regions are approached. 

This process then, of course, requires us to add back in that which was subtracted; 
so our singularities re-enter our calculation, but now in a much simplified form. These 
"subtraction pieces" can now be analytically manipulated in order to extract their 
poles explicitly, these being subsequently cancelled by the addition of virtual and 
counterterm corrections. 

In practice, since two types of poles exist in each original partial-fractioned piece, 
it is usually convenient to allow for two separate subtractions in order that each 
resulting subtraction piece contain only one type of pole. There is a further reason: as 
stated in Section 2.4, we expect the soft limits to reduce to known 2-body subprocess 
matrix elements, and the coUinear limits to show up as these times convolutions 
of familiar splitting functions and luminosity factors. Thus for the m = {a, b} 
(initial-state) terms, it will be convenient to have one version which fixes the direction 
of P3 while taking to zero (the soft approximation), and another which takes 
Ps — > Pm (coUinear). Then, to keep our soft singularities in one place and avoid 
double counting, we subtract the soft limit from our coUinear piece. 

The above is only a rough statement for the moment, however, as we did not leave 
ourselves with p^ and ps (the energy and direction) as free variables, but instead the 
set {p3,2/3,03}. We must be more precise. 

First of all, we eventually intend to integrate over y^, so we must be clear about 
the limits of this integral. Allowing y^ to go to infinity in either direction, without a 
corresponding decrease in the magnitude {p^l, would at some point violate momentum 
conservation, as the momentum fractions of the incoming partons cannot exceed 1. 
In other words, when we take initial-state soft and coUinear approximations (in which 
the transverse momentum p^ necessarily tends to zero), we can do so only if these 
factors of p-^ stand alone; if p^, multiplies some divergent function of y^, we must take 
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additional care. So far, we have not explicitly included these bounds, but we do so 
now in terms of restrictions on y^. 

From equations 6.31 and 6.32, we have the momentum fractions 



V o 

where 



Pi = \IV2+Vl + 2p2P3 cos(02 - 03) • (6-38) 
Restricting xi < 1 is thus equivalent to the restrictions 



< In 



- Pie 



P3 

VS - pie^y^ - p2e 



-y2 



ys > (6.39) 

P3 

However, since the factors of p^ in pi do not multiply exponentials of y^, these can 
be dropped in the low-ps limit, leaving 



{l-Xa)VS 

ys < In 

y, > -,„(l^i±l^, (6.40) 
Pa 



in which Xa, Xb are the usual two-body fractions 

xa = P^(e^y^ +e^y^) . (6.41) 

In practice, we can ignore one or both of these limits if the relevant integrand is 
not singular as ys approaches infinity from that direction. Our initial-state soft pieces 
must therefore retain both limits, but our initial-state coUinear pieces can ignore one 
of the limits, and in the case of final-state pieces, both limits are ignored. 
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So, although naively expecting our soft limit to be adequately represented by the 
conditions "p3 — > 0, 1/3 and 03 fixed" ignores the fact that = p3 cosh 1/3 does not 
necessarily go to zero if is taken large enough, now that we have finite limits on 
ys, these conditions make sense, and we will write our soft limits in the form 

where 6a and Ob express the 2/3 limits above: 

Ob = e{ys + \n[{l-Xb)^/S/p^]) . (6.42) 

Unfortunately, the same notation won't work for the initial-state collinear terms. 
What we wish to express is that as p-^ becomes parallel with p^, for example, 
Ps-Pa ~ Pse"^^ andp3 alone go to zero, while ps-pb ~ Pse^^ does not. The descriptions 
we will use are: ^ 

^bbse^^ = 0,P3e-'^',PT3 = 0)eb , (6.43) 
the soft limits of which are 

^a(P3e^^ = 0,p3e-^' = 0,PT3 = 0)0a 

MPse""' = 0, P3e-y' = 0, PT3 = 0)9b , (6.44) 

and will be called soft- collinear limits. These will be subtracted from the collinear 
terms. 

■^Moving to light-cone coordinates pg = {p^ ,pT3}, in which p^ = (pg + p'^}/\/2 = p3e^^/%/2 
and pg = (pg — p3)/\/2 = p3e~^^/\/2 would allow us to write our collinear limits more cleanly, but 
at the expense of additional notation. We will have to be content with the above descriptions. 
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For m = a, then, we have in the end these pieces: 



4[123]x(e) 



OiP2 - P3)0a0b^a{P3 = 0,2/3,03) 



(6.45) 



PsPa ■ P3 




4[123]x(e) 



9{P2-P3)e, 



a 



PsPa ■ P3 



,P3e 



= 0,PT3 = 0) - J^aiPse^' = 0,p^e-y^ = 0,pt3 = 0) 



(6.46) 




X J^a{yi,P2,y2,(p2,P3,y3,h) - ^aiPae^^pse = 0,prp^ = 6)9, 

+ Ta{p3e'^'' = 0,P3e-^' = 0,PT3 = 0)^a - ^a(P3 = 0, 2/3, ^s)^,^^ 



a 



(6.47) 



Note that d a^^™*^ is finite everywhere, so we can in fact take e — > in this 
term, thus replacing (ie[123]x(e) with the notations (io[123]xo and D = 4 — 2e 
4. Meanwhile, d^a^"'''' has a collinear singularity but no soft ones, and d^a^"-^^ 
has both soft and collinear divergences. The sum da"^"-^^ + rfcr™" + rfa-^*™*^ ~ 
•^(2/1,^2,2/2,02,^3,2/3,03), which is what we started with. 

It will become evident in Section 6.2.5 that significant simplification of our 2/3 
integrals arises if we take advantage of the further breakdown daa = dcTab + dcFa2 
for the soft pieces {{mn} = {al,61,21} do not occur via the arguments of section 
6.2.1). In addition to rfcrf*™*^, then, we have just defined three subtraction pieces, 
two soft (labeled {mn} = {ab, a2}) and one collinear. Analogous terms are similarly 
constructed for m = b. 

The final-state pieces have their own organizational requirements, but thankfully 
are simpler as regards notation. For m = 2, we subtract the soft-collinear piece from 
the soft term instead, leaving it with no collinear poles. 
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jDsoft 

(1 C 2 = 

P3P2 ■ Pi 



d{P2 - Ps) 



X 



^2{P3 = 0, Vs, h) - ^2{P3 = 0,y3 = y2, h = <t>2) 



(6.48) 



.Daou _ 4[123]x(e) 
a an = 



P3P2 ■ P3 



0{P2- P3)^2{P3,y3 = y2,' 



(6.49) 



^4 mite 



cio[123]xo 

P3P2 ■ P3 



d{P2 - P3) 



^2{yi,P2,y2,4>2,P3,y3, 



-^2(^3, z/3 = y2, 



+ ^2{P3 = 0, ys = y2, 4>3 = <p2) - ^2{P3 = 0, ys, h) ■ (6.50) 

Again, it will be useful to break the m = 2 soft piece into two partial-fractioned 
terms {mn} = {2a, 26}. 

For m = 1, we will worry about a collinear piece only, as it can be shown that 
there are no soft poles: 



a cTi = 

Pi ■P3 



0{P2- P3)^l{P3,y3 = yi,' 



(6.51) 



4 finite 



rfo[123]xo 

Pi ■P3 



0iP2 - P3) 



X J^i{yi,p2,y2,(p2,P3,y3,h) - ^i{P3,y3 = yi,h = (l>i) ■ (6.52) 

Note that, as promised, we have switched back to Pi ^ P2 in the phase space of 
^D^fimte order to take advantage of one common event generator in our Monte 
Carlo program. 

In all, we have ten subtraction pieces; six correspond to soft singularities and are 
labeled {mn} = {ab,ba,a2,b2,2a,2b}. The other four correspond to the legs with 
which may be collinear, and these are labeled m = {a,b, 1,2}. The next section 
is devoted to performing the approximations herein defined. 
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6.2.4 Perform Approximations 

Since the Tm-, like tlie squared matrix elements from which they are derived, are 
given in terms of invariant dot-products "Pi -pj, we must first see what these look like 
in terms of our free variables {yi,P2, ^2, 02,P3, ^3, ^s}- 



2pa 


■Pb 


= X1X2S 


2pa 


■Pl 


= XiVSpie''^'^ 


2pa 


■P2 


= xi\fSp2e~^^ 


2pa 


■P3 


= Xiv^pse'^^ 


2pb 


■Pl 


= X2\fSpie^^ 


2pb 


■P2 


= X2\fSp2e^'^ 


2pb 


■P3 


= X2\fSp3e^^ 


2pi 


■P2 


= 2pip2[cosh(yi 


2pi 


■P3 


= 2piP3[cosh(yi 


2P2 


■P3 


= 2p2P3[cosh(y2 



in which 



y2) 


— cos(i 


h - h)] 






— cos(i 


h - 03)] 




ys) 


— cos(i 


h - 03)] , 


(6.53) 



2 



Pl = yV2+Vl + 2P2P3 COS (02 



tan- P2sin02+P3sin03 _ ^g ^^^ 

P2 COS 02 + P3 COS 03 



In the soft hmit, P3 — > while and 03 remain fixed. We then obtain pi p2 

and {xi,X2} — > {xa,Xb}, where 

= ^{e^y^ + e^y^) , (6.55) 

leading to the soft limits expressed in Table 6.1, in which the ab 12 subprocess 
invariants s, t, it are given (for massless particles) as 
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Table 6.1. Soft limits of invariants. 



2pa 


■Pb = s 








2pa 


■ Pi = -i 


2pb 


■pi = -u 




2pa 


■ P2 = -U 


2pb 


■ P2 = -i 


2pi ■P2 = s 


2pa 


■P3 = 


2pb 


•P3 = 


2pi ■P3 = 2p2-p3 = 



S = 2pa-Pb = 2pi ■ P2 

= XaXbS = 2pl{l + cosh(yi - 112)) (6.56) 

i = -2pa ■ Pi = -2pb ■ P2 

= -XaVSp2e~y' = -pI{1 + ey^~y') (6.57) 

U = -2pa ■ P2 = -2pb ■ Pi 

= -XaVSp2e~y' = -pI{1 + e^'-y^) . (6.58) 

In fact, for 2-body subprocesses, these invariants are not independent. For massless 
participants, t and u can be written in terms of s and a single angular variable v. 



t = -(i-v)s 

u = —vs 
1 

V = -(l + cosei) = . (6.59) 

In this way, we can show that, as expected, the soft limits J^miPs = 0, ys, (ps) reduce 
to products of 2-body matrix elements T'^j{v) and luminosity functions C{xa,xi,)- 
The denominators will retain certain y^, 03 dependences after approximation, and 
we will integrate over these in what follows. 

As an example of how the collinear limits are taken, we look at the case in 
which a gluon is radiated from one of the incoming legs of a subprocess, as in the 
diagram below. As the gluon (here pg) becomes collinear with the parent quark p^, its 
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momentum, as well as that of the daughter quark which enters the subprocess) 
are related to the parent in terms of a simple fraction z. That is, P3 = (1 — and 




Figure 6.2. Pd ■ Ps = collinear kinematics. 

The functions to be approximated, JF^, are in general functions of the invariants 
Pi -pj, where {i,j} = {a, b, 1, 2, 3} for a three-body scattering p^+Pt ~^ Pi +P2 +P3- 
^ In the collinear limit of massless particles, however, these invariants reduce to 
functions of z and the subprocess invariants (see Table 6.2). 







Table 6.2. 


Pa 


■ P3 = collinear limits of invariants. 


2pd 


■Pb 


s 






2pd 


■Pi 


= -f 2p, 


■Pl 


= —u 


2pd 


■P2 


= -f 2pfe 


■P2 


= -i 2pi-p2 = s 


2pd 


■P3 


= 2pb 


■P3 


= 2pi ■ P3 = -^i 2p2 ■ P3 = 



The fraction z is not arbitrary; its value in terms of our free variables is, as it 
should be, given by the {pse^^ , p^e'^^ = 0,^3 = 0} limit of the ratio Xa/xd- 

"^Labeling of the partons involved in a collinear splitting can be confusing. Having started with 
the 3-body labels ab 123, it would seem appropriate to leave the parent label unchanged (e.g. 
a), and add a tilde to the daughter (e.g. a). The convention, however, is that the parent always 
carries the tilde. 
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One can easily verify that in the soft-colhnear hmit, when pse"^^ itself goes to zero, 
z 1 and the invariants in Table 6.2 reduce to those in Table 6.1. 

For the m = 1 collinear contribution, becomes collinear with p^, and their sum 
balances P2- A new momentum fraction z is thus defined via pi = z{pi +P3), and 
the momenta of the incoming partons becomes the {2/3,03} {yi,0i} limit of xi, 
namely 

(see equations 6.31, 6.36). For m = 2, we have the {^3, (p^} {y2, ^2} limits 



P2 

z ~ 



xa = i^l±^(e±2^i +e^3'2) . (6.62) 



P2 +P3 
2+P 

In this way, all collinear Tm can be shown to factorize into two functions, a 
splitting function P[- which depends on only the relevant momentum fraction z, and 
a two-body matrix element T'^^{v) characteristic of the subprocess but independent 
of the earlier splitting. 

Adding up the contributions to each d^u^^^ and d^a'^'''' from all our 3-body 
matrix elements, and grouping by 2-body subprocess i (see Appendix A for a list), 
we can finally write out our subtraction pieces as follows. With the definitions (see 
equations 6.33, 6.38, 6.41, 6.42, 6.59): 



b 



P2 



+ e 



,yi 
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Xi 



Pi 



Pb 



4[123] 



p\^p\^ IP2P1, cos( 



x.^d 1 



•^a II 

—P-a 

Xi 

-Pl = x,^{l -1 
X2 " 2 

9{\n[{l-Xa)VS/p3-y3]) 
9{ys + \n[{l-x,)VS/ps]) 
dyidp2dy2d(pl~'^''dp3dy3d(j)l~'^^ 



aal^^" / P2 



27; \ 7^ I P2P3^Cf 



our soft pieces become, for subprocess i: 



(6.63) 



da': 



soft _ 4[123]x(6 



da] 



"'^"^ Pa-P3 

soft ^ 4[123]x(6 
4[123]x(e 



7 soft 
"'^a2{i) 



7 soft 
^^2a(i) 



J soft 
"'^2b(i) 



Pa-P3 

4[123]x(e 



Ph-Pz 
4[123]x(e 



P2-P3 

4[123]x(e 



P2-P3 



i^mx,,x,)K'T'^^iv) 

i)l^C{Xa:Xb)K'T'^j{v) 
^a2^Xa, Xb)K'T'^j{v] 
^b2^Xa,Xb)K'T'j,j{v) 
^m^a, Xb)K'Tjj:{v] 
'llj2b^Xa,Xb)K'T'j,j{v) 



Pa-Pb 



iPa-P3+Pb-P3 

Pa ■ Pb 
iPa-P3+Pb-P3 

Pa-P2 
iPa-P3+P2-P3 

Pb-P2 
iPb-P3+P2 ■P3 
P2 - Pa 



9adb 
Oa^b 
Oadb 
Gadb 



P2 



-Pa ■P3+P2-P3 P3 

P2 -Pb P2 

Pb-P3+P2-P3 Psi 



in which the [K'T'^j{v)] are given in Appendix A and 



(6.64) 
(6.65) 
(6.66) 
(6.67) 
(6.68) 
(6.69) 



= [€mi = (2Cf - Nc)S^,6^, + Neil - 6^,6^,) . (6.70) 

As an example, for the subprocess qq — > •yg, ipab — "^L — 2Cf — Nc, as both m = a 
and n = h refer to quarks, while '^\a2,2a,b2,2b} equal Nc for this subprocess. 



With the further definitions (see equations 6.36, 6.60, 6.61, 6.62) 
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Xr. 



Zb 



S 



Pi 



4[213] 



Pl +P3 

dy2dpidyi(f~^^(j)idp3dy3(f"'^^(f>3 



-2t 



f P3 



z = 



x^ 



{2^xyS\2^l^l. 

P2 

P2+P3 

Vs 



-2e 



PiP34Cf 



we find the collinear contributions: 



X 



Pa-P3 ' 
^l C{Xa/Za,Xb) -I , ^ 2C{a)Za 



da 



C{X(i^ Xb) 
coll 4[123]x(6 



b{i) 



b 



Pb-p-i 

C{Xa,Xb/Zb) 



1 - Za 

[Cixa,Xb)K'T'^jiv)ieb 



6„ 



dal 



coll 
(i) 



^{Xai Xb) 

4[213]x(e 



2C{h)zb 

I- Zb 



g3 



d(y^2(i) 



Pl -Ps 
4[123]x(g) 

P2 -PS 



[C{Xa,Xb)K'T'^j{v)]^Y.P' 



79 V 



[C{xa.Xb)K'T'Av)]XP'i. 



Here Sg^ = 1 if tlie radiated tliird particle is a gluon; otlierwise Sg^ = 0. 



(6.71) 



(6.72) 

(6.73) 
(6.74) 
(6.75) 



6.2.5 Perform ?/3-Integrals. 

We will now integrate analytically over the rapidity y^. This process will introduce 
further new notation, but as the functional form of these integrals is shared by more 
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than one subtraction term, the net increase in clarity should render the abbreviations 
forgivable. 

To start, the notation we're currently using for our free, unobservable phase-space 
differentials is (ie[123]. After integration over ^3 it may be easiest to simply write 

4[123 - ys] = dy,dp2dy2dcj)l-^'dpsd(j)l-^' . (6.76) 

Next, we'll see that for our initial-state subtraction terms, the limits on y^ (as shown 
in the functions 6a and 6b) will appear in the final result. These limits will be written 

as 

ylt^\n[{l-Xa,b)VS/p3] , (6.77) 

and will provide us, eventually, with the logs of Qx we expect. 

The integrals we'll need are performed in Appendix B; the results are as follows: 





jVZm 


ha 








112 {X) 




62 






= ±ln 



IL (x) 

2b 



-'2 



n = 



dyz 



p±y3 



g±{3/3-y2) 



- In x) ± - In 
2 ^ 2 



1 + x^^e-^y-i- 



1 + x^^e-'^y-i 



(6.78) 



I- cosh(i/3 - 1/2) - cos(</)3 - 02) 

±1 —OC.a 



^1 + 2xf\-^^y^3rn-y2) + 1 

^(l + 2x)^'e-2(?'L+3'2) + l- 
2{n - a{x)) cot a{x) - ln(l + 2x) + 2(?/|„ t Vi) 

dyz 



(6.79) 



^cosh(?/3 - ^2) - cos(03 - 

X 



xe^y^ + cosh(i/3 - 1/2) - cos(03 - ^2)-- 
ln(l + 2x) + 2(7r - (03 - ^2)) cot(03 - ^2) - 2(7r - a{x)) cot a( 



dy-i 



2(n 



00 cosh(y3 - ^2) - cos( 



dy?. 



PsZa 



X. 



H{Xa/Za,Xb} 
H{Xa,Xb) 



sin(03 - 02) 
2C(a 



x) 

(6.80) 
(6.81) 
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1 dz, 



H{xa/za,Xb) 5, , ^ 2C(a) ^ 



. "1/3- 



1 dz?, 



H{xa.Xb/zb)f,, , , 2C(6) 



H{Xa,Xt) 



1- Zb 



H{xa,Xb/zb) f^, 2C{b) 

L — Zb 



Zb [ H{Xa,Xb) 

in which we've abbreviated 
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cos a a; 



COS(^</)3 



VTT2x 
H{xi,X2) = 2X1X2S C{xi, X2) . 



The subtraction terms then become: 



7 soft 
^^ab{i) 


= de 


123 


- ysM^ 


^ soft 


= d^ 
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6.2.6 Map to Photon-Jet space 

It's now time to translate to a common set of observables {y^, pj, yj, (pj, Qt, 4>q} so 
that Born, virtual, counterterm, and subtraction terms may all be added and poles 
cancelled. As stated in Section 6.2, the poles in e will not cancel unless this is done 
first. 

For the initial-state contributions m = {a,b}, we simply have {yi,P2,y2,4>2} 
{y^, Pj, yj, (f)j} as in the two-body case, but as we have a non-zero p^ which contributes 
transverse momentum to the pair, we take ps — Qt and 03 — ^ 0g — vr. As this affects 
the photon transverse momentum, we also list the smeared dependent variables 



The final-state soft pieces {mn} = {2a, 2b} will not be resummed, so we 
could approximate the kinematics as well as the weights, and simply set Qt = 
here. Unfortunately, this would require an analytic integration over 03 (namely 
/ (i^~^^03(7r — a) cot a), and our lack of a closed form for such an expression makes it 
difficult to explicitly show the cancellation of poles. As these terms describe radiation 
which, in general, falls outside the jet cone, and thus contributes to the Qt of the 
photon-jet pair, it will be more natural to treat these terms in the same way as the 
{mn} = {a2, b2} pieces above. 

For the m = 2 collinear term, we unquestionably have zero Qt, and due to 
the collinear configuration of particles 2 and 3, they are counted as a single jet of 
transverse momentum p2 + P3- With the definition of z used earlier, we can take 



while {yi, y2, 02} {y-y, yj, 4>j} and p^ = pj, 0^ = (pj + n. As we now have purely 
2^2 kinematics for this piece, 0^ is undefined, and so we (as in the two-body case) 




(6.95) 




(6.96) 
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include a factor d}~'^'^(j)q/ J d}^'^^(j)q. The integrations over 03 and z must be done 
analytically. 

Finally, for the m = 1 collinear term, we do have a nonzero Qt = P3, and pj = p2 
as in the initial-state pieces. However, as we have switched Pi ^ P2 in the phase 
space for calculational simplicity, p2 is no longer a free variable. Instead we have its 
approximation pj = pi + ps. Thus we write pi{= p-y) Pj — Ps = pj — Qt- The 
splitting functions will continue to be expressed in terms of the momentum fraction 
z, where now z = {pj — QT)/pj- As for the rest of the variables, {1/1,112} — ^ {y-y^Vj} 
and 01 (= (p-y) (pj + TT, while 03 ^ 0^ — vr as in the initial-state terms. 

After performing the above operations, we can split each x and x factor into three 
pieces. The first, will be common to all contributions (2-body and subtraction), 
the second, (3^, common to all higher-order contributions, and the third, Xiacj will be 
expanded through the balance of the relevant contribution. With the definitions: 



d^l'jj] = dy.ydpjdyjd^ '^''cpjdQxd^' 

Oi(y.s Pj ( Pj \ An \ri^ .^M^'n^i 



2tx\M) ) r(l -2e) 



r2x ^ 



in which Mj is the factorization scale, we can now write all our subtraction pieces in 
the form: 



da(i) = 4[7i]«^iAr(i) , (6.98) 

where 
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psoft 
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(6.107) 
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6.2.7 Extract Poles in e 

Toward this end, we notice tliat for most contributions, the poles reside purely 
at Qt = 0, and can be extracted with the help of: 



Qt \Qt 
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dQi 



1 

Q^ 



2e 
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fiQr 
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I Qt 



2e2 



Qt 

KQt 



whereas for the m = 2 coUinear piece, there are poles at z = 1 and 03 
then need: 



(6.109) 

2e . 
(6.110) 

(hi + TT. We 



d 2 

111 



1 3 



27r2 



13 
Y 



■ 1 11 


/27r2 


67Y 


+ Nf 


rl 


23 1 


.~e ~ ¥ ~ 


^<3 


~ "9"/. 


.3^ 





Nc 



16^' 1^(1 -e) 



(6.111) 
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e r2(l-2e) 

Finally, the m = 1 coUinear piece has nonzero Qt (implicit in our definition of z) 
and therefore dependence upon a defined pair angle (pq. Unlike the m = 2 coUinear 
case, then, integration over the function /| is not independent of our free variables, 
and we must extract its (t>n 



d 



'j + vr pole in the sense of a distribution as follows: 
[9.(71 - Arh'] 1I 

(6.114) 
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with A0 = MOD(||0g - 4 
This brings us, finally, to: 
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2e . 
(6.117) 
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where the d in Z{d) is the type of parton (g or g) on leg (i (for subprocess ab —* ^d), 
and 
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Note that in T^abla a2 b2 ) there exist comphcated logarithms of the form 



In 



1 + xe-'^y^^ 



1 + ie-2?/L 



Since these tend to zero as Qt — ^ 0, we can drop them for the sake of simplicity. 
This amounts to a redefinition of our asymptotic form, and is valid as long as we get 
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rid of these terms in the finite piece also. This has been done in the Monte Carlo 
program. 



6.3 Add Born, Virtual, and Counterterm Contributions. 

We may now add in the two-body pieces from Section 6.1: 



de[lj]fii5iQT) 

de[l3]KAT^^) , 



(6.126) 
(6.127) 
(6.128) 
(6.129) 
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and the virtual parameters A^, Al, A'^, B"" are given in Appendix C. 



(6.131) 
(6.132) 

(6.133) 

(6.134) 



6.4 Cancel Poles and Take e ^ 

For simplicity of notation, we allow the virtual and counterterm pieces to share 
{mn} designations of their own, and write the sum 



Er 



mn{i) 



p I psoft I psoft I psoft I psoft 

J- Virt(i) + J- afe(i) + J- ba{i) + ^ a2{i) + ^ b2(i) 



TSOft 

2a{i) 



TSOft 

2b{i) 



114 



(6.135) 



which should now be finite. The full NLO result is then 

rfa = rf(T*^"3 + c/a^^"^*% (6.136) 
with given by the finite corrections of equations [6.47, 6.50, 6.52], and 



5(gT)+/9.Er 



mn{i) 



(6.137) 



After pole cancellation, of course, we can take e — > 0, and 
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With the exception of the m = 1 Bremsstrahlung pieces, which we will handle 
separately, each Tmn can be written in the form (subprocess label i suppressed) 
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(6.139) 



the coefficients of which are given in Appendix C. Note (there) that the sum of all 
2^mn coefficients is zero, as required, as is the sum of iTmn coefficients. Also note 
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that logs of Qr/Vj have been expanded into logs of Qt/Q plus logs of Q/pj in order 
to anticipate comparison with the expansion of the resummed form. 

6.5 Extract Resummation Parameters 

Each resummation formalism described in the text has, as parameters, a set 
of functions which depend on the particular process under study. In the 6-space 
formalism, these parameters are A, B, Ca/a, and C^^^, and we have seen in Chapter 
5 that these same parameters can be used for resummation in fc^-space. Now that 
we have a set of finite terms which together constitute the cross section to a fixed 
order (here NLO), we have the material necessary to calculate these parameters. 

This can be done in a couple of ways, each of which involves an expansion of 
the resummed form to the same order in as our perturbative result, followed by 
identification of A, B, Ca/a, and C^^i in the expansion with corresponding coefficients 
in that perturbative result. 

However straightforward this procedure may sound, it is not without its compli- 
cations. Although the 6-space and /cy-space resummed expressions are formally finite 
at Qt = (even without the non-perturbative parametrizations), this convergence is 
inextricably tied to the all-orders nature of the Sudakov exponentiation. Attempts 
to expand either result directly, to a fixed order, will lead to expressions which are 
unregulated at zero Qt- 

For the 6-space formalism, one way to resolve this is to "meet the expansion 
halfway" by comparing the Fourier transform of the NLO result to an expansion 
of just the 6-space integrand. Our starting point is equation 4.35, minus the finite 
correction: 



/TT J 

The C parameters live inside da^ , and can be expanded to NLO in as- Since 
we seek a form in which the partons are sampled at the same energy as in the NLO 




(6.140) 
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expression, we must also expand the parton distributions fa,b{bo/b) into fa^^Mf) plus 
appropriate convolutions over splitting functions P"*". The resulting expansion of da^ 
takes the form: 



da 
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a,b 



fa/A{Xa, Mf)fh/B{Xb, M j) 

1 dz 
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(6.141) 

The A and B coefficients live inside the Sudakov exponential S{b,Q), and can 
also be expanded to NLO in Performing the integral over /i in S, while keeping 
as{fi) — as{Mf) to first order, yields logarithms of Qb/bo: 



-Sib,Q) 



27r 



2A(^) In^ 9h + 25(1) 1^ qb 



(6.142) 



Multiplying these two expressions and keeping only terms of order 1 or gives an 
expansion which can be compared with the Fourier transform of the NLO result. It is 
perhaps easiest to go back a few steps in our NLO derivation to directly transform the 
regulated Qt^'^'^'' In™ Qt/Q pieces, but the transforms over plus-distributions can be 
done as well. Either procedure yields the required logs of Qb/bo, so that identification 
of the resummation parameters can proceed. The following integrals are useful in 
this regard: 
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(6.144) 



Since our plus-distributions have been defined with the upper limit pj, the Fourier 
transforms we would need are T^{bpj) and T}^{bpj). These give logs of bpj/bo, which 
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can be expanded into the required logs of bQ/bo, plus logs of Pj/Q- Since any finite 
upper limit to our Qr-integrals is immaterial, the latter logs should cancel similar 
logs in the NLO result (those created when the plus-distributions were introduced). 

That said, the actual method we will use to extract the resummation parameters 
is even simpler, and it relies on the fact that in the /cy-space formalism, the transform 
from 6-space back to transverse momentum space has already been done. We still 
cannot directly compare a finite-order expansion, but we can integrate both the 
resummed form and the NLO expression up to some arbitrary limit, and then expand 
the former. Beginning with equation 5.26 (without the remainder term TZ), we 
integrate on up to some p'^: 
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(6.145) 

Since the Sudakov exponent goes like S ~ X^m In™" Q/Qt, taking kills the 

second term above, and we're left with only the first term, which, when expanded. 
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Pt Pt\ 

To this we must compare the integral of the NLO expression. As can be seen in 
equation 6.139, we'll need only three integrals to perform this task. They are: 
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Note that the coefficients of the plus-distributions we've used are constant, and these 
distributions were defined with an upper hmit pj. The integral from to pj is 
therefore zero; the logs appearing in equations 6.147 arise solely from the integral 
from Pj to Pt- 

Applying these expressions to equation 6.139, followed by comparison with 
equation 6.146, yields the following resummation parameters: 
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Note that delegation of the Amn, -Bmn, and sCmn coefficients to one or the other 
leg (a or h) is arbitrary; we have chosen to assign half to each leg. We choose in this 
dissertation to resum only initial-state {IS) pieces, and for these, aDmn = bDmn = 1, 
so the In M|/p^ term drops in the above. Final-state {FS) pieces are left in the NLO 
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form shown in the last section. From here on out, we will make a distinction among 
the following four types of contributions: 

7 7 sinq I 7 finite 

da = dajg ^ + dajg 

+ daps' + daif. (6.149) 

dg-smg (regulated) asymptotic approximations, da-^*"**'^ are the finite correc- 

tions (compare equation 6.136). Resummation will be performed on only the first 
term, and afterwards we will call it daf^^"^^. The initial/final distinction cannot be 
applied to virtual terms; we make the choice to include them in dafg^, where they 
show up in the CSS C^^^ pieces. 

Note also that in this transverse momentum space comparison, we again obtain 
the logs of Q/pj mentioned above. As the sCmn pieces already contain logs of pj/Mf, 
addition of the new logarithms serves to produce pj-independent C^^-'-parameters, 
as expected. The new DmnP^^'^{M'j/pJ) terms also serve to cancel corresponding 
terms in aCmn and hCmn- 

Finally, a word on free variables: references to 6-space and /c^-space resummed 

expressions usually assume the set of free variables {<5, Qt, ^g, cos 6'css, 0css}) in 

which and (f)q are the mass and azimuthal angle of the observed particle pair, 

and {6*038! 0css} the direction of one of the particles in a specific 2-body rest frame. 

We, of course, are using a different set of variables, namely {y^,pj,yj, (f)j,QT, 4>q}, 

but given the existence of a nonsingular Jacobian between these two sets of variables 

(as there must be), there is sufficient justification to write the resummed result with 

either set of variables. ^ 

^Once the CSS coefficients are found, Monte-Carlo coding of the resummed piece proceeds 
as given by the CSS formahsm, in which the upper Hmit on Qt is the jj pair mass Q (which 

= y^2p|(l + cosh(j/^ ~ Vj)) ill om" variables), not the limit given by the theta function 9(j)j — Qt) 
at NLO. The question of whether this constitutes double-counting, or is physically justified given 
the multiple emission being described, will be discussed in a future work. 
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6.6 Matching Resummed and Fixed-Order Results 

Depending on the quality of the approximations used in defining the resummable 
asymptotic piece daf^^, there is no guarantee that it will continue to be a useful 
quantity at high pair Qt, either before or after resummation. In fact, it may even 
become negative: looking at the resummed fc^-space expression of equation 5.29, one 
can see that as Qt gets bigger, both the Sudakov and non-perturbative exponentials 
tend to one, while their derivatives tend to zero. This leaves only a term proportional 
to the derivative of the parton distributions with respect to scale. According to the 
Altarelli-Parisi equations, this contribution (and thus da^g^^"^) will be negative for 
all but the smallest momentum fractions (see figure 6.3). 

On the other hand, the 3-body perturbative piece (that is, the unregulated NLO 
result) is fine in the high-Qr region. The act of resumming its logs of Q/Qt should 
only make a difference in the Iow-Qt region. 

One is thus forced to make a switch from the resummed piece d(jf§^^^ to the 
perturbative 3-body piece at some value of Qt (see figure 6.4). The best method is 
still a matter of debate. Remember that in daj^g^^^ there live both the perturbative 
3-body result and a 3-body asymptotic approximation thereof. The resummed piece 
is a regulated and exponentiated version of the latter, and if the resummation 
procedure leaves its high-Qr behavior significantly unchanged from the asymptotic 
approximation on which it is based, there is hope that cancellation will occur between 
these pieces, leaving only the perturbative piece. If not, then this cancellation will 
require a little help. 

Kauffman [81] proposed inclusion of a function /{Qt) to smoothly force this 
cancellation as Qt becomes greater than some arbitrary matching value Q^^^'^^. That 
is. 



da ~ /{Qt, )(Resummed - Asymptotic) + Perturbative , (6.150) 
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Figure 6.3. Scale dependence of parton distributions, 
where, it was suggested, an appropriate function might be 
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(6.151) 



This allows all contributions to do their job at low Q^, but shuts off the resummed 
and asymptotic contributions at high-Qr- 
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Figure 6.4. Fixed-order perturbative result better at liigli Qt- 

Of course, it remains to determine Q™^*'^'^, and it is by no means certain that 
one particular value will suffice for all circumstances. If one's desired output is a 
single QT-distribution at fixed pair mass and rapidity, it may be easy to determine 
the optimum value by trial and error; better would be to histogram each piece 
separately and, after the full run, pick a Qr-value at which the resummed and 
perturbative pieces most closely agree, then drop all higher-Qr resummed and 
asymptotic contributions. 
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The procedure we will use in this dissertation allows the kinematics of the 
particular incoming state to determine the proper cutoff (as opposed to a fixed 
Qmatch-j^ However, it does so only by allowing one distribution to be binned at a 
time. 

That is, we assume an output consisting of a single one-dimensional histogram, 
the ordinate of which we will here label x. Given a domain and step size in x, we 
have then a defined number of bins in x. 

1. For each x-bin, we construct a two dimensional array in Q and Qt (up to 
\fS jl each). Thus in effect we have a 3-dimensional array in x, Q, and Qt- 

2. Proceeding with Monte-Carlo evaluation of rfo"'''^*"™', we fill this 3-D array with 
weights. 

3. We then collapse (integrate over x) this array into a 2-D array in Q and 
Qt only. For each Q-value, we start at Qt = and go up in Qt until we see the 
x-integrated weight dip below zero. This will define a cutoff Q^''^^'^^{Q) for each value 
of Q. That is, for each Q, all bins above this Qr-value are set to zero in the original 
3-D array, for all x. 

4. The original 3-D array may now be integrated over Q and Qt to determine 
the proper dafg^^"^ for each x. 

5. The array of cutoffs Qt^^'^^^Q) is now used in the subsequent evaluation of 
the 3-body weights (dcr-^*"**"^). As each 3-body point is generated, Q and Qt are 
calculated for this point, and if Qt > Qt^^'^^ for that Q, the asymptotic weight is 
not calculated, leaving only the perturbative. 

This procedure, while it has its limitations, nevertheless respects the approxima- 
tions in their regions of accuracy, while not forcing a unique matching Qt- 

However, for the sake of relative simplicity, we're not using the optimum cutoff 
criterion mentioned above, which would be to compare the resummed cross-section 
to the perturbative, and perform the cutoff at the Qt value at which they most 
closely agree. This would require additional arrays and introduction of some sort of 
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tolerance parameter. It could be done, but the results we've been able to achieve 
with the present system don't reveal any need for the additional overhead. 

In fact, at high enough hadronic-to-partonic energy ratios (low momentum 
fraction), the resummed piece goes smoothly to zero (without going negative) at high 
Qti and matching procedures of any sort are not required. The above procedure 
alters nothing in this case; by contrast, a fixed-Q™'^*'^^ procedure would kick in 
anyway. 
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CHAPTER 7 



RESULTS AND CONCLUSIONS 

7.1 Photon pt Distributions 

The direct photon spectrum is probably the easiest test of our formahsm, as 
there are data, and no correlations with jets to consider. We begin with results at 
collider energies, and end with lower energy fixed-target results. We use the fcy-space 
resummation formalism alone, at first, and use the nonperturbative parameters 
a = 0.3GeV~^ and kTHm = 4GeV for each subprocess. A discussion of these 
parameters, and why we may want to use different fits for different subprocesses, 
will follow in the next section. 

In figure 7.1, we show four curves: The LO contribution, the NLO results, and 
finally resummed results using the formalism, both with and without matching. 
Data are from the Fermilab DO experiment [82] . The process studied is pp — > 7+X at 
= ISOOGeV, with the photon rapidity constrained at |?/^| < 0.9 and factorization 
scale Mf = p^, where p^ is the transverse momentum of the photon. There is also 
an isolation cut on the photon, which rejects events with a jet of energy pj > 2GeV 
within R = \J {y^ — yj)^ + (0^ — 0^)^ < 0.4 of the photon. 

In figure 7.2, we show the same four curves, this time at the lower energy 
= 38.7GeV, and for the quantity E^d^a/d^p^, the invariant cross section, as 
opposed to d'^a/dpydy^f. Data are from the Fermilab E706 experiment [83]. The 
process is pBe — ^ 7 + X, with the photon rapidity constrained at —1.0 < y-f < 0.5. 

Notable here is the effect of our matching prescription (described in Section 
6.6), which cuts out the poor high-Qr behavior of the resummed and asymptotic 
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pieces. At ~ 3GeV, the incoming partons' momentum fractions are on the order 
X ~ 2p^/ \fS ~ 0.2, but as one approaches ~ 6.5 ~ 12.0GeV, these fractions reach 
X ~ 0.4 ~ 0.8, and in this range, the derivative of the parton distributions (which 
is dominant in the resummed piece at high-Qr) is negative. In contrast, matching 
at the higher energy ^fS = ISOOGeV is not an issue, as the momentum fractions 
concerned rarely go above x ~ 0.13. 

For completeness, we repeat the trials at the still lower energy = 31.5GeV 
(see figure 7.3). Data here are also from the Fermilab E706 experiment [83]. The 
process is pBe — > 7 + X, with the photon rapidity constrained at \y^\ < 0.75. For all 
these curves, the factorization scale Mf = is used, and the parton distributions 
are those of the CTEQ Collaboration, specifically CTEQ5M [84] . 
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7.1.1 Nonperturbative Parameter Choices 

As the lower-energy resuhs are much more sensitive to changes in the nonper- 
turbative parameters, we focus exclusively on comparison with the E706 results at 



In the plots discussed above, the /c^-space resummation formalism was utilized 
exclusively. In accordance with the work of Ellis and Veseli [77], this includes a 
nonperturbative factor of the form 



Besides a, there is also a second nonperturbative parameter, kTiim, which controls 
the boundary between perturbative and nonperturbative regions. We use a value 
kTiim = 4GeV, but changes of up to 50% in this parameter have little effect, as can 
be seen in figure 7.4. Thus, for the balance of this argument, we will think of the 
/cT-space method as having effectively a single-parameter nonperturbative form. 

Starting with a value of a = 0.3GeV~^, figure 7.5 shows the effect of raising and 
lowering a as kTiim is held fixed. As can be seen here and in the previous plots, our 
initial value works well with available data, but the results are significantly dependent 
upon the precise value at low v^. 



= 31.5GeV. 
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This being the case, it is important to compare these /c-r-space results with the 
predictions given by the 6-space formahsm, the parameters of which are, ostensibly, 
better known. Using the fe-space CSS formalism of Chapter 4, and the Sudakov 
parameters calculated in Chapter 6, figure 7.6 shows the Born, NLO, and resummed 
theoretical results compared with data. For the resummed curve, the 2-parameter 
nonperturbative model of Landry, Brock, Ladinsky, and Yuan (LBLY) [71] is used, 
with their best-fit parameters gi = 0.24GeV^ and g2 = 0.34GeV^. 

Figures 7.7 and 7.8 show, respectively, the result of altering gi and g2 by 50% 
each. As can be seen from the parametrization 



raising either gi or g2 dampens out more of the high-6 region, thereby broadening 
the Qt spectrum, and consequently raising the average recoil momentum. 




(7.2) 
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p Be — > y + X at root(S)=31 .5 GeV, -0.75 < < 0.75 
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Unfortunately, none of these curves agrees as well with the data as does the 
single-parameter, Q-independent /cy-space result. The first question to be answered, 
then, is whether the difference lies entirely in the Q-dependence of the 6-space 
parametrization. To do this, we first look at the approximate analytical relation 
between each parameter set and the average Qt "kick" given to the system. 

For example, in the 6-space model, to first order in we are resumming only 
the Born terms, and the 6-dependence of the Fourier integrand is simply that of the 
nonperturbative exponential, which we write as simply exp(— (76^). Then, 



da 2Q 



2tx J 



dQr Itx 

= 2Qt I dhhJo{hQT)e-3^^ 



9Ze~QV^9 ^ (7.3) 
9 



which has a maximum at < Qt >= \/2g. 
In the /cy-space model, we have 



da ^ ^ d 



dQr ^ ^'^^dQ^ 



1-e 



T 

= 2QTde-~^^'T , (7.4) 

which peaks at < Qt >= 1/V25,. Equating these two results gives a simple relation 
gd = 1/4 which we can use to test the effect of the Q-dependence in g. 

Figure 7.9 shows the result of a /c^-space run with a = 1/(4(7) superimposed upon 
a 6-space run of the canonical two-parameter form g{Q) = g2 In + gi. One can 
see that adding the Q-dependence to a reproduces the 6-space result. Conversely, as 
shown in figure 7.10, removing the Q-dependence by fixing g = l/(4a) also results 
in agreement between the b and /cj-space formalisms. 
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The next question is whether a Q-dependent form is truly required by the existing 
data, and for that to be answered we must look at the application of resummation 
methods to intermediate vector boson production, e.g. Drell-Yan processes, which 
comprise the bulk of the data upon which the current nonperturbative fits have been 
based. 

By simply swapping out our direct photon CSS parameters {Aj, Sj, Cj} for those 
calculated by Ellis, Ross, and Veseli [73] or Balazs and Yuan [85], we can begin 
comparing resummed theory to data for dimuon production, as was done in references 
[70-72] (see Section 4.2). Here we include only the resummed piece (which gives the 
bulk of the cross section in this region), not the finite corrections. Figure 7.11 shows 
an invariant cross section as a function of Qt for various Q-bins. The data are 
from the E605 experiment [86], and the theory curves are the 6-space resummed 
predictions using the two-parameter Q-dependent nonperturbative form. The values 
gi = 0.24GeV^ and g2 = 0.34GeV^ are the result of the global fit performed by 
Landry, Brock, Ladinsky and Yuan (LBLY) [71] , which considered data only out to 
Qt = 1.4GeV. 

However, as shown in figure 7.12, we find that using a single-parameter, Q- 
independent form, with gi ~ O.SGeV^ and g2 = O.OGeV^, produces a better fit 
to the shape of the cross section at every range of Q, especially in light of the 15% 
normalization uncertainty of the data. This is true also for the dimuon data of the 
R209 experiment [87], as shown in figure 7.13. These data have a normalization 
uncertainty of 10%, and were also used in the LBLY fit. 
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As seen in the previous section, kjspace resummation did well in describing 
direct photon data, as long as a was not too far from 0.3GeV~^. The value of d 
which corresponds to a 6-space parameter g = O.SGeV^, however, is a ~ 1/(4(7) — 
0.5GeV~^. Drell-Yan dimuon production, however, has only quark annihilation 
subprocesses as Born terms, while direct photon production has gluon-induced 
Born-level terms, so it is possible that a color-structure-dependent nonperturbative 
form is required. 

To check this, we repeat the direct photon predictions of the previous section, 
using single-parameter, Q-independent forms for both b and /c^-space resummations. 
For the qq contributions, we use g = O.SGeV^ and a = O.SGeV^^, and for the qg and 
gq contributions. 



g = (^£+^o.5GeV2 =0.83GeV2 
and 

a = l/{4g) = 0.3GeV^2 . (7.5) 

The results, for E706 data at ^/S = 31.5GeV and 38.7GeV, and for DO data 
at \^ = ISOOGeV, are shown in figures 7.14, 7.15, and 7.16, respectively. The 
factorization scale is Mf = for all plots. All show improved agreement with data 
for both b and /c^-space resummation procedures. 
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7.2 Other Parameter Choices 



7.2.1 Scale Dependence 

In a leading-order calculation, the Born terms are combined with energy- 
dependent couplings and distribution functions to produce a Leading-Log (LL) ap- 
proximation. As discussed in Chapter 2, when next-to- leading order contributions are 
added, the dependence on renormalization and factorization scales should diminish. 
We check this by looking at photon distributions in pBe collisions at a fixed target 
energy of 530GeV. We set the renormalization and factorization scales equal to each 
other, and run the LL and NLL calculations for 3 < p-y < 12GeV at three values 
of this scale. When we're done, for each bin, we plot the difference between the 
low and high scales in ratio with the value at the middle scale. The NLL calculation 
should show smaller percent differences, and does, as shown in figure 7.17. 

7.2.2 Jet Cone Width 

The jet cone definition we have used counts two partons as belonging to the same 
jet if they are within an angular radius R of each other, this radius defined by [80] 



However, our asymptotic pieces do not depend on this radius, as we have taken a 
fixed jet definition for each (see Chapter 6). The only possible dependence we have 
on this parameter comes via the three-body perturbative contribution, and then 
only through those observable quantities which are defined differently depending on 
whether or not the cone condition is met. For example, if we are looking at a photon 
p^ spectrum, with cuts on the photon rapidity |/^, we should see no dependence on 



the cone radius R, since P7 = ypi + + 2p2P3 cos(02 — ^3) and = yi everywhere 
in phase space. However, if we are looking at a Qt distribution, say, then raising R 




(7.6) 




149 




J \ I L 

4 fi 8 10 



Figure 7.17. Scale Dependence. 

should steepen the curve, since more weights are being binned as if the photon were 
recoiling against a single jet, i.e. at Qt = 0. 

Figures 7.18 and 7.19 show the cross section vs. Qt for those pieces associated 
with initial and final-state singularities. In each, the perturbative, asymptotic, 
2-body, finite (perturbative - asymptotic), and total are displayed. Here R = 0.2, 
6.5 < Q < 7.5GeV, -0.75 < < 0.75, -0.3 < yj < 0.3, and the reaction is 
pBe —>■ 'J + j + X . 
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p Be — > y + j + X at root(S)=31 .5 GeV 
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At present we know of no data for the Qt distribution of a photon/jet system, 
which as just as well, as there are caveats associated with predicting Qr-distributions 
in our current calculation. 

As can be seen in Chapter 6, in order to extract the logarithms of Q/Qt for 
resummation, we needed to analytically integrate over the rapidity of the third 
final-state particle. However, our jet definition depended upon ^3. We were forced 
to approximate the kinematics as well as the matrix elements, thus removing this 
?/3-dependence. For the initial-state pieces, we took Qt = ps everywhere, even though 
radiation inside the jet cone would have resulted in Qt = 0. This doesn't cause much 
of a problem for small cone width. For the final-state pieces, we treated the soft terms 
differently than the coUinear, and this can cause a mismatch even outside the cone. 
For example, if a 3-body event comes along for which pg is soft and just outside the 
cone, we have the following scenario: 



Table 7.1. Anomaly at low Qt- 



Piece 


...looks like: 


...but 


gets binned at: 


3-body Pert. 


SC 


Qt = 


Ps 


m = 2 (Soft - SC) 





Qt = 


P3 


m = 2 Coll. 


SC 


Qt = 






Here SC stands for "Soft-Collinear". The 3-body perturbative piece and the m = 2 
asymptotic coUinear piece both give roughly the same weight, but the asymptotic 
collinear is always binned at Qt = 0, while the perturbative piece isn't. If the bin 
width in Qt is larger than this ps, the cancellation will take place. If not, it won't. 

In figure 7.20, we show a comparison of Qt distributions for the m = 2 finite 
(perturbative - asymptotic) piece, one curve corresponding to the approximated jet 
definition we were forced to use, the other curve showing what would result were we 
able to use the non-approximate version. The latter tends nicely to zero as Qt 0, 
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The result is a non-continuous Qt distribution, and so further work is required in 
order to deal properly with these final-state terms. Eventually, of course, we'd like 
to resum these pieces, and to do so in a way that is independent of the particular jet 
definition in use. For the present, integration over Qt-, as in the preceding inclusive 
p^-distributions, gives a good, unambiguous result. 

7.3 Conclusions and Future Improvements 

We have seen the means by which soft radiative corrections can be resummed to 
all orders, and the significant effect such contributions have on the direct photon 
spectrum. In principle, any steeply-falling spectrum is a candidate for this procedure, 
including most single-particle inclusive transverse momentum distributions. As 
contributions from all orders are included, resummation improves agreement with 
data as compared to fixed-order calculations. Thus there is the potential for a 
significant impact on the prediction of a wide variety of processes. 

In this work, we have also seen that direct photon data can play a role in 
constraining the nonperturbative function; in particular we have seen evidence of a 
dependence on the color structure of the subprocess, in the form of mass-independent 
Gaussians for each structure. The nonperturbative parametrization is thus far from 
pinned-down, and it is likely that further exploration of perturbative methods can 
aid in probing the nonperturbative region. 

Proper matching of the resummed to fixed-order result has been shown to be 
important, especially at fixed-target energies. The matching procedure outlined 
in this work accomplishes this task in a fiexible and automatic manner, without 
infiuencing the result in any way when not required. 

The FORTRAN program written concurrently with this dissertation, and from 
which the results of this chapter were produced, is available upon request. It is 
capable of producing either NLO or resummed output, as the user desires. 
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Without significantly altering the adopted procedure, it should be possible to 
avoid problems associated with jet definitions by looking at photon plus pion final 
states, for which there are a good amount of data [83]. Furthermore, Qt and pout 
distributions exist for these data, providing a more sensitive base upon which to 
examine nonperturbative parametrizations. 

Inclusion of threshhold effects is also possible, as indicated by the recent work 
of Laenen, Sterman, and Vogelsang [88]. This method involves both Fourier and 
Mellin transforms, and looks, though complicated, to be quite powerful. 



156 



APPENDIX A 



DIAGRAMS AND MATRIX ELEMENTS 



FOR NLO PHOTON 



HADROPRODUCTION 



In this appendix, we include Feynman diagrams of the contributions to photon 
production by hadrons, first to leading order, then to next-to-leading order. 

Figure A.l shows the basic couplings involved. The first is a QED-like photon- 
quark coupling of strength cqed- The second is the analogous QCD gluon-quark 
coupling, this time of strength gg. The third, a gluon-gluon coupling, has no QED 
analogue. In general, all non-Abelian theories (such as QCD) will have self-couplings 
among the gauge fields. Note that in QCD there is another gluon coupling (with four 
attached propagators), but as it is of order g"^, it will not be needed here. 



To find the leading-order contributions to photon hadroproduction, we need to 
put these vertices together (along with suitable propagators and external legs) in all 
possible combinations consistent with the following conditions: 





Figure A.l. q^y, qg, and gg vertices. 
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1. There must be one final-state photon. 

2. There must be two incoming particles, each either a quark or gluon. 

3. Each diagram must be simply- connected and maintain momentum conser- 
vation. 

4. Each diagram must use as few vertices as possible. 

5. Each diagram must be topologically distinct. 



These conditions constrain us to the diagrams shown in figure A. 2. 



I . I D .1 



qg-* jq 





g q 1! q 



Figure A. 2. Leading order photon hadroproduction diagrams. 



Now we go to next-to-leading-order (NLO). We first seek the basic simply- 
connected Green's functions (that is, tree diagrams) which include one photon-quark 
vertex and two vertices from the set {qg,gg}- Given the "connected" constraint, we 
rule out two gg vertices, and arrive at the following four diagrams (fig. A. 3). 
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Figure A. 3. qq'ygg (a, b and d), and qq'yqq time-independent Feynman diagrams. 

Note tliat we have not yet sliown any time direction. That is, we have not assigned 
any of the external lines to "incoming" or "outgoing" legs. Figures A. 4 and A. 5 show 
the range of such associations, given that the incoming particles (on the left) must 
be either quarks or gluons (this is hadroproduction, after all), and the photon must 
be outgoing. In addition, we will have next-to-leading order virtual contributions, as 
shown in figure A. 6. 

Finally, there will be Bremsstrahlung contributions, which arise when one final- 
state parton in a purely QCD subprocess fragments and a photon is produced. The 
signature for such an event is a non-isolated photon (one with hadronic energy 
deposited in a small cone about the photon axis) and a jet in the opposite direction. 
The purely QCD 2-body subprocesses are as follows; diagrams for these, as well as 
the derivation of their matrix elements, can be found in Ellis and Sexton (1986) [89]. 
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A.l Two-Body Matrix Elements 

For dimension D = 4 — 2e, and with the definitions 



K' 



s 
i 
u 

qq 



K' = K' 

qg gq 



^2Q 



{Pl + P2f 

iPl-P3? 

{Pl - Pif 
1 

3 

1 

2(47r)'CFa,V 



we can write the two-body matrix elements as [89-91]: 



A. 1.1 q{pi)q{p2) l{P3)g{Pi) 



= K',,Qlu;,,{l - e)[{l " + j 



2e 



A. 1.2 q{pi)g{p2) l{p?>)q{Pi) 



t s 
s t 



2e 



A. 1.3 g{pi)q{p2) l{Pci)q{Pi) 



nM\' = -K',,gW(l _ e)f(l - e)(5 + ; 



u s 
s it 



-2e 



A.1.4 q{pi)q{p2) ^ g{Pz)9{Pi) 



S|M| 



K' 



qq 



gg^gg 



^2q(i - e; 



Cfj:: - Nc 
tu 
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A.1.5 q{j>i)gij>2) q{p3)g{p4) 



m\"=-K',,uj,Qii-e: 



— e 



Cf^ - Nc 
su 



A. 1.6 g{pi)q{p2) q{Pz)9{PA) 



£2 + 



— e 



st 



A. 1.7 g{pi)g{p2) q{Pci)q{PA) 



E|M| 



^2Q — e 



ut 



A. 1.8 g{pi)g{p2) g{P3)g{pA) 



E|M| 



99 



99—99 2 



QNcUJ2Q 



tu SU st 
P P V? 



A. 1.9 q{pi)q{p2) ^ qipMPi) 



S|M| 



99—99 4 



^2Q 



t2 n2 



A. 1.10 q{pi)q{p2) q{Pz)q{Pi) 



E|M| 



99—99 2 



^2Q 



2 I "2 -2 I X2 



2e- 



2(l-e) /n 



A. 1.11 q{pi)q{p2) q{p3)qiP4) 



S|M| 



rS2 + P £2 ^^2 



99 



99—99 2 



'^2Q 



+ 
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A.1.12 q{pi)q'{p2) ^ q{P3)q'{Pi) 



2 

qq'~*qq' 



K' 



11 , 



(1 



US 



A. 1.13 q{pi)q'{p2) ^ q'{Pz)q{PA) 



K' 



m\qq'^q>q = ^-^2Q[{l-e) 



A. 1.14 q{pi)q{p2) q'{P3)q'{PA) 



S|M| 



11 



qq~*q'q' 



^2Q 



(l-)-2t 



Note that in the text, we will often use the definition 



^2Q 



(A.14) 



(A15) 



(A.16) 



(A 17) 



in which v = —u/s. Use of either K oi Tjj without the prime refers to the e 
limit thereof. Here and in the text, Cp = 4/3 and Nq = 3. 



A. 2 Three-Body Matrix Elements 

With the definitions 



Sl2 






tl3 






tl4 






^15 








= ip'2 




^24 


^ ip'2 




hb 






S34 




+p'.f 
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S35 

K' = K' 

UJ3 



= ipl+p',)" 
1 

3 



1(1 -e) 



2 ,,6e 



4(47r)^Ciraa>' 



(A.18) 



we can write the 3-body matrix elements as follows [91]: 

A.2.1 q{pi)q{p2) l{P3)9{p4)g{p5) 
We need: 



M, 



QI199 



X 



,■512, ^135 ^14, ^15; ^23, ^24, ^25; ■534, S35, S45) 



+ 



+ 



6(1-26) 



S12 , Nc tu , Nc t24 



2 ^15^45 
+2 I +2 



+ 



C 

2 / ^15^25 

^13 + ^23 _^ ^14 + ^24 _^ ^15^25(^15 + ^25) 



2 t25^S45 
(-2 



^14^24 ^13^23 
■(tl4 — ^24) (^15 



^13^23^14^24 



6(2 + 6 
26(4 - 6 



r^l3^23 



^14^24 



^14^24 



^25) _|_ (^13 "~ ^23) (^15 ~ ^25) 



+ 



^13^23 
^15^25(^13 ~ i23)(^14 



t24) 



+ 



''15''25 



^13^23^14^24 



^13^23 ^13^23^14^24 
^ _l_ ^15^25 _|_ ^15^25 



tl3t 



13''23 



14 ''24 



(A. 19) 



from which we obtain 



- 2 " 

S|M|g-^^^^ = 2 Qq^3Mgggg{Si2, ^13, ^M, ^15) ^23) ^24, ^25) -534, -535) ■S45) • (A. 20) 



A.2.2 q{pi)g{p2) l{p3)q{P4)g{p5) 



E|M| 



-K' qgQgU^Mqgggi^tu, 1 , S 12 , ^15, ■S34, ^24, ■S45; ^23; ■535; ^25) • (A. 21) 
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A.2.3 g{pi)q{p2) lij>z)qij)i)gij>5) 



^l^lgq^^gg — ~K' gqQgUJ3Mgggg{t24:, ^3 , ^ 12 , t25, S34, ti4, S45, ti3, S35, tis) . (A. 22) 

A.2.4 g{pi)g{p2) l{p-i)q{Pi)q{Pb) 



- 2 

S|M|^^^^g- = 777^—7:; 7;-<35'^3M}qgg(s45! ^35, ^25) ^15) -5345 ^24, ^14, ^23, ^13) S12) ■ 



A. 2. 5 q{pi)q{p2) l{P3)q{Pi)q{Pb) 
With 



,^ ~Si2 ^25 . ^24 1 ^14 S45 ^ ti5 

ikZo — "; — ; r h h 



(A.23) 



^13^23 ^ ^23"S35 ^^23534 ^13"S34 535534 tl3S35 
j^^^ _ '^12 + ^25 + ^14 + "^45 



_ ^12 + ^24 + ^15 + '^45 



^25^14 

M3 = (.12 + .45)%^ 

J24tl5 

M4 = (.12 + .45)%^ 

^25^14 

-^5 = (.12 + ■^45) (2.12.45 ~ 2^25^14 ~ 2^24^15 ~ (^25 + ^14) (^24 + ^15)) 
+ (.12 — .45)^(2(^25 + ^14)(^24 + ^15) + ^25^14 + ^24^15) 
+ ((^25 + ^24)^ — (^14 + ^15)^) (^15^^14 — ^25^24) 

— ^15^24(^24 ~ ^15)^ ~ ^14^25 (^14 ~ ^25)^ ~ ^15^25(^15 + ^25)^ ~ ^24^14 (^14 + ^24)^ 
+ 8^25^24^14^15 + (^24 — ^25) (^15 ~ ^14)(^24^25 + ^M^is) (A. 24) 

and the further definition 

Mqqqq (.12, ^13, ^14, ^15, ^23, ^24, ^25) .34, .35, .45) 

= |2Mo(Mi + M2 -2e) 
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+ 



X 



X 



+ 



4e 



Vtl3'534 23 -535 
26 ' ■ 



V^13^23 ^23-535 ^13^34 "S35S34/ 



|(Si2 — 545)^ + Mti4 — t 



25 J 



-S12 



^24^15 
^24 S45 



( 



'ti3t23 t23Su S35S34 
2(512 — S45) + 2(^15 ~ t 



+ 



- + M1 + M3 

^15 N 



2) 



^13'535 ' 
n2 



24 J 



ty.^t 



+ M2 + M4 - 2) 



25''14 



.Mo 



{Si2 + ■545)(Si2S45 — ^14^25 ~ i^^l5^24) 
^25^24^14^15 



+ 4e 



2e [S35(si2S34 — ^23^14 " ^24^13) _|_ ^34(512535 — ^25^13 ~ ^23^15) 



c 



^25^23^13^15 



^23^24^14^13 



^23(~-S45^13 + ^14^35 + ^34^15) ^13(~^23'S45 + ^25'S34 + ^24^35) 
^25^24534535 

e M0M5 \ 

Nc ^25^24^15^14 I 



^14534535^15 



(A.25) 



we can write 



E|M| 



K' 



-CF<5gCi^3^gggg(5l2, ^13; ^14, ^15; ^23; ^24, ^25; 534, 535, 545) . (A. 26) 



A. 2.6 q{pi)q{p2) ^ l{P3)q{Pi)q{Pb) 



S|M| 



99 



CpQ'^UJ^Mggqg^tir,, ^13, ^14, 5i2, 535, 545, ^25; -334; ^23; ^24) • (A. 27) 



A. 2. 7 q{pi)q'{p2) ^ l{P3)q{PA)q' {p^) 
We begin with the definitions 



Mo 
Ml 



-S12 



24 



^13^23 ^235^34 535534 
2(512 — 545) + 2(^15 ~ ^: 



545 _|_ ^15 



^13535 
\2 



24; 



^25^14 
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^25^14 



tl3'534 123*35 
M, = (^l^ + ^^5)(t24 + ti5) ^^_2J 
^25^14 



Mgpgp (Sl2, tl3, ^14, ^15; ^23; ^24, ^25; S34, S35, S45) 

= 2(M2 - 2e){QgQg'Mo + M3) - 4e(M2 + M4 - 2)Ms 

- 2e(Mi + M2 + M4 - 2)QgQg,Mo . (A.29) 



and arrive at 



2 

S|M|gg,^^gg, = — -— Ci?U;3Mgpgp(si2, tl3, tl4, ^15; ^23; ^24, ^25; ■534, S35, S45) . (A. 30) 



A. 2.8 q{pi)q'{p2) ^ l{vMPi)(l' iVb) 



- 2 

S|M|^-,^^g-, = CFUJ3Mqpgp(ti5, ti3, tu, S12, S35, S45, t25, S34, ^24) • (A. 31) 



A. 2.9 q{pi)q{p2) ^ l{Pz)q' {PA)q' {Pb) 



S|-M|^-^^^,-, — —^CpUJsMqpqpitii, ti3, S12, ^15, -534, ^24, S45, ^23, ■S35, ^25) • (A. 32) 
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Figure A. 4. Time-ordered Feffi&ian diagrams from (a) and (b) 



Dia^riiTis fror"i (c) 



11 ::■ y 1 1I 



'--lLv.--.v-.v.i..: 
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THe longEst veiTicsl hnr indicates the pnrticlE referred tc 'Extra" in the text. 

168 

Figure A. 5. Time-ordered Feynman diagrams from (c) and (d) 




+ ghost loops 



Figure A. 6. Virtual Contributions. 
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APPENDIX B 



USEFUL INTEGRALS 

Here we collect the details of the and (f)^ integrals used in Chapter 6. Certain 
sub-integrals have been found in the tables of Gradshteyn and Ryzhik (G.R.) [92]. 



B.l /, 



ab 



dy3 



xe 



-2?/3 



2/3 



ln(l 



xe 



-22/3 ^ 



Vsm + + o In ■ 



xe 



1 + xe^^sr, 



(G.R. 2.313.1) 



(B.i; 



Similarly, 



V3n 



1 



1 + xe'^y^ 



1 , 1 + xe-^y's'^ 

- m — 

2 1 + xe^y^m 



(B.2) 



B.2 Ia2, h2 



ysr, 



dy-i 



xe'^^y + cosh Ay — cos Acj) 
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^3m 



dAy 



with Ay = 1/3 - 1/2, Acj) 



y-im-y^ xe~^y + cosh Ay — cos A(j) ' 
= 03 — 02- With the change of variable t] 



(B.3) 



(^71,0 = e j^^^ _ gi/3^ (.g^j^ write this as 



VHI 



2r]dr] 



VLO 



rf — 2?7 cos 



2x) 
(i?7 



2r]dri 



VLO r/2 + (l + 2x) 



2ry 



r/^ - 2?7cos A0 + (1 + 2a;) r]"^ + {I 



2x) 
(B.4) 



where the new limits are justified because in either the soft or coUinear regions, 
{^Sm' Vsm^ v} go to infinity, and the second integrand above goes to zero. 
The first term is easy: 



drf 



r/2 + (1 + 2x) 



In 



e^ivsnr-y^) + (1 + 2x) 
e-^iyL+y^) + (1 + 2x) 



(B.5) 



while the second term, after defining cos a = cosA0/vl + 2x, can be partial- 
fractioned as 



drj 



2r] 



2i sin Q!a/1 + 2x 



1 



7] — 



(B.6) 



-1] — + 2x 7] + i^Jl + 2x^ 

The first set of fractions has the form z — z* = 2iIM.{z) for complex z, while the 
second set looks like z + z* = 2TZ8{z). In this way, we can rewrite the integral as 



drj- 



/o \Jl + 2x sin a 
where 



XM 



V 



r] - e*"vr+2x 



dri2nS 



7] — iy/l + 2x 



(B.7) 
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IM 



IM 



1+ 



?7 - e*"Vl + 2x 



7] - e*"Vl + 2x, 



(B. 



and so 



(i?7 



VI + 2xsi 



sm a 



7] - e*°Vl + 2x 



2IM 
2IM 



oo 



(i?7 



sm a Jo 7] — e*"Vl + 2x 

ia 



sm a 



In {ri - e^°Vl + 2x) 



Meanwhile, 



dr]2nS 



Tj — ii/l + 2a;. 



27^£ 



In (?7 — iVl + 2a;) 



So, since 



cos a + z sm a 



In (r/ - e*"Vl + 2x) 



In y (r/ - e*"vr+2a;)(r/ - e-*"Vl + 2x) 
1 



+ ztan 



— "v/l + 2a; sin « 



?7 — vT+2a; cos a. 



we have 



oo 




(B.9) 



(B.IO) 



fB.lll 



XM 



e^" ln(r/ - e^"Vl + 2x) 



= sin a In Y + (1 + 2x) — 2?7VT+~2a' cos a 
+ cosatan"^ 



-Vl + 2xsina 



r/ — Vl + 2x cos a 



2ns ln(r/ - i VI + 2x) = 2 In V (?] - iVl + 2a;)(r7 + iVl + 2x) 

= ln(r?2 + (1 + 2a;)) . (B.12) 



Putting it all together, we obtain 
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In 



77^ + (1 + 2x) - 2?7V1 + 2x cos a 



7^2 + (1 + 2x) 



cot atari 



2 cot a [vr — a] 



+ 2x sin a 
T] — a/I + 2x cos a. 



and so 





B.13) 



Ia2{x) = In 



r e 



2(3/3^-3/2) 



:i + 2x) 



-e2(-yL-s/2) + (1 + 2x) 
Ib2ix) will be similar. We here have 



2 cot a[7r — a] 



(B.14) 



Ih2(x) 



dAy- 



-Ay 



lyi+yl^ ' xe^2/ + cosh Ay - COS A0 ' (B.15) 
which, by comparison with equation B.3, can be obtained via the negative of Ia2 with 
the signs of the limits reversed. We then get 



Ib2{x) = In 



eMm+y2) + (1 + 2x) 



-e2(-3/3™+3/2) + (1 + 2x) 
Again, cos a = cos A0/V1 + 2x. 



+ 2 cot aliT — a] 



(B.16) 



B.3 / 



2a,2& 



dAy- 



X 



- 1 



cosh Ay — cos A0 ^xe + cosh Ay — cos A0 
To start, consider the first term, which can be partial fractioned as follows: 



(B.17) 



haix) 



dAy 



00 
00 



2 / dr] 
Jo 



cosh Ay — cos A0 xe + cosh Ay — cos A0 
?7 ?7 



77^ — 2?7 cos A0 +1 rf — 2ri cos A0 + (1 + 2a:) 



[B.18) 



in which we've defined 77 = e^^. Each of these terms can itself be partial-fractioned 
and simplified as follows: 
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V 



V 



— 2?7 cos A0 + 1 (?7 — e*'^'^) [f] — e 



2i sin A0 
sm A0 



?7 — e*^*^ rj — e' 



-iA<f) 



■T] — e 



while, with the designation cos a = cos A0/a/1 + 2x, we also have 



(B.19) 



772 - 2?7cosA0+ (1 + 2x) {r] - y/TT2xe'°'){r] - y/TT2xe' 



V 



V 



2i sin ay/1 + 2x^7] - y/1 + 2a;e*° 77 - a/1 + 2x6" 



7/ 



sin aVl + 2x '-77 — Vl + 2xe* 
We're now left with the simple integrals 



(B.20) 



dt]- 



V 



drj- 



[t]- e^^<^) 



(r^-e^^-^) + e'^^ln(r/-e 



(?7 — Vl + 2a;e' 



VI + 2a;e'^'^ln(r/ - Vl + 2xe 



ih.4>\ 



(B.21) 



the imaginary parts of which we now evaluate (before taking limits). Of use will be 
the identities 



COS ± 2 sin - 
1 



In 2; = - In zz* + ztan ^ 



IMiz) 



TZSiz) 



(B.22) 



which lead to 



17 A 



— sm 



+ cos A0tan 



- sin A01n(?7^ 
— sin A0 



2?7cosA0 + 1) 



7] — COS 



(B.23) 



IM 



a/1 + 2xsinQ; 



ln(r7^ - 2?7a/1 + 2x cos a + (1 + 2x)) 



— Vl + 2x sin 



a. 



Vl + 2a; cos atan 



sma 



•x] — \J\ + 1x cos a 



[B.24) 



Adding these two results, with the proper coefficients from equation B.19, and taking 
the indicated hmits results in: 



haix) = ln(l + 2x) + 2 cot A0(7r - A(/)) - 2 cot a{7r - a) 
From equation B.17, we see that we still need the integral 



(B.25) 



hnix) 



dAy 



oo 
oo 



cosh Ay — cos A0 



^2 



T]^ — 27] COS A0 + 1 

°° dr] r 1 



2zsinA0 
^ IM 



p—iA(j> 



sin A0 
2 

sin A0 
2 

sin A0 
2 

sin A(j) ' 



dr] 



T] — e 



r] - e*^<^ 



XA^ ln(?7 — e 



oo 




-tan 



TT 



sin A0 



7] — cos 



A4>) 



(B.26) 



This then makes 
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= ln(l + 2x) + 2 ^"""^ 7 ^ V - - 2 cot a(7r - a) . (B.27) 
sm 

Since J26 is the same as /2a with Ay ^ —Ay, but both are symmetric with respect 
to this transformation, we have l2b{x) = l2a{x). 



B.4 /| 



d 



sin A0 



r2(7r- A0)(1 - cosA(/)) 2(7r - A</)) cos A(/) 



sin A0 



sm 



fB.28) 



where the last term is 

j ^i-2.A^^ 2(7r- A0) cosA0 
7 sin A0 

Now, since cos A0 is odd around -7r/2, we have 



2 / dVLi_2, j d^<j) sin"^-^' A0 cos A0(7r - A0) 



(B.29) 



Tx \ dA(j) sin"^^^' A0 cos A0 = , 
Jo 



while 



/■TT 

/ c?A0sin"^"^' A0COS A0A0 
Jo 

can be integrated by parts to give 



A0 _,_2. 



-2e 



sin"^' A0 



dAd) 

^ -2e 



2e 



2er 1 - e 



Since 



(B.30) 



(B.31) 
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d^l2-2t 



it follows that 



T(l-e) 



(B.32) 



2(7r - Acf)) cos Acf) 



sin A(j) 

and thus, in the sense of a distribution, 



d 



l-2e. 



fB.33) 



3-'2isy3; 



d 



l-2e. 



2(7r- A0)(l-cosA0) 1 



(B.34) 

sinA(/) 

It will be useful to perform this integral fully as well. Concentrating on the first 
term, we can choose a change of variable w = {1 + cos A0)/2, which leads to 



and the integral 



sin A(j) = Aw{l — w) 

dw = —-sinA(f)dA(j) 



A(f) = -2i\n{V^ + iVl - w) , (B.35) 



2(7r- A0)(l-cosA0) /"i , 

d^ ^^(j)3— / \ — = / dwA ^ (l-w) 



sin A0 Jo 

X f7r + 2iln(v^ + iVl-M;)) . (B.36) 



The log term can be integrated by parts. With 



u = ln(v^ + iVl — w) 
dv = dww~^~''{l — w)~'' 



[B.37) 



we find both du 



idw 



and 
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V = dww "(1 — w) 



w ^(1 — w)' 



dww "(1 — w) 



(B.38) 



the latter term of which is / dww ^ "^{l — w) ^ + ^w "^(l — w) by symmetry. We 
then have 



w "(1 — w)' 



so that 



/ dww ^ "{1 — w) ^ \n{^/w + — w) 
Jo 



w ^(1 — w)' 



w 



This gives, finally, 



d 



l-2e. 



2(7r- A0)(1 -cosA0) 
sin A0 



from which 



(B.39) 



— / dww 2 '^(1 — w) 
2e Jo 

_^7r r(l-2e) 

2e r2(l -e) ■ 



(B.40) 



16"T^(1 - e) 



r2(l-2e) Jr(i-e 



2tt 



l-e 



16~T4(1 - e) 



r2(l-2e) 



d'-''<Psm,) = 
follows by addition of equation B.33. 



16"^ T\l - e) 
e r2(l-2e) 



jl-2e. 



(B.42) 



B.5 I2a{x)-K 
This integral is not used in the text, but may be useful to those interested in 
treating all final-state pieces with a common Qt = definition. The crucial integral 
here is 
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J (f ^^A</)(7r — a) cot a = J d^li^2t J (iA</)(7r — a) cot a sin ^^A</) 

= h{h-2eh) (B.43) 



to 0(e), where 



2 rfl-e) 4-^r2fi 



I ^ (47r)'r(l -2e) tt r(l - 2e) ' 



r7T 

I2 = dA(f)(7r — a) cot a 



= —h{x) , where 
2 /""" 

h{x) = - (iA0(7r - a) cotalnsin A0 . (B.44) 

TT Jo 

lo can be evaluated as follows: 



sin A(f)dA(f) 

d cos a = — sin ada = 



sinA0 = y 1 - cos^ A0 = Vl + 2xy sin^ a 
a(A0 = O) = cos~^^ ^ 



2x 



l + 2x 



.VT+2x. 

a(A0 = 7r) = cos"^( ^^_L- ) , (B.45) 



and so 



cos-i(-l/VT+2^) (ia{7T - a) 



COS a 



i(i/vT+2^) Jsin^a- 



(B.46) 

l+2x 

Now, with TT — a = | + (^ — a), this can be broken up into two terms, the first of 
which is odd about 7r/2 and the second even. Since our limits of integration straddle 
7t/2 evenly, integration over the first term gives zero, while the second term gives the 
same integral from '^os~^{^^/=^) < a < | as it does from | < a < cos^"'^(;^7=j). 
Thus 
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da4^±2^ . (B.47) 
cos-i{i/VT+2^) y'sin^ a - 

The first term here can be evaluated using the change of variable y = sin a, which 

yields 



dy 

Vo = Z/(«o) 



cos ada 



sm cos 



and thus 



2x 



vrT2^ Vl + 2x 



(B.48) 



TT 



dy 



nln2(^y+ ^Jy^ - y^J 



yo 



Trln (B.49) 



2x 



by G.R. 2.261. The second term can be obtained via G.R. 3.842.2 with u 



cos^i(l/Vl + 2x): 



'^/^ daa cos a 



sm^ a 



2x 
l+2x 



IT ln(l + cos u) = TT In 



1 + VTT2X 
VI + 2x 



Thus 



, 1 + Vl + 2a; , 1 + VI + 2x tt , 1 + 2x 

io = TT In ^ TT m , — = — m 

V2x VI + 2x 2 2x 



and so, to C(e) 



(B.50) 



(B.51) 



d^~'^'A(f){n - a) cot a = 4" 



r2(l-e)rl l + 2x 



In- 

2 2x 



r(l-2e) 

Now, from equations 6.80 and 6.81, we find that 



eh{x] 



jl-2e. 



(B.52) 
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d^-^'A(j) ln(l + 2x) -2 



(tt - A0)(1 - cos A0) 



In (1 + 2x) 



1 16^^ r^(i 



sin A0 



2(7r — a{x)) cot a;(x) 



e r2(l-2e) 



-4- 



X 



In 2x - In 16) + e[2h{x) + In 4 In 



r(l-2e) 
2x 



In 



2x 



2x 



In^ 16 



2x 



IT 

~3 



2eh{x] 
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APPENDIX C 



DETAILED NLO RESULTS 

If we write out our Tmn factors from Chapter 6 grouped by QT-dependence (and 
including poles): 



= S{Qt 
+ S{Qt 

+ E 



2^mn l^r. 



+ S{QT)sCmn 



^ dz H{Xa/z,Xb) , . 

/ , I ~ z77Z _ \ d'~^mn[Z) 



b 

4A 



dz H{Xa,Xb/z) 



xt Z H(Xa,Xb) 



bC'm,n{z^ 



r In Qt/Q] 




■ 1 ■ 


[ Qt \ 


+ 


[Qt\ 



+ 



Qt 

+ bDmn E / 



aDmn ^ ^ / 
- Jx 



1 dz H{Xa/z,Xb) . 
Z H{Xa,Xb) 



1 dz H{Xa,Xb/z) 

Z H{Xa,Xb) ''^''^ 



(C.l) 



then each Tmn is specified by writing out {2T,iT,sC,aC,iC,A,B,aD,bD} for each 
{mn}. On the next few pages, we'll lay out our findings for these parameters. 
As usual, Cp = 4/3, A'^^ = 3, and Nf is the number of participating flavors 
at the appropriate renormalization scale /i. In general, each Tmn will depend 
on the subprocess under consideration, labeled in the text by the index i, or 
more specifically by indices {a, b, d} over the types of partons participating in the 
subprocess ab 'yd, for instance qg — > 7g. For the virtual contributions, each r(j) 
is presented independently. The rest are expressed in terms of constants C(/) and 
7(/), where 
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C(g) = Cp 

C{g) = Nc 

7(g) = \c, 

7(.) ^ ^^^V^ (C.2) 



For example, C(a) + C(6) — C((i) for qq — > 7(7 is 2Cf — Nc- 

C.l r^/jri 



2-'- Virt 



^Virt 



a'^Virt 
.^(95^79) 

n (99-^79) 

(99-^79) 



-(2Cf + Nc) 
_3C _(ilNc-2Nf) 



6 



bDyirt 













-L n 



TT 



(C.3) 

(2C.-Nc)ln^ + Ncln^ (C.4) 

(C.5) 
(C.6) 
(C.7) 

(C.8) 
(C.9) 
(C.IO) 
(C.ll) 



-(2Cf + Nc)- 
o 



where 



A (99^79) 



h - In"' 

6 2s 



(llNc-2Nf), s 
— m — 



6 



3C 



In 



M' 



/ 



6 til 

2C,-^^)!^-7Cp + CpfW|^|+l„=|^ 

2/3 V s s 
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Virt 
Virt 



^ Virt W 
Virt W 



U) 
i). 



(C.15) 
(C.16) 



C.2 r-ft 



2 Tab 

Aab 
Bab 

sCab 
aCab 
iCab 
aDab 
bDab 



C(a) + C(6) - Q.{d) 





2 




C(a) 


+ C(6) 






2 




C(a) 


+ C(6) 


- C(d) 




2 




C(a) 


+ C(6) 


-C(d) 




2 




C(a) 


+ C(6) 


-C(d) 










In 



In- 



M2 



Pj Pj XaPj 6 



(C.17) 

(C.18) 

(C.19) 

(C.20) 

(C.21) 
(C.22) 
(C.23) 
(C.24) 
(C.25) 
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c.3 n^j 



C(a) + C{b) - C{d) 





2 




C(a) 


+ C(6) 


-C(d) 




2 




C(a) 


+ C(6) 


- C{d) 




2 




C(a) 


+ C(6) 


-C{d) 




2 




C(a) 


+ C(6) 


-C{d) 



In- 



2 Xbpj 



+ 










C.4 r-ft 



C(a) - C(6) + C{d) 
2 

C(a) - C(6) + Cid) 



In 



P] (2 - 



- 2(7r - aa) cot + 2yj 



C(a) - C(6) + C(d) 
2 

C(a) - C(6) + Cid) 



In ,^"^'^,1^, ''^ + 2(vr - « J cot - 2y, 
(2 - v)Q'' 

C{a) - C(6) + C{d) 



-In 



2 



In^ 



In 



(2 - 



+ 2(7r - aa) cot aa - 2yj 
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lCa2 



aDa2 
bDa2 










where v = —u/s, aa = 0^( 2(1-1,) ) '^Vj — 



Xa[l-v) 



XbV 



C.5 ri^ 



(C.40) 
(C.41) 
(C.42) 
(C.43) 



2^b2 
l^b2 



Am 

Bb2 
5Cb2 

aCb2 
bCb2 
aDb2 
bDb2 



-C(a) + C{h) + C{d) 
2 

-C(a) + G{h) + C{d) 



-C(a) + C(6) + C{d) 



2(7r — ab) cot ab — 2yj 



-C(a) + C(6) + C{d) 



-C(a) + C(6) + C{d) 



(1 — X;,) Sv 

m tt^t; — \- 2[n — ab) cot + 2yj 



(1 + tOQ' 



(C.44) 

(C.45) 
(C.46) 
(C.47) 



-In^ 



In 








Ml 



In 



(1 - XbfSv 
(1 + v)p'] 



Ml 
p] 



+ 2(7r — ab) cot ab + 2yj 



TT 



where ab = a(^) and 2y 



In 



Xa{l-v) 
XbV 



(C.48) 
(C.49) 
(C.50) 
(C.51) 
(C.52) 



C.6 rif 



2T2a = 
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(C.53) 



2a 



A2a 
B2a 

aC2a 

aD2a 
bD2a 



X 



C(a) - C{b) + C{d) 
2 

, 1 — V (tt — A(h) , . ^ , 

\n- h 2^^ — - cosA0) + 2(7r - a^) cota^ 



2-v 



sin A0 



X 





C(a) - C(6) + C{d) 
2 

, 2-v (tt - A0) , ^ ,x , X 

m z 2 — : — :— — (1 — COS A<p) — 2[n — cot aa 



1-v 



sin A0 



X 



C(a) - Cjb) + Cjd) pj_ 
2 M| 
, 2-t; (tt- A0) , 

m z 2 — : — :— — (1 — cos A0) — 2[tt — cota^ 



1-v 



sin A0 










(C.54) 
(C.55) 

(C.56) 

(C.57) 
(C.58) 
(C.59) 
(C.60) 
(C.61) 



where = and Acj) = MOD(||0g — (pj — 7r||/7r) 



C.7 rff* 



2b 



2i- 2b 







ir2fe 

A2b 
B2b 

&C2b 



-C(a) + C{h) + C(rf) 



X 



In^ h 2^^ — - cos A0) + 2(7r - at) cotaj, 



1 + w 



sin A0 







-C(a) + C(&) + C{d) 
2 

1 + (tt - A0) , A ,N / N 

m 2 — ^ — ^- — (1 — cos A<p) — 2(7r — cot 



sm 



X 



-C(a) + C{h) + C(ci) p|_ 

2 M2 

1 + (tt - A0) , A ,\ / \ 

m 2 — : — :— — (1 — cos A0) — 2(7r — cot ai, 



sin A0 
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(C.62) 

(C.63) 
(C.64) 

(C.65) 



(C.66) 



dC2b 

aD2b 
bD2b 










(C.67) 
(C.68) 
(C.69) 
(C.70) 



where ab = a(^) and A0 = MOD(||0q 



IT IT 



C.8 r 



coll 



2Ta 
iTa 
Aa 
Ba 

sCa 

dCa 
iCa 



aDa 



bDa 





7(a) + 2C(a)lni^-^ f' ^^^^§^l^l^P+(z) 



Xn 







-7(a) - 2C(a) In 



In 



2 r 



7(a) + 2C(a) In 



1 - 




1 




(C.71) 
(C.72) 
(C.73) 
(C.74) 

(C.75) 

(C.76) 
(C.77) 
(C.78) 
(C.79) 



C.9 Tf^ 



2Tb 

iTb 
Ab 
Bb 





7(6) + 2C{h) In -Y.f -^^wr^p^S^) 

Xb r J xi, Z iiyXaiXb) 







-7(6) - 2C(6) In 



I - Xb 
Xb 
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(C.80) 
(C.81) 

(C.82) 
(C.83) 
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a^b 


= 


(C.103) 


riCT 

b^b 


= 


(C.104) 




= 


(C.105) 


b-l-'b 


= 


(C.106) 



C.12 



2^2 

A2 
B2 
5C2 



aC*2 
b^2 

bD2 



C{d) 

lid) - C{d) 







A0) 



sin A0 



Mf 

(1 - cosA(/)) - In^ 



-lid) 

Cid) 
1 



2 










2^ --^(1 -cosA0) -In^ 

sinA0 Pj . 

(7r-A0), . 

^ --^ 1 - cosA0 In^ 

sm Acp Pj 

27r2 1-^^ ^2- 



13 
Y 



TX 



where A0 = MOD(||</)g 



TT /TT 



Pi 



+ 5, 



(C.107) 
(C.108) 
(C.109) 
(C.llO) 



dg- 



(17Nc - 23Nf) 



(C.lll) 
(C.112) 
(C.113) 
(C.114) 
(C.115) 



One can easily verify that the poles all cancel; that is, the sums Y^mn 2^mn 



and X^mn iTmn both equal 0. Additional cancellations will occur upon summing 
the remaining pieces, and further simplification arises upon integration over (pq. In 
particular. 



rf0,2^— — P(l-cosA0) 
sm A0 
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27rlnl6 



(C.116) 



2tt 



d(pq2{n — a{xi)) cot a{xi) 



27rln 



l + 2xi 
2xi 



(C.117) 



Thus, for example, terms involving cota( 



2(1-1,). 



will cancel corresponding ln(2 — v) 



terms, and cota(^) terms will get rid of ln(l + v) pieces. 
As it is, we have the initial-state and final-state sums: 



A- ■ 



init 



Av^rt + + A,a + A,^ + ^62 + Af' + Af' + + ^ 



A 



final 



C(a) + C(6) 
A2a + A2b + a: 



coll 



(C.118) 
(C.119) 



B 



init 



Bvirt 

C(a) 

m 

C{d) 



Bah + Bhn + Bal + B, 



5 62 



- In h In — + (vr - a„) cot a„ - (tt - ab) cot 

,Z Z — V I — V LJ^ 

1^2 — vl — v^s , , , , 

- In h In — - (tt - aa) cot + (tt - ab) cot ab 

2 1 + V V Q"^ 

1 , f (1 — v) 

- In ^2 - t>)(l + t>) + ~ "f*) + (vr - afe) cot at 



(7(a) +7(6)) 



(C.120) 



Bfr, 



final 



B2a + -82;, + B2 



coll 



= C(a) 

+ m 

+ C(d) 



(tt — aa) cot + (vr — «;,) cot ab 



1, 

- m 

.2 1 + 1 - 

l-l+t"! — t", , , , 

- m h (vr — aa) cot — (vr — ab) cot 

1, (2-v)(l + v) , , , , 

- In r (vr - aa) cot - (vr - ab) cot ab 

.2 v[l — v) 



2^ --^ 1 -cosA0 

smA0 ^ ' 



(C.121) 



init 



C(a) 



2 p] 



2 p 



In h 2 In ^ 

2 — 1 — p. 
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final 



+ 2(7r — aa) cot — 2(7r — ab) cot at 



+ 



+ C(6) 



-In^ 



M 



/ 



M 



2 ~ o In 2 



In 



2-vl-v 

1 + V V 



+ 21n^ 
Pj 



— 2{ti — a„) cot aa + 2(7r — a^) cot 



2 p] 



2 r 



In 



f (1 — v) 



+ 2(7r — aa) cot + 2{n — a^) cot 



Mi 



B 



+ (7(a)+7W)ln::^ + Ao. + ^-C(rf)^ 

To 6 



C(a) M,^ 



2 r 



2 P,^ 
0(5) ^^M^ 
2 p 



In 



1 + f 1 — f 



2 r 



In- 



V V 

V V 



1 + V 1 — V 



(C.122) 

2{n — aa) cot — 2(7r — ab) cot at 
2(7r — aa) cot eta + 2(7r — a?,) cot 



+ C{d) 



M 



-In „ 

2 P,^ 



In 



f (1 — v) 
{2-v){l + v) 



2(7r — tta) cotOa + 2(71 — ab) cotofo + 4-^-^ — 7rT^(^ ^ cos A0) 



M 



- In „ 

2 P] 



4 



- l{d) 



sm 
(tt - A0) 



(1 - cos A0) 



sin A0 

iin^^; ^ 

2 V 3 



(1 — cos A(f)) — In — ^ 

sin A0 



+ 5, 



/ . ^27r 
17Nc - 23Nf 



13\ 



(C.123) 



aCinit ~l~ bCinit 
aC final ~l~ bC final 



^ aa^ bb 



In _f p+ + p+ 



(C.124) 
(C.125) 



a Ami = 1 (C.126) 

bDimt = 1 (C.127) 

aDfinal = (C.128) 

bD final = (C.129) 



Calculation of the CSS {A, B, C} parameters from these results is outlined in Section 
6.5. 
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For the m = 1 Bremsstrahlung contributions, we derive (to order 1): 



'2(7r- A0)^ 
sin A0 



cos A</)) + 2 In 



z{l- z)pj 
Mf 



(C.130) 



where z = {pj — QT)/Pj, A0 



7r. The poles cancel. 
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APPENDIX D 



THE MONTE CARLO METHOD IN HIGH 
ENERGY PHYSICS 



D.l Basic Theory 

Monte Carlo is, at its essence, a numerical method of performing integrations. 
Compared with other methods, it has certain features which make it particularly 
well-suited to the types of problems encountered in high energy physics. The 
fundamental theorem is based upon the recognition that the average of a function 
over a domain V can be expressed as the average of a set of discrete samples of the 
integrand, in the limit that the number of samples goes to infinity [93]: 

Here Xj is a particular vector in the domain V of dimension n, and N is the number 
of samples, each generated randomly by computer. 

In practice, the generally available random number generators choose values on 
a — s> 1 scale, and the function to be integrated is then re-expressed in terms of 
variables Ui which have been thus scaled. Now (for each variable), 

/ f{x)dx becomes / f{y)dy, where (D.2) 

y = "^'"'"^ and the volume V = f[[ f' dy,] = 1" = 1 . (D.3) 
The theorem is then expressed as 
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1 

jryf{y)=rim-Y.f%). (D.4) 

Often there will be variables i* in / that are not integrated over, but are not 
fixed. They define a desired distribution / dyf{y,z). Instead of stepping through 
values of zj and calculating / dyf{y, zj) at each, we can define bins in ^ and expand 
our definition of the volume V to include the domain of z. Each evaluation fXyj, zj) 
is then added to the appropriate bin in z, and the result / dyf{y, {z}fj for bin k is 
calculated by dividing by the number of evaluations A^. 

Further integrations can be done from here. For instance, if z consists of two 
variables Zi and Z2, we have a two dimensional grid of bins in space: 



T 

Zl 



Z2 ^ 

Integration over Zi then consists of summing the contents of each column of bins into 
a new one-dimensional set of bins in Z2: 



Z2 — ^ 

While evaluating an integral, Monte Carlo can also be used to compute derivatives 
of this integral with respect to a variable. If Ax is a region in the domain V, then 
the derivative with respect to x is found by binning in x and dividing by the width 
of the bin [94]: 



195 



where A J is the sum of those weights Vf/N which fall in the range Ax. In practice, 
we find the approximate derivative as we can't take the limit. In fact, decreasing 
Ax increases the uncertainty since fewer points fall into each bin. Note that we can 
differentiate with respect to any variables, not just the ones being used to evaluate 
the integral, as long as the former are expressible in terms of the latter. 

D.2 Advantages of Monte Carlo for High Energy Physics 

Among the many quantities of interest in high energy physics, surely the one 
most frequently evaluated through numerical integration is the cross section a for 
a particular reaction. The integrand is a differential cross section ^j^-j which, 
depending on the order of the calculation, may be a function of many four-momenta 
{Pi}, and of course each four- vector has four components. Integrating over all of these 
variables would yield one number - the total cross section; but often we will want to 
leave one or more uninintegrated so as to see the distribution in these variables. Also, 
some variables may be fixed by the associated experiment (like the center-of-mass 
energy), or there may be complicated relations among the variables which fix one 
of them with respect to others. These things would all tend to reduce the number 
of integrations required, but in practice it is not always clear at the outset which 
variables will be of interest. Furthermore, some of the integrations may be impossible 
analytically, either due to the iterative nature of the integrand (as in the case of 
parton distributions) or to the complicated, discontinuous domains of integration 
that can arise when experimental constraints are applied. It is of necessity, then, to 
seek the fastest, simplest, and most flexible tool for the job. 

The fundamental theorem of Monte Carlo integration allows us to think of the 
situation in the following way: We have a space V of randomly selected points x, 
each point carrying a contribution Vf{x)/N to the total integral, and each section of 
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our program, from selection of the points, through apphcation of the cuts, evaluation 
of the weight, and finally binning, can carve up this space in its own way, with its 
own variables. The following considerations should make this advantage clear. 

Speed. For integrals over many variables, Monte Carlo converges faster than 
other methods. For example, to do an n dimensional integral with A^-point Gaussian 
Quadrature, the integrand would need to be evaluated A^" times. Using randomly 
distributed vectors, by contrast, requires only evaluations. 

Variable Transformations. Analytically, if one has calculated a distribution 
(say ^^^^^^ ) , and one now wants this as a distribution in the related set dydr, one 
must calculate the Jacobian function for the transformation: 



da da 



dy dy 

dxi 8x2 

St dr 

dxi 8x2 



(D.6) 

dydr dxidx2 ^ ' 

In Monte Carlo, one simply bins in the new variables, modifying the bin sizes if 

desired. That is, the weights a may still be evaluated using random {xi, X2} values; 

one need now only insert the calculation of {y,T}, and use these values instead to 

decide on the bin. 

Discontinuous Domains. There are times when an integral we wish to evaluate 
has a domain that is not analytically expressible, is discontinuous, or is more 
easily expressed in variables other than the ones which provide the simplest form 
for the integrand. In a typical high energy cross section, it may be that each 
Lorentz-invariant piece is simplest in a different frame, and experimental cuts may 
be most easily expressed in yet another frame or coordinate system. There may be 
regions of the detector which are blind by design or damage, and which therefore 
should be similarly removed from the theoretical calculation. Creating these types of 
"holes" in the domain would violate the principles upon which Gaussian quadrature 
and Simpson's Rule are based; in Monte Carlo the evaluation of the integrand at 
any two points is divorced from continuity constraints between these points, and the 
integral becomes possible. 
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The trick is to use a larger volume V which contains V, and simply ask, for each 
random vector x in V, whether or not it is in V. This prescription allows us to 
perform the random selections Xi within a conveniently continuous space, within one 
subroutine, and leave until later (a different subroutine) the question of whether a 
particular point x is in the desired domain V. The condition (or cut) can be described 
using any quantities calculable from x, with the result being a decision on inclusion 
of the point's weight in the integral: 

lixeV then WEIGHT = fix)V/N 

lix^V then WEIGHT = 
To minimize the uncertainty, the volume V should be as close to V as possible. 

Estimate of Precision. In contrast with Gaussian Quadrature, Monte Carlo 
integration directly produces an uncertainty estimate. For large (the number of 
sampled points), the variance is given by [94]: 



a 



\ — 



(D.7) 



1 N 1 N 

Where {w) = -Y.m), {w') = -^^fm ■ (D-8) 

1=1 1=1 

Having a numerical value for the statistical variance of an integral is a first step 
toward trying to decrease it, as we will see. 

Event Generation. Monte Carlo can be used as an event simulator in particle 
physics. The fact that a particular randomly chosen event Xi yields an evaluation 
f{xi) proportional to its probability of occurrence can be turned around - once 
we know which regions of phase space give the greatest weights, we can generate 
sets of four-vectors (events Xi) with frequency f{xi)AV [95]. Such simulators 
(e.g. PYTHIA, HERWIG, ISAJET) are helpful in determining background rates 
for unwanted processes. This "cross-referencing" of the weights to the points usually 
turns out to be necessary anyway, for reasons discussed below. 
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D.3 Importance Sampling 



Monte Carlo with uniformly chosen ordinates has one major drawback: the 
frequency with which a region AV is sampled depends only on the size of AV, 
not on the size of the integrand there. Thus just as much computing time is spent 
generating contributions in regions of low weight as in high-weight regions. Avoiding 
this inefficiency is especially important for most particle physics cross sections, which 
are rapidly falling functions of their arguments. If the function / were simple enough, 
certain changes of variable might sufficiently flatten it [96], but particle physics cross 
sections are rarely simple, and for good results some means of importance sampling 
is usually required. 

If we select our points x with a different density p{x) in each region of space, 
we will be more efficient; but if we expect to get the same answer for the integral, 
then those regions with more points should count less. In other words, normalization 
will require that the size of the volume element AV^ about point x go as l/p{x), 
and in the above formulae/ (x) is replaced with f{x)/p{x). It can be shown that the 
uncertainty is minimized when the density has the same shape as the absolute value 
of the integrand: 



but there are certain means of approximating the best density on an iterative basis. 
If one divides the integration volume into smaller volumes, then as implied above 
there are two ways of modifying the density of selected points - either have equal 
subvolumes with a different number of points in each, or have varying subvolume 
sizes each with the same number of points. The latter method, called stratified 
sampling, is used in most of the common algorithms, such as VEGAS and SHEP 
[97]. Initially, the subvolumes are of equal size. At each iteration, a two-point Monte 
Carlo integration is performed in each subvolume, generating a contribution to the 



Best 




(D.9) 



Of course, if we could calculate this, all our work would be done already; we can't. 
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total integral and to the variance. This information is then used to redefine the 
subvolume sizes for the next iteration, such that concentrations of small subvolumes 
are built up in regions where the variance was initially largest. The overall variance 
is thus decreased at each iteration until no significant change occurs. The better 
algorithms are able to use the integral information from each step, not just the last. 

D.4 Program Structure and Example 

Here we will outline the main sections of a Monte Carlo program, using hadronic 
two-photon production as our example. At Born level, a quark from one hadron 
annihilates with an antiquark from the other, to emit two photons whose momenta 
we detect. At higher orders, we may have additional particles in the final state, 
for example gluons that are radiated from the incoming quark legs. We'll stick with 
one-gluon emission here, and ask the question, "How much transverse momentum will 
this extra gluon impart to the two-photon system?" Without the gluon, the photon 
system has zero transverse momentum, as there is none coming in (the quarks are 
assumed coUinear). With the additional kick, the question makes sense, but so do 
many others, so we want to keep as many of the photon variables unintegrated as 
possible. In other words, we look at: 

da{AB ^ 27 + X) = f dxA f dxB[qf{xA)qf{xB) + {f ^ f)]dHff 27) , 

(D.IO) 

where 



dd-iff ^27) ~ J d pgd p^^d 

X + p^;- - - p^^^ - p^) ^ \M{ff 77^?)!' • (d.ii) 

Note that we are integrating only over the gluon, as we want all photon variables to 
be free. 
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Let's see what we can do to get this ready for Monte Carlo integration. We 
know Monte Carlo doesn't like delta functions, so we get rid of these first. Out of 
the 20 variables we started with (five 4-vectors, one for each of the two incoming 
quarks, the outgoing photons, and the gluon), we now are left with 16. Additionally, 
we are assuming the partons and photons to be massless, which means that one 
component of each 4-vector is fixed relative to the others. This brings our variable 
count down to 12. Fixing our l-axis along the beam line also helps, as the incoming 
partons now have no transverse momentum, which gets rid of two variables apiece. 
Finally, we work at fixed total hadronic center-of-mass energy v^, so we are left 
with seven variables. Since it takes six to specify the two photons, and these we 
wish to leave free, we will actually need only one "forced" integration (over one or 
the other momentum fraction x). Integration over the others (the photons) will be 
left for the user to decide upon through binning. Whatever is left determines the 
desired distribution. That is, our total integration volume V is 7-dimensional, as we 
add the free variables to the volume as discussed in Section D.l. For the purposes of 
generating the random points in this space, and their weights, we need to choose a 
set for which we know the kinematic bounds. For binning of these weights, the user 
may choose any set he or she wishes. 

With a bit of analytical work, then, we obtain a simplified expression, each piece 
of which we will write as a function of the simplest parameters for that piece: 

da{AB ^ 27 + X) ~ L dx[qf{x, Q')qf-{x, Q') + (/ ^ f)]P{Q',pT', y, M 

f 

^ ^ ^ ^ ^ ^ 2 
X X! 1-^(^12' ^13, ^14,^15,^23,^24,^25)1 • (D.12) 

Here x is a known function of x, given by the massless gluon constraint. Q^, pt"^, y, 
and are, respectively, the (squared) energy, transverse momentum, rapidity, and 
azimuthal angle of the two-photon system in the hadronic cm frame (which we take 



201 



to be the lab frame - at a collider). P is a function which depends purely on phase 
space. 

The matrix element for the subprocess is most easily written in terms of Man- 
delstam variables stj and Uj. The astute reader will notice that there are seven of 
these, so they cannot all be independent of the four we've already mentioned. With 
a little extra work, we could choose a frame and rewrite the matrix element in terms 
of Q^, pt'^, y, 0275 X, and two more variables, but the result would not be pretty. 
And indeed there is no need to. All that matters is that we be able to calculate the 
4-vectors of each particle in whatever frame is easiest for that particle, boost them 
all to the same frame, and calculate the Mandelstam variables. 

The best frame for the photons is the rest frame of the two-photon system, in 
which their momenta are equal and opposite, and their energies are simply Q/2 each. 
The direction of their "momentum line" in this frame can be anywhere from to vr 
in theta, and to 27r in phi. Those three variables completely specify the photon 
4-vectors in that frame, and now we're done seeking variables. We have seven, and 
know their limits: 



<p,^ <i|±^-Q' 



<Q' <S 

_2 ^ 

4S'cosh^ y 

ymin — y — Vm 

<(f)o^ < 2tt 



<9^ < vr 
< 0^ < 27r 
-^e^ <x < 1 (D.13) 

Note that in reality will never approach S (unless the hadrons disappear 
altogether!). We are taking a larger volume (as discussed in Section D.2) for lack of 
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a better limit. Here ymin and Umax are assumed to be given by the user, based on 
the viable range of the detector. We are now ready to begin writing the code. 

In the following, variables entering a subroutine will be denoted by lowercase 
names, outgoing variables will be in uppercase, and variables which are both inputs 
and outputs will be italicized. To avoid confusion, we'll show the passed variables 
only in subroutine definitions, not in calls. 

Program MAIN 

Initialize fixed values. 
Initialize bins 
Loop {i = 1...N) 

Call EVENT 

Call WEIGHT 

Call BINIT 

Write LABVECS, WEIGHT (optional) 
end Loop 

Here we read our fixed values (e.g. S, ymin, Umax, etc.), make sure our bins 
are empty, and then process N randomly selected events. For each event, we find 
the four- vectors (EVENT), calculate this event's contribution (WEIGHT), put the 
weight in the appropriate bin (BINIT), and optionally, send the event out to a file. 

Subroutine EVENT (fixed, LABVECS [5,0:3], UNSCALED, VOLUME) 
Random selection of scaled vars (e.g. XQ2=RAND (SEED) ,XPT=RAND (SEED), etc.) 
Scaled vars ^ Unsealed vars (e.g. Q2=S*XQ2) & VOLUME 
Find photon cm 4-vecs (using Q^, 9^, 0^) 
Find boost parameters (using Q^, pr^, y, 4'2'y) 
Boost photons to lab frame 
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Find parton 4-vecs in lab (from S, x, x) 
Find gluon 4-vec {= p'j + p'j — pi^i — 
return 

Remember that we can't randomly choose (what I call) the "unsealed" variables 
[Q"^ 1 Pt^i y, 027, df, 07, x] directly; our random number generator picks between the 
limits and 1. So we define a "scaled" variable (e.g. XQ2 in the code) for each, choose 
it randomly, and then calculate the unsealed variables from them. The volume V is 
the product of the unsealed variables' domains, and this will be multiplied by the 
integrand in evaluation of the weight for this event. Once we find the four- vectors 
for all the particles of the subprocess, we'll be able to calculate the matrix element; 
we also send out the unsealed variables, as they may be more directly useful for the 
phase space factor, cuts, and binning. 

Subroutine WEIGHT (fixed, labvecs, unsealed, volume, WEIGHT) 

Implement CUTS (using labvecs, unsealed), e.g. IF(P5(0) .LT.0.75) GOTO 
REJECT 

Calculate MANDELSTAM vars tij (using labvecs) 

Calculate MATRIX ELEMENT (using Uj) 

Calculate PARTON DISTRIBUTIONS qfix,Q^) 

Calculate PHASE SPACE FACTOR P (using unsealed) 

INTEGRAND = (partons) x (phase space) x (matrix element) 

WEIGHT = (volume) X (integrand) /(iV) 

return 

REJECT: 

WEIGHT = 

return 
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The first section here (CUTS) is open to user modification. The user has access to 
all of the variables that describe the event, and may invent new variables as needed, 
all in order to define IF statements which either accept or deny this event as valid. If 
denied, control is passed to the REJECT line, where the weight is set to zero. This 
CUTS section may be put in the binning routine instead, but this just wastes time, 
as the weight will be unnecessarily calculated. 

Subroutine BINIT (fixed, labvecs, unsealed, 1:100], WEIGHT) 
Decide upon and calculate binning variable(s) (e.g. the pair PT in our case) 
BIN# = INT((PT-LO)/STEP) + 1 
BIN(BIN#)=BIN(BIN#) + WEIGHT/STEP 
return 

Again, the first section here is user-accessible, for the purpose of defining whatever 
variables are of interest. The user may wish to look at the transverse momentum of 
the photon pair (as we do), rapidity differences between the photons, or single photon 
spectra (in which case each photon might contribute to a different bin). Here we 
assume a one-dimensional bin array, corresponding to but the generalization to 
higher dimensions is trivial, and affects no other subroutine. Note that our WEIGHT 
corresponds to da, so if we want we must divide by the bin width (STEP) - see 
the discussion on differentiation in Section D.l. 

In figure D.l we present results for two-photon production in Ti~p collisions, to 
first order in a^. The fixed quantities are: S'=526.33 GeV^, Aq(7£)=0.2305 GeV, 
STEP=0.08 GeV. The cuts that were used are the following: 

1. The photon of higher px must have px > S.OOGeV. 

2. The photon of lower pt must have px > 2.75GeV. 

3. The quantity zz = —{pE^ ■ pr^l ) / max (^^^,1 , PTt,2 ) must exceed 2.75. 

205 



and the code that took care of these cuts: 

C — CUTS ON LAB VARS — 

PT3=DSQRT (P3 ( 1) **2 . D0+P3 (2) **2 . DO) 

PT4=DSQRT (P4 ( 1 ) **2 . D0+P4 (2) **2 . DO) 

ZZ=-1 . DO* (P3 ( 1 ) *P4 ( 1 ) +P3 (2) *P4 (2) ) /DMAXl (PT3 , PT4) 

IF (DMAX1(PT3,PT4) .LT.3.D0) GOTO REJECT 

IF (DMIN1(PT3,PT4) .LT.2.75D0) GOTO REJECT 

IF (ZZ.LE.2.75D0) GOTO REJECT 

On the graph, the pt distribution of the photon pair is shown. The dashed hne is the 
perturbative 1-gluon result; the sohd hne is the result of resumming the soft gluons. 
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Figure D.l. vr p ^ 77X pair px distribution. 
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